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Abstract 



Frustration refers to the inability of a system to simultaneously satisfy all interactions. 
The principal ingredients of frustration are (i) competition between interactions within the 
system, and (ii) conflicting constraints established by symmetry, or higher energy interac- 
tions. The behavior of frustrated systems is generally difficult to predict, and often results 
in novel cooperative states and phenomena. Frustration is common in disordered systems; 
in this dissertation, however, we focus on frustrated systems without quenched disorder. 
Effects of frustration tend to be pronounced when coupled with weak connectivity. Weak 
connectivity refers to a topology in which the order in one part of the system does not con- 
strain the order in other parts. We have used neutron scattering to study structurally and 
magnetically frustrated systems, both strongly and weakly connected. 

In ZrW 2 8 and Sc 2 (W0 4 ) 3 , anomalously large Negative Thermal Expansion (NTE) 
has been observed over a wide temperature range. Analysis indicates that the NTE effect 
in these materials is driven by folding of rigid units consisting of strongly bonded oxygen- 
metal polyhedra. Weak connectivity is manifested by structural under-constraint, which 
enables soft optical phonons. Using neutron scattering, we observed an unusually high 
density-of-states at low energies, and mapped the dispersion relations of the most intense 
phonon modes. While in conventional materials soft phonon modes typically signal an 
imminent structural phase transition, in NTE materials the transition is frustrated by the 
incompatibility between the crystal symmetry and the symmetry of the soft phonon modes. 
This results in the unusual NTE effect. 

RbFe(Mo0 4 ) 2 is a frustrated, quasi-classical (5 = 5/ 2), quasi-two-dimensional Heisen- 
berg antiferromagnet on a triangular lattice, with easy-plane anisotropy. In RbFe(Mo0 4 ) 2 , 
frustrated units (triangles) share edges, and are strongly connected. Although a long range 
magnetic order is established below T N = 3.5 K, approximately 1/4 of the magnetic mo- 
ment remains in a fluctuating state, even at the lowest temperature that we accessed (1 .7 K). 
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Neutron scattering measurements indicate that the ordered state consists of a so-called 120° 
planar structure. Weak coupling between the planes causes an incommensurate ordering 
along the c-axis. Application of an in-plane magnetic field results in a sequence of phase 
transitions, and we mapped a rich H-T phase diagram. Excitations consist of spin waves, 
and we measured the spin wave dispersion relations, both in-plane and along the c-axis. 

QS-ferrite (Ba 2 Sn2Ga3ZnCr 7 022) is a frustrated, quasi-two-dimensional, S = 3/2, 
antiferromagnet on a kagome sandwich lattice. Frustrated units (tetrahedra) are linked 
through corners only, and the system can therefore be characterized as weakly connected. 
QS-ferrite does not have a magnetically ordered structure at any temperature (> 0.1 K), 
despite the strong magnetic interactions (6cw = —312 K). We use inelastic neutron scat- 
tering to study spin dynamics in the system, and observe low temperature spin freezing, 
with a timescale dependent onset and magnitude. The relaxation rate distribution displays 
uj/T scaling, indicating proximity to a quantum critical regime, and anomalously broad- 
ens below 1 1 K. The wave- vector dependence of the excitation spectrum is consistent with 
the notion that the low temperature dynamics is driven by residual interactions between 
satisfied spin simplexes. 

Work in this dissertation was performed in the Department of Physics and Astronomy 
at the Johns Hopkins University, under supervision of Professor Collin Broholm. 
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Preface 



Mirum enim, quam belle ab ipsa consideratione Naturce ad ccelestium rerum contem- 
plationem disponitur animus, & ad ipsum Divinum tantai molis Conditorem assurgit, 
infinitam ejus Potentiam, Sapientiam, Providentiam admiratus, que erumpunt undique, 
& ubique se produnt. 

— P. Rogerio Josepho Boscovich, S.J., Theoria Philosophise Naturalis (MDCCLXIII). 1 

Science is a search for Beauty, be it in "fire, or wind, or the swift air, or the circle of 
the stars, or the violent water, or the lights of heaven." And the Beauty lurks at unexpected 
corners, and makes itself seen only by the patient and humble. And they search for Her 
diligently, "because the things are beautiful that are seen." 

Unfortunately, science is today rarely seen as a search for Beauty. Frequently, it is a 
blind conquest of knowledge. Ever since the sun became the center of our universe, we 
became increasingly marginalized in our own world. The often devastating consequences 
of this approach have been witnessed over the past few centuries. 

It is increasingly unpopular to fuse the scientific work with ethics, the science of the 
ideal of human character. Ethics of a modern scientist often applies only to his conduct in 
the lab, not to the end result of his work. It is increasingly unpopular to call wrong wrong. 

It is a responsibility of those fortunate enough to admire nature's beauty and have her 
secrets revealed, to offer humanity not only fruits and tools, but also the purpose of their 
work. Every scientific work has a potential for misuse. Indeed, even a silence can be 
misused, and how much more so a work with no strings attached, without a guide for the 
user. 

It is my desire that the work presented in this dissertation be used in a benevolent 
manner for all humanity, and never to harm anyone. 



1 For it is marvelous how exceedingly prone the mind becomes to pass from contemplation of Nature herself 
to the contemplation of celestial things, & to give honour to the Divine Founder of such a mighty structure, 
lost in astonishment at His infinite Power & Wisdom & Providence, which break forth & disclose themselves 
in all directions & in all things. — Josip Ruder Boskovic, A Theory of Natural Philosophy, dedication (1758). 
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Chapter 1 
Introduction 



The term "frustration," as used in this dissertation, refers to the inability of a system 
to simultaneously satisfy all the local interactions. Frustration is a universal phenomenon, 
manifesting itself not only in physical systems, but also in technological, biological, eco- 
logical, economic, and social systems. Indeed, nearly any real system of sufficient com- 
plexity will contain a degree of frustration. The principal ingredients of frustration are 

(i) competition between interactions within the system, and (ii) conflicting constraints es- 
tablished by symmetry, or higher energy interactions. This categorization is somewhat 
arbitrary, but can be justified by the following examples. The triangular antiferromagnet is 
a frustrated system because of (i) the competing local antiferromagnetic interactions, and 

(ii) the symmetry imposed by the triangular lattice. If either ingredient is removed, the 
frustration is relieved. The triangular ferromagnet removes competition between interac- 
tions, and is not a frustrated system. On the other hand, a phase transition in which some 
interactions become favored over others (as in ZnCr 2 4 at T = 12.5 K), eliminates the lat- 
tice geometry needed for frustration. One can find similar examples in many other systems; 
in socio-economic systems, for instance, (i) competing personal interests, and (ii) imposed 
socio-economic structure. 

For a system to be frustrated, the energy scale which fixes symmetry must be much 
larger than that of local interactions; otherwise, the frustration can be relieved by altering 
the symmetry. The behavior of a frustrated system is generally difficult to predict, and 
often leads to novel states and phenomena, characterized by composite degrees of freedom 
(CDF). Although omnipresent, frustration is analyzed only in the context of physical solid 
state systems in this thesis. 

Composite degrees of freedom are common in many-body systems. They are usually 
characterized by an energy scale far below that associated with the constituent elements. 
Examples of CDF can be found at all scales in the universe, from nuclei which consist 
of protons and neutrons, to galactic clusters made of galaxies. The simplest and perhaps 
most common CDF in solid state systems is a phonon: an excitation involving many atoms. 
While the energy scale characterizing interatomic bonds is in the eV range, phonon energies 
lie in the meV range. CDFs need not be extended in space, as is the case with phonons; in 
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frustrated magnets, for instance, they can be local. 

We shall start this introduction by briefly describing un-frustrated systems, and then 
contrast them to frustrated systems. A local interaction is satisfied when the degrees of 
freedom involved are in their minimum energy configuration. The ability of a system to 
simultaneously satisfy all the local interactions, generally yields a ground state exhibiting 
long range order, and a small (or zero) ground state entropy. Examples of un-frustrated sys- 
tems are numerous, particularly in solid state physics. Nearly all ordered, periodic struc- 
tures, such as perfect crystals, can be described as un-frustrated. Conventional magnets, 
such as those in which ferromagnetic (FM) or antiferromagnetic (AFM) long range or- 
der minimizes the energy of every pairwise interaction between spins (i.e. constitutes the 
ground state), are not frustrated. Un-frustrated systems are generally easier to understand 
and describe theoretically than frustrated systems. 

Magnetism presents an ideal setting for studying frustration. The unique advantages 
of magnetic systems include (i) choice of size and dimensionality of individual degrees of 
freedom (spins), (ii) choice of size, dimensionality, and topology of the underlying struc- 
ture (magnetic lattice), (Hi) variability of spin-spin interactions (type, range, (an)isotropy), 
and (iv) suitability for experimental study. Additionally, an order parameter is generally 
(though not always) easy to construct in magnetic systems; for instance, the staggered mag- 
netization in antiferromagnets. Disorder can be introduced, and fine-tuned via doping. An 
applied magnetic field provides an external parameter to which spins couple directly. For 
these reasons, magnetic systems have been central to scientific work on frustration, and 
most of the examples used in this section will be based on magnetic systems. 

Frustration is common in systems with quenched disorder. Perhaps the best examples 
of frustrated, disordered systems are found among spin glasses, where frustration is caused 
by randomness in (i) the position of spins, and/or (ii) the sign and the strength of mag- 
netic interactions. Spin glasses form a class distinct from conventional magnets. They 
establish no long range magnetic order at any temperature. However, there can be a well- 
defined transition temperature, T g , below which the spins freeze into a highly irreversible, 
meta-stable state. Spin glasses have a highly degenerate, but discontinuous, ground state. 
(Discontinuity here refers to the existence of energy barriers separating different states.) 
Low temperature hysteresis of the response function (i.e. the magnetic susceptibility) is 
another distinctive feature of spin glasses. Some of these properties, such as a highly de- 
generate ground state, arise from frustration. Discontinuity in the ground state, on the other 
hand, probably arises from randomness. (Although this is still an open question, some re- 
search suggests that this is indeed the case. For instance, the ground state of a frustrated, 
but structurally ordered, antiferromagnet on a pyrochlore lattice has a highly degenerate, 
but continuous, ground state [1].) The frustration in spin glasses is caused by, and therefore 
inseparable from, randomness. It would, however, be interesting to separate the two effects, 
and study the effects of frustration alone. Luckily, there is a class of magnets, referred to as 
geometrically frustrated magnets (GFM), in which frustration does come without structural 
disorder. 

Geometrically frustrated systems combine competing interactions with a perfectly or- 
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Figure 1.1: (a) Un-frustrated antiferromagnet (AFM) consisting of Ising spins on a square 
lattice. Ground state is the Neel state (shown), (b) If one of the antiferromagnetic inter- 
actions is replaced by a ferromagnetic (FM), the system becomes frustrated, as the lower- 
right spin cannot simultaneously minimize both of its interactions. If FM interactions are 
randomly introduced, the system becomes spin glass, (c) A simple frustrated system, con- 
sisting of Ising spins on vertices of a triangle, with AFM interactions. The top spin cannot 
simultaneously minimize both of its interactions, (d) A frustrated system consisting of 
Ising spins on vertices of a tetrahedron, with AFM interactions. Two top spins cannot si- 
multaneously minimize all of their interactions. Triangles and tetrahedra are the building 
blocks (frustrated units) of geometrically frustrated magnets. 
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Figure 1.2: Geometrically frustrated magnets, (a) Triangular lattice, (b) Kagome lattice, 
(c) Face Centered Cubic (FCC) lattice, (d) Pyrochlore lattice. Interactions between the 
spins are antiferromagnetic. 
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Figure 1.3: Classical ground state for antiferromagnetic Heisenberg spins on vertices of (a) 
a triangle, and (b) a tetrahedron. The ground state corresponds to zero total spin, 5 = 0, 
and is three-fold degenerate for a triangle, and five-fold degenerate for a tetrahedron. 



dered underlying structure. A spin chain, for instance, in which both the first and the second 
neighbor interactions are antiferromagnetic, and of equal strength, is a geometrically frus- 
trated system. This model is, however, hard to realize in a real material for experiments. 
Instead, most real frustrated systems are based on so-called frustrated units, such as trian- 
gles and tetrahedra. It is easy to see that the antiferromagnetic interactions between spins, 
located at the vertices of a triangle or a tetrahedron, cannot all be simultaneously satisfied 
to the fullest extent. Frustrated networks can hence be built from such frustrated units. 
(Ref. [2, 3] provide a nice and thorough review of geometrically frustrated magnets.) 

Geometrically frustrated systems can generally be divided into two broad categories, 
according to connectivity; namely, into (i) weakly connected, and (n) strongly connected 
systems. Weak connectivity refers to a lattice topology in which order in one part of the 
system does not constrain order in other parts of the system. Strong connectivity refers 
to the opposite case. Weakly connected systems typically consist of frustrated units that 
are connected at corners. Representative of a weakly connected frustrated system in two 
dimensions is a kagome lattice, which consists of a network of corner sharing triangles, 
while in three dimension the standard example is a network of corner sharing tetrahedra, 
as in the pyrochlore lattice. On the other hand, the two-dimensional triangular lattice, 
consisting of edge sharing triangles, is a strongly connected system. In three dimensions 
strongly connected system would be an FCC lattice, a network of edge sharing tetrahedra. 

In this thesis we study structurally and magnetically frustrated systems. In chapter 4 we 
present our research on Negative Thermal Expansion (NTE) materials zirconium tungstate 
(ZrW 2 8 ), and scandium tungstate (Sc 2 (W04) 3 ). In these materials, soft phonon modes, 
that usually signal the proximity of a structural phase transition, are frustrated by the lattice 
symmetry. The result is an anomalously large negative thermal expansion over a very wide 
temperature range. In NTE materials, weak connectivity is manifested by the small number 
of links (i.e. low coordination number) between rigid units, which participate in lattice 
dynamics. We refer to this feature as the openness of the structure. Our research indicates 
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> Dynamics driven by 
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Table 1.1: Properties of conventional and quantum magnets. 



that both the frustration, and the weak connectivity are important for the anomalous NTE 
effect in these materials. 

In the remainder of the thesis, we study geometrically frustrated magnets. In strongly 
frustrated magnets, quantum fluctuations are significant even at T = 0. Therefore, they are 
good candidates for studies of quantum critical (QC) phase transitions. Magnetic systems 
in which spin fluctuations are significant at T = are referred to as quantum magnets, as 
opposed to conventional magnets in which spin fluctuations in the ground state are neg- 
ligible. Table 1.1 lists some of the properties distinguishing quantum from conventional 
magnets. A quantum critical transition is a phase transition that occurs at T = 0, as the 
variables characterizing the system (spin, dilution) or its environment (pressure, magnetic 
field) are varied. Chapter 6 gives a brief introduction into quantum phase transitions. We 
have already mentioned a large variability of properties that make magnetic systems suit- 
able for the study of frustration. The same properties make them ideal for studying quantum 
phase transitions. Table 1.2 arranges the three magnetic materials we studied according to 
dimensionality and connectivity. Figure 1 .4 shows a probable location of these materials in 
the quantum critical phase diagram. 

Chapter 5 presents our research on rubidium iron molybdate (RbFe(Mo04) 2 ), a frus- 
trated antiferromagnet on a triangular lattice. The frustrated units (triangles) share edges in 
a triangular lattice, and are strongly connected. Therefore, RbFe(Mo0 4 ) 2 is an example of 
a frustrated system without weak connectivity. Although there are still large spin fluctua- 
tions at low temperatures, RbFe(Mo0 4 ) 2 does establish long range magnetic order, and has 
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Figure 1.4: Quantum critical phase diagram, with probable placement of geometrically 
frustrated magnetic systems studied in this dissertation. 
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Table 1.2: Classification of frustrated magnets studied in this dissertation. 



excitations in the form of spin waves. Accordingly, this system is relatively far from the 
quantum critical point. This indicates that the effects of frustration are diminished when 
the connectivity is sufficiently strong. 

Closer to the quantum critical point, zinc chromite (ZnCr 2 4 ) is a three-dimensional 
frustrated antiferromagnet based on a pyrochlore lattice. In ZnCr 2 4 , frustrated units (tetra- 
hedra) are linked only through corners, and the system can be described as frustrated and 
weakly connected. Our research indicates that the low temperature dynamics in this, and 
perhaps in other frustrated and weakly connected systems, may be driven by the composite 
degrees of freedom (CDF). Zinc chromite has a structural phase transition at T = 12.5 K, 
which relieves the frustration [4]. This prevents the study of the low temperature quantum 
dynamics. Our research on ZnCr 2 04 is reported in Ref. [5]. 

In chapter 6, we discuss our research on QS-ferrite (Ba2Sn 2 Ga 3 ZnCr7022X another 
frustrated and weakly connected magnet. The magnetic structure of QS-ferrite is quasi- 
two-dimensional, and based on a kagome sandwich lattice (derived from the pyrochlore 
lattice). QS-ferrite is interesting as it remains frustrated (has no structural phase transi- 
tion) down to at least 1 .7 K. This enables analysis of dynamics of a frustrated and weakly 
connected system to much lower temperatures. 
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Chapter 2 

Experimental Techniques 



2.1 Introduction 

The experimental techniques employed in our research can be generally divided into 
two main categories: (i) sample preparation, and (ii) sample property measurements. Al- 
though the ultimate purpose of every experiment is a measurement of some physical prop- 
erty, only good samples will yield good results. Indeed, often it is the quality of the sample 
that determines the extent of physics that can be studied. Therefore, sample preparation 
constitutes a major part of experimental condensed matter physics. The measurements 
techniques employed for research presented in this thesis can be divided into (a) thermo- 
dynamic, and (b) spectral measurements. Thermodynamic probes enable study by mea- 
surement of statistical quantities, such as heat capacity, susceptibility, or transport proper- 
ties. Among thermodynamic probes, we used a SQUID magnetometer for magnetic sus- 
ceptibility measurements. Spectral measurements probe physical systems in reciprocal, 
momentum-energy space, by measuring correlation functions (static or dynamic, structural 
or magnetic). Among these, we used X-ray and neutron scattering. 

In this chapter, we briefly introduce these experimental techniques, with the exception 
of neutron scattering. Research presented in this thesis relied so heavily on neutron scat- 
tering, that the technique merits its own chapter (as well as a reference in the title of the 
thesis). 

2.2 Sample preparation 

As remarked in the introductory section, the quality of experimental work can be only 
as good as the samples on which it is done. While we have not synthesized the samples 
ourselves, we were very fortunate to receive excellent samples from our collaborators; it is 
a pleasure to here acknowledge their work and thank them. For each of the samples, we 
provide a reference which details the sample preparation technique. Powder samples and 
single crystals of scandium tungstate (Sc 2 (W0 4 ) 3 ), as well as those of zirconium tungstate 
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(ZrW 2 8 ), came from the group of L. Schneemeyer, G. Kowach, and A. P. Ramirez of Bell 
Laboratories, Lucent Technologies; see Ref. [6]. Zinc chromite (ZnCr 2 4 ) single crys- 
tal samples were synthesized by W. Ratcliff and S.-W. Cheong from Rutgers University; 
see Ref. [7]. Rubidium iron molybdate (RbFe(Mo0 4 ) 2 ), numerous single crystals and a 
ceramic powder sample, were synthesized by A. I. Smirnov, L. E. Svistov, and L. A. Pro- 
zorova, from the P. L. Kapitza Institute for Physical Problems, Moscow; see Ref. [8]. A 
QS-ferrite (Ba2Sn2Ga 3 ZnCr 7 022) powder sample came from the lab of R. J. Cava, from 
Princeton University; see Ref. [9]. 

2.3 X-ray diffraction 

X-ray scattering is a spectral probe particularly well suited for static measurements. 
X-rays are photons, quanta of electromagnetic radiation, with wavelength between A = 
1CT 9 m and 1CT 12 m. The technique has many similarities to neutron scattering. The differ- 
ences in use can be traced to the unique properties of X-rays and neutrons. While a neutron 
with wavevector of the order of interatomic distances in solids, k = 2n/d ~ 1 A -1 , has 
energy comparable to those of phonons and magnetic excitations, E = h 2 k 2 /2m ~ 2 meV, 
an X-ray with the same wavevector has an energy of E = hck ~ 2 x 10 6 meV, far higher 
than typical energies of collective excitations in solids. Accordingly, X-ray scattering is a 
better probe of structural (static) than dynamic properties of condensed matter. Two great- 
est advantages of X-rays over neutrons are (i) larger intensity of the X-ray beam, and (it) 
portability of X-ray sources. The two advantages, however, do not come in a "single pack- 
age:" portability comes with the tube sources, while the intense X-ray beams require use 
of dedicated synchrotron facilities. For our research, we used both tube and synchrotron 
sources. 

Tube sources 

We used the local X-ray source primarily for orientation of single crystals: princi- 
pal axes determination and alignment. We used a Philips X'Pert MRD Diffractometer, a 
four-circle (and three translation axes) diffractometer, with ceramic X-ray tube and copper 
target (Cu Kc^: X = 1.5405 A). This instrument's angular resolution is 0.003°. Although 
optimized for thin film scattering, we successfully used it, and found it invaluable, for 
alignment of single crystals and sample preparation. 

Another X-ray instrument, a Brukner SMART system, with Molybdenum target (Mo 
Kai. X = 0.7093 A), was used for structure refinement of RbFe(Mo0 4 ) 2 . We have not 
used this system directly, but through external collaboration; the details and the results of 
these measurements are given in Appendix B.l (pp. 139 et seq.). 
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Synchrotron X-ray sources 

We used the National Synchrotron Light Source (NSLS) at the Brookhaven National 
Laboratory (BNL) for our single crystal X-ray experiments on zinc chromite (ZnCr 2 4 ). 
The NSLS has two storage rings, the larger ring at energy 2.8 GeV is for production 
of X-rays, and the smaller ring at 800 MeV for UV radiation. The X-ray ring has ap- 
proximately 60 beamlines, specialized for crystallography/diffraction, spectroscopy, imag- 
ing/microspectroscopy, nanoscience characterization, etc. We used beamlines X22A and 
X22C, which are optimized for single crystal X-ray diffraction. 

2.4 Magnetic susceptibility 

A Superconducting Quantum Interference Device (SQUID) is a device capable of mea- 
suring minute changes in magnetic flux. A conventional magnetometer consists of an in- 
ductive (and conductive) coil, at which ends a potential difference is generated, and then 
amplified by standard electronic means. In contrast, a DC SQUID consists of a supercon- 
ducting ring with two (parallel) Josephson junctions. Upon entering the ring, the super- 
conducting wavefunction of the electron pairs carrying the current passes through each of 
the junctions, and acquires a phase shift which results in quantum interference. The crit- 
ical supercurrent is the maximum current that passes through the loop without generating 
a difference in potential. The interference effect causes a modulation in the supercurrent, 
dependent on the magnitude of the magnetic field passing through the loop. The modu- 
lation occurs because the magnetic field changes the phase of the wavefunction, and the 
resulting supercurrent oscillates between zero and the maximum value, the period being 
the flux quantum, $ = nh/e = 2.06 x 10~ 15 Wb. (By measuring dV/d§, the fractional 
changes in magnetic flux can also be measured.) During the measurement, the sample is 
passed through a small pick-up coil (typical diameter of a pick-up coil is about 1 cm), 
and the response measured via changes in the supercurrent. These are fitted to extract the 
magnetization of the sample. 

We used a Quantum Design MPMS SQUID Magnetometer to measure magnetic sus- 
ceptibility of RbFe(Mo0 4 ) 2 single crystals. The magnetometer has a n H = 5 T super- 
conducting magnet, an accessible temperature range T = 1.7 — 400 K, and a magnetic 
moment sensitivity m = 1 ^emu= 10~ 9 Am 2 . The results of these measurements are pre- 
sented in Sec. 5.3 (p. 73). The magnetometer is a part of the Materials Research Science 
and Engineering Center (MRSEC) at the Johns Hopkins University. 
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Chapter 3 
Neutron Scattering 



3.1 Introduction 

Most of the experimental information presented in this thesis has been, and can only 
be, acquired using various neutron scattering techniques. Today an indispensable probe of 
condensed matter systems, the neutron was discovered in 1932 by James Chadwick [10]. 
Two years later, Enrico Fermi, believed to be the greatest expert on neutrons of his time, 
discovered that slow neutrons interact strongly with almost all nuclei [11], thus indicating 
that they could be used for studying condensed matter. Yet, in contrast to the closely related 
technique of X-ray diffraction, which by then had already been used for decades as a probe 
of atomic scale structures, the development of neutron scattering was hampered by the very 
low intensity of the early neutron sources. The situation changed dramatically after Fermi, 
following the discovery of nuclear fission by Otto Hahn in 1938 [12], demonstrated the 
first self-sustaining nuclear chain reaction in 1942. Soon after, a first nuclear reactor was 
built in Oak Ridge, Tennessee. Using neutrons from this reactor, a 3.5 MW "Clinton Pile" 
(CP-1) [13], as it was called, and a two-axis spectrometer they built, Ernest Wollan and 
Clifford Shull performed the first extensive neutron diffraction experiments [14]; and so 
the neutron scattering technique was born. 

The next big impetus came with the development of the triple-axis spectrometer by 
Bertram Brockhouse in Chalk River, Canada (1956) [15]; which opened yet another dimen- 
sion, energy, to neutron scattering. The following decades witnessed the development of 
new and ingenious instruments, such as the time-of-flight, neutron diffractometer and spec- 
trometer, spin-echo, various focusing techniques, and polarized neutron scattering. These 
advances were enabled by, and also encouraged improvements in, instrument parts (de- 
tector, monochromator/analyzer). Though the neutron sources remained relatively few in 
numbers, the ever increasing flux of these sources — from 10 7 cm _2 s _1 of the CP-1 to over 
10 15 cm _2 s _1 in modern reactors, and even higher pulsed flux in spallation neutron sources 
— contributed to a growing use of the technique in condensed matter physics, as well as its 
spread into other fields of science and technology. 



14 



Mass 1.674 927 16(13) x 1(T 27 kg 
1.008 664 915 78(55) u 
1.001 378 418 87(58) m prot0 n 
1 838.683 6550(40) m e i ectron 
939.565 330(38) MeV/c 2 
Charge 



Spin 



Magnetic dipole moment -0.966 236 40(23) x 10 _26 JT- 

-1.913 042 72(45) /z N 
Lifetime 885.9(9) s 



Table 3.1: Basic properties of neutron [16]. 



The properties that make neutrons an ideal probe of condensed matter can be summa- 
rized as follows: 

(i) The de Broglie wavelength of thermal neutrons is comparable to the interatomic 
distances in solids and liquids. 

(ii) The energies of thermal and cold neutrons are of the same order of magnitude as 
many collective excitations in condensed matter. 

(in) Because they carry no charge, neutrons can penetrate deeply into the bulk of materi- 
als. 

(iv) Neutrons interact strongly with atomic nuclei via short range nuclear forces. 

(v) The neutron has a magnetic moment, and interacts with unpaired electrons in mag- 
netic atoms via electromagnetic dipole-dipole interaction. 

No other particle has all of these desirable properties. Protons carry charge and so can- 
not interact but with the surface of a material. Electrons are also limited only to surfaces. 
X-rays, though they penetrate deeper (not as deep as neutrons), are still very sensitive to 
the shape of the surface. 1 Furthermore, X-rays and electrons can only have wavelength or 
energy comparable to those of most condensed matter excitations, but not to both simulta- 
neously. Neutrons "see" light atoms even in the presence of heavier ones; specifically, light 
hydrogen is nearly transparent to X-rays, but interacts strongly with neutrons. Since they 
interact with nuclei, neutrons can distinguish atoms with similar atomic numbers, and even 

J Here, I would like to mention yet another probe, perhaps not as widely known, that I became familiar 
with during my past research: the helium atom. Used in surface physics, beam of helium atoms is ideal for 
studying interfaces, surface phonons and other related phenomena. [17, 18] 
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different isotopes of the same atom, which is not easily done with other probes. Finally, 
neutrons interact relatively weakly with matter, making them a non-destructive probe of 
physical systems. 

The major disadvantage of neutron scattering techniques is of a practical nature: Pro- 
ducing a bright (high flux) beam of neutrons requires specialized, large scale facilities, as 
the thermal energy released per neutron is extremely high (190 MeV/neutron for fission, or 
30 MeV/neutron for spallation). 



3.2 Neutron sources 

Free neutrons are commonly produced in fission reactors, or more recently through nu- 
clear spallation by energetic protons. In the later case, a beam of highly energetic protons is 
incident on a heavy target, and various particles are produced in the collision, among them 
neutrons which then propagate to neutron scattering instruments. Spallation sources create 
neutrons in pulses (at a frequency of 10-50 Hz, with about 10 14 neutrons in a single burst of 
about 1/is). Clever designs are implemented to use the time structure of this pulsed source, 
which has peak neutron intensity higher than that reactor sources. Since none of the data 
in this thesis have been obtained using a spallation neutron source, I will presently end its 
description, adding only that spallation sources have a great, perhaps dominant, future for 
neutron scattering, due to some unique advantages unparalleled by reactors. 

Nonetheless, the reactor sources, with their high continuous neutron flux (~ 10 15 
cm~ 2 s _1 ) will not become obsolete anytime soon. Using 235 U as a fuel, fission reactors 
produce on average 2.5 neutrons per reaction, such as: 

235 U + n — > 144 Ba + 90 Kr + 2n + 170 MeV, or 
235 U + n — > 94 Zr + 90 La + 3n + 197 MeV. 

The NCNR reactor, running at 20 MW power, would thus produce about 10 18 neutrons 
each second. However, only a fraction of these neutrons can be used for experiments, 
the majority being absorbed by the shielding, control rods, or in maintaining the chain 
reaction. The peak flux of the NCNR reactor is 4 x 10 14 cm~ 2 s _1 . Since the energy 
released is so much greater than the desired energy, these neutrons must be "cooled down" 
in a moderator consisting of deuterated water (D 2 0), maintained at some desired, low, 
temperature. Thermal neutrons are cooled by the heavy water that cools the reactor itself, 
and is maintained slightly above room temperature. There, in a series of collisions, neutrons 
lose energy until they are in thermal equilibrium with the moderator, and have a Maxwellian 
distribution of speeds: 

peaked at about 35 meV For lower energies, a cold source is used, and at NCNR it con- 
sists of liquid hydrogen held at a temperature of 25 K, thus giving a peak neutron intensity 
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E (meV) 


A (A) 


k (A" 1 ) 


v (m/s) 


/ (THz) 


Cold neutrons 


1.000 


9.044 


0.6947 


437 


0.242 




5.000 


4.045 


1.5534 


978 


1.210 


Thermal neutrons 


25.00 


1.809 


3.473 


2187 


6.045 




50.00 


1.2791 


4.912 


3093 


12.090 



Table 3.2: Typical values for cold and thermal neutrons. 



2 

E[meY] = = 3.956 v 2 [km/s] 

2 

= hu) = 0.6554 a; [K^V 1 
2n 2 H 2 2 



si.si/a 2 [a; 



mX 2 

h 2 k 2 , , r = -i 

= 2.072 k A 

2m L 

k B T = 0.08617 T [K] 



Table 3.3: Some often used neutron relations. 



around 5 meV. From there on, neutrons are guided through special neutron guides, tubes 
coated with neutron mirrors (Ni, for instance), using total internal reflection to minimize 
the loss of neutrons. Thence, neutrons are delivered to individual instruments. 

Thermal neutron are commonly considered those with temperatures between 60 and 
1000 K, while those with energies above or below are referred to as hot or cold, respec- 
tively, the boundary not being precisely defined. Thermal neutrons at room temperature, 
T = 293 K, have average kinetic energy E = 38 meV, and typical de Broglie wavelength 
A = 1.20 A. These values are quite desirable for studying condensed matter. Furthermore, 
such thermal neutrons move at the speed of v ~ 2200 m/s, speeds achievable by mechan- 
ical devices, thus enabling construction of instruments, based on neutron time-of-flight 
measurements, to determine energy and wavelength. 
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3.3 Neutron scattering theory 

The quantum nature of the neutron, its wave and particle character, is well illustrated in 
a scattering experiment. The neutron propagates as a wave, as evidenced by the observed 
diffraction patterns, whereas it interacts with the nucleus practically as a point particle. 
The nuclear interaction can be well described in terms of a contact potential, known as the 
Fermi pseudopotential: 

V( r ) = — 6<5(r). (3.2) 
m 

The quantity b in this equation is known as the scattering length, and is in general a 
complex constant. The real part of b is related to the scattering cross-section: 

dcr ,,, 

= (Re "> 

^totai = 47r(Re6) 2 , (3.3) 

and the imaginary to the absorption of neutron by a nucleus. For most nuclei, the scat- 
tering length is a real constant. Elements such as cadmium or gadolinium have scattering 
lengths that are predominantly imaginary, and correspondingly materials containing these 
elements are virtually opaque to neutrons. Such materials can be used for neutron shielding 
purposes. 

While the interaction of a neutron and a single nucleus results in a simple, spherically 
symmetric, scattered wave, a regular array of nuclei, as are found in crystals, will produce 
revealing diffraction patterns. At this time it should be remembered, as was mentioned 
earlier, that neutrons differentiate not only different atoms, but also atomic isotopes. Since 
the structure of a crystal is determined mainly by the electron interaction, in order to mini- 
mize the free energy atoms arrange themselves in a periodic lattice. However, the nuclei of 
these atoms need not be, and in general are not, of the same type; indeed, different isotopes 
of the same atom are randomly distributed throughout the crystal. This, now, implies that 
neutrons in a crystal will see its periodic structure, giving rise to sharp diffraction peaks, 
as well as the random distribution of isotopes, causing a more uniform scattering. Math- 
ematically, these two contributions can be separated, and the former is named coherent 
scattering, while the later is called incoherent scattering. (More precisely, incoherent scat- 
tering arises from the random distribution of the deviations of the scattering lengths from 
their mean value. Even if the sample consists of a single isotope, nuclei in different nuclear 
states can have different scattering lengths, and so cause the incoherent scattering.) The 
total coherent and incoherent scattering cross sections for an atom are: 

^coh = 47ib 2 C0h = An(b) 2 (3.4) 
a incoh = 47r6f ncoh = 47r((6 2 )-(6) 2 ). (3.5) 

Neutron scattering is described as elastic if the incident energy equals the scattered 
neutron energy, and as inelastic if it does not. The incoherent scattering just described 
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Atom 


Z 


M 




^incoh 


Ccoh 


"mcoh 


°scatt 


Cabs 


H 


1 


1.0079 


-3.7390 


— 


1.7568 


80.26 


82.02 


0.3326 


1 H 


1 


1.0078 


-3.7406 


25.274 


1.7583 


80.27 


82.03 


0.3326 




1 


2.0141 


6.671 


4.04 


5.592 


2.05 


7.64 


0.000519 


V ) 


Q 
o 


1 S QQQ4 


^ RO^ 




4 IX) 


0008 
u.uuvo 


4 IX) 


0001 Q 


Al 


13 


26.9815 


3.449 


0.256 


1.495 


0.0082 


1.503 


0.231 


Sc 


21 


44.9559 


12.29 


-6.0 


19. 


4.5 


23.5 


27.5 


V 


23 


50.9415 


-0.3824 




0.0184 


5.08 


5.1 


5.08 


Cr 


24 


51.9961 


3.635 




1.66 


1.83 


3.49 


3.05 


Fe 


26 


55.845 


9.45 




11.22 


0.4 


11.62 


2.56 


Zn 


30 


65.39 


5.680 




4.054 


0.077 


4.131 


1.11 


Ga 


31 


69.723 


7.288 




6.675 


0.16 


6.83 


2.75 


Rb 


37 


85.4678 


7.09 




6.32 


0.5 


6.8 


0.38 


Zr 


40 


91.224 


7.16 




6.44 


0.02 


6.46 


0.185 


Mo 


42 


95.94 


6.715 




5.67 


0.04 


5.71 


2.48 


Cd 


48 


112.411 


4.87-0.70i 




3.04 


3.46 


6.5 


2520. 


Sn 


50 


118.710 


6.225 




4.871 


0.022 


4.892 


0.626 


Ba 


56 


137.327 


5.07 




3.23 


0.15 


3.38 


1.1 


Gd 


64 


157.25 


6.5-13.82i 




29.3 


151. 


180. 


49700. 


Pb 


82 


207.2 


9.405 




11.115 


0.003 


11.118 


0.171 



Table 3.4: Scattering cross sections for some atoms of interest. Columns: (1) atom; (2) 
charge number; (3) mass (in u ); (4) coherent scattering length (in fm); (5) incoherent scat- 
tering length (in fm); (6) coherent scattering cross section (in barn); incoherent scattering 
cross section (in barn); (7) incoherent scattering cross section (in barn); incoherent scatter- 
ing cross section (in barn); (8) total scattering cross section (in barn); (9) absorption cross 
section for v = 2200 m/s neutrons (in barn). (1 fm=10~ 15 m. 1 barn=10~ 28 m 2 .) Adapted 
from [19]. 



19 



arises from the static crystal. Crystal dynamics results in inelastic coherent and incoherent 
scattering. The former stems from the time correlation between the positions of different 
nuclei, and the incoherent scattering from the time correlation between the position of 
the same nucleus. Whether elastic or inelastic, we can generally say that the coherent 
scattering produces interference effects, and arises from regularities in the crystal (in time 
and space), while the incoherent scattering comes from randomness, thus producing only 
self-interference effects. 

The categories being so defined, we can classify all neutron scattering into the following 
four groups, indicating some of the applications: 

(i) Elastic coherent. Used for measurement of static properties, such as nuclear or mag- 
netic crystal structures. 

(ii) Inelastic coherent. Extremely valuable, as it is often the only experimental technique 
by which crystal dynamics, phonons or magnons, can be studied in detail. 

(iii) Elastic incoherent. Since it is diffuse, and depends only weakly on the scattering 
geometry or crystal structure, and is rather easy to measure. Incoherent elastic scat- 
tering is often used for normalization, instrument alignments, etc. 

O'v) Inelastic incoherent. Enables studying of dynamic phenomena in a more compre- 
hensive (albeit less detailed) manner; for instance, a density-of- states measurements 
(incoherent approximation). 

Elastic coherent scattering, from the nuclei in a crystal, is governed by the following equa- 
tions: 

'%) =A^£*(k-)lW (3.6) 

/ coh el T 

where 

F N (k) = b ~d exp(ik • d d ) exp(-W d ) (3.7) 

d 

is the nuclear structure factor, and 

s 

is the Debye- Waller constant (the quantity e~ 2W is called the Debye- Waller factor), which 
quantifies the thermally induced motion of atoms around their equilibrium position. This 
motion causes an exponential decrease of neutron scattering intensity at high momentum 
transfers, and usually does not play a significant role in the diffraction. In the above equa- 
tions, (da/dQ,) is the differential cross section, the number of neutrons scattered in a given 
direction per second per unit solid angle; N is the number of unit cells in the sample; v is 
the volume of the unit cell; r a vector in the reciprocal lattice; k momentum transferred; b d 
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coherent scattering length of atom indexed with d, whose fractional coordinate within the 
unit cell is d, and mass m d . The polarization tensor e, of a phonon mode indexed with s, 
with energy co s , is derived from a phonon eigenvector, e T = (e x e y e z ), as: 



e e 1 



e e* 

E X 


e e* 


e e* 


e e* 


e e* 


e e* 


e e* 


e z e* 


e z e* z 



and is normalized so that: 

^2^2( e ds)* a (e d s>)p = $s, s > <W; a,/3 = x,y,z. (3.9) 

d a 

Finally, (n s ) is the number of phonon quanta at energy Hlu s in thermal equilibrium at tem- 
perature T, according to the Bose-Einstein distribution ((2n s + 1) = coth(hu s /2kBT)). 

The cross section for the inelastic coherent scattering of a phonon indexed by s is given 

by: 



k f (2tt) 3 N 



h 2v uj s 

coh inel 



e 2ik-d e ik T -?k e -2W (i (k) 



+ (3.10) 




where (n s ) = l/[exp(huj s /k B T) — 1]. 

As mentioned earlier, elastic incoherent scattering has a rather simple dependence on 
the scattering geometry: 

= £E(^coh),e- 2 ^W. (3.11) 

incoh el d 

Since in most scattering experiments the term with the Deybe-Waller factor is close to 
unity, and, being elastic, neutrons so scattered are easily found, this scattering is often used 
for normalizing the measured intensity into physically meaningful units. 

The cross section for neutrons scattering incoherently and inelastically from a phonon 

m) |e— + ,3,2, 

>. / incoh inel d 

Comparison with (3.10) reveals that the main difference from the coherent inelastic scat- 
tering is the lack of the diffraction term, e lk d . This has a profound effect on the scattering 
intensity, which now varies slowly and continuously with k. In a powder sample, the inten- 
sity can be written as a product of purely radial and angular functions, and all the phonons 
on the energy shell contribute to the intensity, regardless of their wavevector. This enables 
experimental measurement of the phonon density of states in a powder sample: 

—) = k lzNY (6i " coh) ^ e-tW |e kT ^- k |l Z(E). (3.13) 

dQ dE i h ) ^ 2m d 1 1 [ Ejh v ; v 

ndos v d ) 
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As the phonon eigenvectors, e ds , are not known a priori, an approximation (k T • e ds ■ 
k) average ~ k/3 considerably simplifies this result: 

J^-) «^LvV ^^ e -^ e )(^±A Z{E) . (3 .14) 
dtt dE i ki I V 2m <i I E * A 

>. / ndos v « J 

Nearly all the terms within the curly braces can be evaluated using only the structure in- 
formation. (The only exception is the Debye-Waller factor, which can be determined 
from a neutron diffraction experiment, as well as with the X-rays. If such measurements 
are not available, it can be estimated using values for similar materials, without introducing 
significant errors into the equation.) 

Neutron can also scatter from two or more phonons before being detected. Although 
rather small, the multiphonon scattering cross section is observable. There are various 
techniques for identifying multiphonon scattering from single crystals, as well as for sub- 
tracting it from a powder spectrum. Since it is a higher order scattering, multiphonon 
scattering does not reveal as much information as a single phonon scattering, so it is usu- 
ally subtracted as an unwanted background. 

Having a magnetic dipole moment, 

eh 

Mn = ~7A<nO-, A<n = " , (3-15) 

neutron also scatters from electrons. Here, 7 = 1.913 is a dimensionless constant; n p is the 
nuclear Bohr magneton; <y is the neutron spin operator (eigenvalue 1, by definition); e > 
is the elementary charge; and m p is the proton mass. The electron, having an equivalent set 
of equations describing its dipole moment: 

eh 

H e = -2/i B s, fi B = - — , (3.16) 
2m e 

thus interacts with the neutron via magnetic dipole-dipole interaction. (Here, s is electron 
spin operator, the eigenvalues of which is by definition 1/2.) Furthermore, being charged, 
an electron with momentum p creates a magnetic field via which it also interacts with the 
neutron. The total potential of interaction can then classically be written as: 

U n —e — U S pj n + t^orbit 

47r |n| 3 

+T~£ 2 ' n = r n -r e . (3.17) 

Even without the full mathematical analysis, one can infer that the resulting scattering 
cross section should augment the nuclear scattering cross section by the following two 
elements: 
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(z) Unlike the nucleus, which can be regarded as a point particle, the electron wave 
function in solids extends over some finite distance. This will result in an additional 
factor in the cross section, called the magnetic form factor, F M (k). 

(ii) Dipole moment has an established direction in space, and the cross section will 
greatly depend on the dipole-momentum transfer geometry. This will be made ex- 
plicit by the first summation, J2 a p> m equation (3.19), below. 

In many magnetic ions of interest, the orbital angular momentum is zero, being ei- 
ther quenched by the crystal electric field (nearly all 3d ions), or annulled by a half-filled 
shell (Mn 2+ , Fe 3+ , La 3+ , Gd 3+ , Lu 3+ ). A general theoretical analysis, for ions with non- 
vanishing orbital angular momentum (as in most rare earth ions), is quite complicated [20]. 
However, when the mean radius of the (unpaired) electron wavefunction is smaller than 
|k| _1 , a so called dipole approximation can be utilized, and the two cases are described 
by essentially the same set of equations. As this generalization involves the electron wave 
function, one can assume, following the reasoning implied in statement (z) above, that the 
cross section will contain a magnetic form factor; as is indeed the case. 

With this in mind, we now write the elastic magnetic scattering cross section, assuming 
— for simplicity — that the magnetic lattice consists of magnetic ions of the same element, 
as is often the case in practice: 




^ M (k) 



e -2W(k) j 



(3.18) 
(3.19) 



I = Y1 OW-kA) ^2(Sp)(Sf)^- v \ a,0 = x,y,z. 

a,P 1,1' 

In addition to the terms previously defined, we have r = n e 2 /4irm e = 5.38 x 10~ 15 m, 
the classical radius of electron; and (Sfi) is the a component of the average moment of 
atom /'. The magnetic form factor, _F M (k), can be written as: 



^M(k) = - 
9 



9sJo(k) + g L (jo(k) + J 2 {k)) 



where g is the Lande gyromagnetic ratio, given as: 

5(5 + l)-L(L + l) 
J(J + 1) 

1 L(L + 1) -5(5 + 1) 

2 + 



9s 
9l 
9 



2J(J + 1) 
_3 5(5 + 1) - L(L + 1) 
+ 2 J(J + 1) 



(3.20) 



(3.21) 
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Ion 


S 


n 


,4 


a 


5 


b 


C 


c 


D 


Cr 3 + 


3/2 





-0.3094 


0.274 


0.3680 


17.0355 


0.6559 


6.5236 


0.2856 






2 


1.6262 


15.0656 


2.0618 


6.2842 


0.5281 


2.3680 


0.0023 


Fe 3+ 


5/2 





-0.3972 


13.2442 


0.6295 


4.9034 


-0.0314 


0.3496 


0.0044 






2 


1.3602 


11.9976 


1.5188 


5.0025 


0.4705 


1.9914 


0.0038 



Table 3.5: Magnetic form factors. For use in Eq. (3.23). A, B, C, and D are dimensionless; 
a, b, and c are in A 2 ; and k is in A -1 . [22] 

The functions J n {k) are defined as: 

poo 

J n (k)=4ir / drr 2 j n (kr) |^(r)| 2 , (3.22) 
J o 

where j n (kr) is a spherical Bessel function of order n, while |V>( r ) P * s me (unpaired) elec- 
tron's radial probability density, as defined by Born. These functions, J n (k), are usually 
given in an approximate form as a sum of Gaussians: 

J n=0 (k) = Ae- ak2 + Be- bk2 + Ce~ ck2 + D 

J n=2 (k) = (Ae~ ak ' + Be- bk2 + C^ ck2 + D)k 2 . (3.23) 

The constants for many ions can be found in the literature [21], and Table 3.5 lists values 
for some of the ions relevant to this thesis. 

It is, perhaps, not obvious from the equation (3.19), but very important for an experi- 
menter, that the magnetic scattering intensity vanishes if momentum transfer, k, is aligned 
with the dipole moment, (Si) = ((Sf Sf 5f) T ). Thus, we can extend the statement (ii) (on 
p. 22) by adding that the neutron scatters only from magnetic dipoles (or the components 
thereof) perpendicular to its momentum transfer (there generally being two such directions, 
or polarizations), while no scattering occurs for longitudinal polarization. 

As a corollary of the equations above, we can now note the cross section for ions with 
spin only dipole moment. With L = 0, the total angular momentum equals that of the spin 
only: 

J 2 = S 2 + L 2 + 2S-L 

= S 2 , forL = 0, (3.24) 

and the magnetic form factor reduces from equation (3.20) to: 

F M (k) = j d 3 r |^(r)| 2 e ik ' r = J (k). (3.25) 

As the theory given so far is adequate for the purpose of this thesis, we will shortly 
close this section; but not before mentioning that the story of magnetic scattering does not 
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end here. The general equations for the inelastic magnetic scattering are rather lengthy 
for this presentation, and interested reader should consult some of the excellent existing 
literature [24, 25]. Suffice it to say that the contributions include the terms with both the 
phonon and magnon excitations, magneto-vibrational scattering being one such example. 
In practice, equations describing magnetic scattering, such as (3.18), can often be reduced 
to quite simple forms (as will be seen later in this thesis), once a particular magnetic model 
is presumed, and summations carried out. 

Furthermore, we have assumed in the discussion above that electrons are localized 
around their host atoms. While this is quite an accurate picture for ionic crystals and 
rare earths, in transition metals conduction electrons play a dominant magnetic role, and 
their wave functions span the entire sample. The corresponding magnetic scattering equa- 
tions differ considerably from those already presented. Liquids, which are also studied by 
neutrons, present an even more extreme case. We also have not yet mentioned polarized 
neutron scattering, in which the neutron spin state is measured in addition to its momentum 
and energy; but by now the reader should have gained a sense of the breadth of the theory 
needed to describe neutron scattering from condensed matter. 



3.4 Experiment 

It is evident, from the equations listed in the previous section, that in order to use 
neutrons to study condensed matter, an experimenter must know the amount of energy and 
momentum it exchanges with the sample; E and k are present in nearly all the equations. 
Luckily, these quantities are rather easy to determine: One can employ Bragg diffraction, 
and the usual conservation laws — all well known from the realm of general physics — 
and convert from the wave to particle picture and vice versa, as need may be. The later 
is easily done with the usual relations linking the energy, frequency, momentum, velocity, 
and wavelength of a massive particle: 

E = hoo = £- = —— = —[ — ) = - — . (3.26) 
2m 2 2mV27rA/ 2m 

Clearly, knowing any one of these quantities reveals all the others; with an exception that 
only the magnitude of vector quantities, p, k, and v, can be determined from the scalar 
ones. The direction is established experimentally through the use of collimators or position 
sensitive detection. 

In a triple axis spectrometer (TAS) — the most versatile neutron scattering instrument 
— among all the neutrons coming from the guide, those with the desired energy and mo- 
mentum pass through the collimator, are reflected off the monochromator, and go through 
the second collimator on the way to the sample, all others being absorbed or reflected else- 
where. These incident neutrons are then scattered off the sample, and in the same manner, 
if they have the desired final energy and momentum, will pass through the collimator- 
analyzer-collimator path before being counted by the detector. Strictly speaking, we do not 
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Table 3.6: Monochromator d- spacing. [26] 



measure momentum nor energy of a neutron in a scattering experiment; instead, using a 
spectrometer we select only those neutrons that have the desired momentum and energy. 
In a triple axis spectrometer — the components of which we will presently describe — the 
monochromator is used to select the neutrons of some desired energy, the incident energy, 
and have only them interact with the sample. Emanating from a neutron guide, with a speed 
distribution given by Eq. (3.1), neutrons are elastically reflected from the monochromator 
according to Bragg 's law: 

A = 2d hH sin0, (3.27) 

or equivalently: 

29 

r hkl = 2k any, (3.28) 

relating the neutron wavelength and wavevector, X and k, to the scattering angle, 29. (Here, 
d hk i is the spacing between the {h, k, 1} lattice planes, and k hk x = 2ir/d hk i.) 

A good monochromator must have a large coherent scattering cross section, and small 
incoherent and low absorption cross section. It must also be possible to produce large 
single crystals, with a mosaic distribution of order 0.5°. Typical materials are copper, 
silicon, germanium, and a layered material called Pyrolytic Graphite (PG). PG consists of 
hexagonal carbon layers, with highly parallel (00Z) axes, but randomly oriented (hOO) and 
(OkO) directions. Thus there are sharp Bragg peaks only along the (00Z) axis, and powder- 
like reflections elsewhere. With an incoherent cross section of only 2% of the coherent 
cross section, and the peak reflectivity higher than that of single crystals, PG is a nearly 
ideal material for a monochromator. 

By using a collimator, a device consisting of parallel absorbing planes, only neutrons 
with the desired scattering angle, i.e. energy, are directed toward the sample for scatter- 
ing. Ideally, all others are absorbed by the collimator, which commonly consists of steel 
blades coated with cadmium (which has large neutron absorption cross section), inserted 
into machined slots in the neutron guide. By adding or removing blades, the experimenter 
can control the divergence of the beam, which is usually set between 20' and 2° (when 
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no blades are inserted), smaller values giving better defined orientation and energy of the 
beam, and larger increased neutron intensity. There is a collimator before and after the 
monochromator, together defining the precision of the monochromator scattering angle, 
26m- 

The axis on which the monochromator is located is the first of three in a triple axis 
spectrometer (hence the name), the sample being located at the second one. Neutrons scat- 
ter from the sample generally in all directions; orientation of the third collimator, after the 
sample, determines which of these neutrons will be detected. The configuration in which 
the detector is located after this collimator is called two axis spectrometer, and predates the 
triple axis version by a decade. Set like this, the spectrometer detects neutrons with all the 
final energies. Since the elastic scattering is much stronger than inelastic one, it is usually 
said that the two axis configuration detects only the elastic scattering; although not true, 
the distinction is in practice negligible. 

If knowing the final energy is desirable, as is often the case, the neutrons are reflected 
once more, off a device called an analyzer, located on the instrument's third axis. The ana- 
lyzer is essentially identical to the monochromator: the later selects the neutrons according 
to their incident, while the former according to final energy, both using Bragg diffraction. 
Yet another collimator, the fourth one, defines the orientation and the divergence of the 
beam that will finally enter the detector. 

Thus, whenever there is a neutron detected in the spectrometer, we know not only its 
incident and final momentum and energy: 

TIT 

kij = - . , (3.29) 
2sm 6ij 

h 2 k 2 

E if = A (3.30) 

J 2m V 

but also the uncertainties in these quantities, through the mathematical concept of the res- 
olution function. Use of the conservation of energy: 

E = Ef — Ei, (3.31) 

and momentum: 

k = k f - ki (3.32) 

then reveals the energy and momentum the neutron exchanged with the sample, the quan- 
tities we need to interpret the scattering results. Wavevectors in the Eq. (3.32) form the 
scattering triangle, which is in practice solved using the cosine rule for triangles: 

26 

k 2 = k 2 + k 2 f -2kik f cos— , (3.33) 

where 26 s is the angle at which the neutron scatters off the sample. 

Neutrons are collected and counted in a device called a detector, which usually con- 
sists of an aluminum chamber filled with a gas, such as 3 He or 10 BF 3 . Since neutrons 
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only weakly interact with matter, and carry no charge, the detector uses nuclear reactions 
initiated by the neutron, such as: 

n + 3 He — > p + 3 H + 0.77MeV 

n + 10 B — > a + 7 Li + 7 + 2.78 MeV 

to produce charged particles. The energy so released in single reactions is sufficient to 
produce over 10 4 electron-ion pairs, and since it takes place in a strong electric field, ~ 
100 V/cm, a cascading ionization of the gas amplifies this current by a factor of 10 5 , to 
levels at which it can be detected electronically. Since this current must abate before the 
next neutron can be detected, there is an upper limit to the counting rate, and is of the 
order of 10 4 to 10 5 counts per second. [It is, perhaps, interesting to note, that although 
the this nuclear reaction is irreversible, the detector continuously running at this, saturated, 
counting rate, for an unrealistically long time, say a hundred years, would cause only a 
small fraction of contained nuclei, about 10 7 , to participate in the reaction. Therefore, 
the efficiency of the detector remains practically unchanged during its lifetime, and the gas 
inside the chamber needs not be exchanged.] 

It is imperative for an experimenter to be aware that the absorption cross section, and 
therefore the efficiency of the detector, is inversely proportional to the velocity of thermal 
neutrons: 

Detector efficiency ~ cr abs0 rption ~ - ~ -^=. (3.34) 

To counter this detector bias toward the lower energies, spectrometer is usually operated in 
the fixed final energy configuration, where variable energy transfer is achieved by varying 
the incident energy. The situation is, however, complicated by the fact that the reflectivity 
of the monochromator is dependent on the neutron energy, and although the reflectivity 
for different monochromators is known, its dependence on the energy is not trivial. The 
remedy is a device called monitor, essentially a very low efficiency detector (achieved by 
decreasing the pressure of the gas in the chamber), through which most of the neutrons will 
pass unaffected, while those detected reveal the intensity of the neutron beam. Naturally, 
the monitor is placed right before the sample, as it is the intensity of the beam impacting the 
sample that we are interested in. As a result, monitor also serves to correct any unsteadiness 
in the reactor power, although this was more important with the early reactors than is today. 

Let us close this description of spectrometer parts with another important device often 
placed in the beam path: filter. Inspection of Eq. (3.28) shows that if neutrons of wavevec- 
tor k are scattered at an angle 29 M , so will be all which wavevector is an integer multiple 
of k, because T( nhtnktn i) = nr^^i). The intensity of these, higher order, reflections is 
proportional to their nuclear structure factor, and for some lattices this intensity can iden- 
tically vanish. For instance, silicon, which crystallizes in a diamond lattice, has significant 
(111), but no (222) reflection, and therefore no contamination from A/ 2 neutrons, as they 
are usually referred to. Although A/3 (and higher order) neutrons will still be present, their 
energy will be 9 = 3 2 times the energy of the desired, A, neutrons, energies at which the 
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Table 3.7: Bragg scattering filters. [26] 

neutron flux already abates. Unfortunately, not all good monochromators have vanishing 
higher order reflections; for instance, PG reflects strongly both for (002) and (004), two 
most commonly used reflections. Contamination in these cases is eliminated by placing 
a filter in the beam. Most common are Bragg-scattering filters, in which a powder of a 
material with well known lattice spacing, large coherent, and small incoherent and absorp- 
tion cross sections, is placed inside an aluminum container, and cooled by liquid nitrogen 
to minimize inelastic phonon scattering. These low-energy pass filters will reflect, and so 
eliminate from the beam, most neutrons with the wavevector larger than 2ir/d max = r m j„, 
because for most k > T m i n there will be a strongly scattering Bragg peak somewhere in 
the reciprocal space (d max is the largest spacing between lattice planes in the crystal). Two 
such commonly used filters are Beryllium, which transmits neutrons below 5.2 meV, and 
BeO, with cutoff energy 3.7 meV. 

PG, which we already know as a monochromator, is also widely used as a filter; placed 
perpendicularly to the beam, it will reflect back neutrons with wavelength A = 2dpo (this 
follows from eq. (3.27) for 9 = 90°), and so will have deep dips in transmission profile. 
Well known dips occur at 13.7 (with the transmission of only 1CT 4 !), 14.7, and 30.5 meV, 
and often have decisive influence on the choice of neutron energy used in the experiment. 

The components just described are all commonly found in a triple axis spectrometer, 
arguably the most versatile neutron instrument, as it gives the experimenter a lot of control 
over the access to the reciprocal space, size and shape of the resolution ellipsoid, incident 
and final energies, etc. (The resolution ellipsoid is a 4-dimensional ellipsoid in the recip- 
rocal k-u; space, inside which the resolution function exceeds the half of the maximum 
value.) Its prime disadvantage is that it probes the reciprocal space one point at a time, 
so it cannot scan large volumes efficiently. If the classic neutron detector is replaced by a 
position sensitive detector (PSD), which records not only neutron count, but also the loca- 
tion at which they were detected, larger volumes can be probed more efficiently. If even 
larger volumes need be probed, multi-detector time-of-flight (TOF) spectrometers are em- 
ployed, which discriminate neutrons by the scattering angle and the time of detection, thus 
revealing their wavevector and energy transfer. 

To produce the monochromatic beam, TOF spectrometers use choppers, which come 
in two designs: (?) Fermi chopper spectrometer (FCS): very simple in design, essentially 
consisting of a rapidly rotating cylinder with the grove machined in it, curved so that only 
neutrons of selected velocity can pass without hitting its walls, coated with absorbing ma- 
terial. By changing the angular frequency of the cylinder, incident energy can be changed, 
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Figure 3.1: (a) Triple axis spectrometer. Green line shows incident, and blue line scattered 
neutron path, (b) Scattering triangle. 
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Figure 3.2: Thermal neutron instruments at NCNR. (a) Triple axis spectrometer BT-2. 
(b) Filter Analyzer Neutron Spectrometer (FANS). Source: NCNR. Reproduced with 
permission. 




Figure 3.3: Cold neutron instruments at NCNR. (a) High Flux Backscattering Spectrome- 
ter (HFBS). (b) Spin-Polarized triple-axis Spectrometer (SPINS). Source: NCNR. Repro- 
duced with permission. 
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Figure 3.4: Guide hall with cold neutrons instruments at NCNR. Among the instruments 
shown, the following were used to collect the data for this dissertation: (a) Spin Polarized 
Triple Axis Spectrometer (SPINS); (b) Fermi Chopper Spectrometer (FCS); (c) Dies Chop- 
per Spectrometer (DCS); (d) High Flux Backscattering Spectrometer. FCS and DCS are 
time- of- flight (TOF) instruments. Source: NCNR. Reproduced with permission. 
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although in practice a monochromator is placed before the chopper for better energy defi- 
nition, (ii) In a disc chopper spectrometer, a series of rapidly rotating aluminum discs, their 
surface coated with neutron absorber except on narrow "windows," are used to monochro- 
mate the neutrons coming from the beam. (These windows are not actual openings on the 
disks, but only areas without the absorbing coating. Rotating at 20,000 rpm, any asymme- 
try in the disc would create severe vibrations.) As a consequence, choppers produce pulsed 
beam, as opposed to continuous in a triple axis. Knowing the initial energy of the neu- 
tron, and the time it passed the last chopper, energy and momentum transfer can simply be 
calculated from the instrument geometry. By using numerous detectors, large volumes of 
reciprocal space are probed; however, having a fixed geometry (incident energy is the only 
variable for the TOF spectrometers), it cannot be tailored to suit a particular experiment 
or purpose. Therefore, we conclude that the TAS and TOF instruments complement each 
other, both having distinct advantages and applications. The spallation sources, producing 
the intrinsically pulsed beam, favor the TOF instruments. However, by combining the two 
instruments (for instance, by using the time structure of the beam to eliminate the analyzer) 
one may get the best of both worlds. 

In the remaining part of this section, we will only briefly mention two (among many) 
instrument designs, the two that were used to collect some of the data presented in this 
thesis. First of these, the backscattering technique, is based of the notion that the energy 
resolution of the neutrons scattered nearly backwards (i.e. 29 »s 180°), is very small. By 
combining the Eqs. (3.29) and (3.30), we can write: 



h 2 (nT)' 
2msin 2 (2#/2) 



E = ~ • 2 ,L ^ - ( 3 - 35 ) 



Assuming the uncertainties in 9 and r to be small, the energy resolution equals: 



AE = dsfe v cot2(M/2) <M)2 + ? (AT)2 - <3J6) 

We see that in the backscattering configuration, as 29 — > 180°, the energy resolution, be- 
comes very small. Indeed, in the very limit, it is only limited by the lattice imperfections, 
Ar/r. In practice, energy resolutions of less than 1/ieV are achieved, enabling study of 
nanosecond range dynamics. 

We already mentioned one drawback of TAS instruments, namely, they probe the sam- 
ple "one point at a time", and neutron scattering in general, in its relatively slow acquisition 
of data. Focusing monochromators and analyzers are designed to address some of these is- 
sues. Focusing analyzer is a device in which a single analyzer, a PG plate for instance, is 
replaced by a number of smaller ones, each separately oriented to reflect the neutrons of 
the same energy, effectively increasing the size of the analyzer, and therefore the neutron 
flux on the detector. (Essentially the same in design, focusing monochromator increases 
the effective area of the monochromator.) Now, an increased area implies a higher flux, 
but also a greater uncertainty in the neutron scattering angle; and indeed, the uncertainty 
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in the momentum transfer becomes much larger for these instruments. However, by being 
able to independently orient the elements of the analyzer, an experimenter can tune the 
resolution function and split the uncertainty between the momentum and energy, extending 
it into some irrelevant part of the reciprocal space, and in turn more efficiently measuring 
some relevant quantities. Furthermore, by using a position sensitive detector, an additional 
variable (detection position, which translates into scattering angle) can be utilized to en- 
hance mapping of the reciprocal space to "one line at a time." Since here I cannot go 
into great detail describing focusing instruments, I'll conclude this section by expressing 
my fortune to have an advisor involved in the design and construction of two such instru- 
ments: the focusing analyzer for SPINS, which I was fortunate to have used; and the still 
ongoing doubly focusing monochromator project (cold neutron (NG-0) Multi Axis Crystal 
Spectrometer (MACS), at NCNR), in which I shortly participated. 

3.5 Techniques 

While the underlying scattering principles are the same, there are different types of 
instruments optimized for different types of samples and physical problems. Criteria used 
to select the right type of instrument are listed below 

(z) For powder samples, TOF and multidetector techniques are generally preferred over 
TAS, while the opposite is true for single crystals. Powder samples are generally 
isotropic, and then some of the parameters that can be controlled by a TAS become 
irrelevant. Features in a powder spectrum are less sharp than those in a single crystal, 
and usually extend over a greater range in Q, making the use of TOF more preferable. 

(if) If details in the excitation spectrum are sought, one should use a single crystal on 
a TAS. In contrast, powder on a TOF instrument will give comprehensive spectrum 
overview, and is a good initial experiment on a new sample. 

(Hi) If the excitation spectrum lies in higher energy range (say, above 20 meV), ther- 
mal neutrons are preferred to cold, and vice versa. For very low energy transfer 
(> 1 /ueV), backscattering technique should be used. 

[Increasingly, a new generation of spallation source instruments with very large detector 
banks are becoming available. They enable access to extensive regions in reciprocal space 
that can be mapped quickly and precisely. The MAPS spectrometer at ISIS, for instance, 
optimized for thermal neutrons, has a 2.6 steradian detector banks with nearly 150,000 
pixels. In addition, a wide range of incident energies from 15 to 2000 meV, make this a 
versatile instrument with many applications. However, no data in this thesis were obtained 
using such instruments, so I will accordingly exclude them from further discussion.] 

Experiments on powder samples are considerably simpler than those on single crys- 
tals. Due to the isotropy of their reciprocal space, energy and magnitude of the wavevector 
transfer are the only variables. As the features in the excitation spectrum are less sharp, the 
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effects of the resolution function are also less prominent, and simpler to account for. In con- 
trast, scattering space for a single crystal is four dimensional, although usually only three 
are accessible in any particular setting on a TAS. Resolution, extinction, (mis)alignment, 
and others effects all contribute significantly to the measured intensities, and have to be ac- 
counted for. For such experiments in particular, a compromise must be struck between high 
beam intensity, which would decrease counting time, and high resolution, which makes the 
measured features sharper and more informative. 

Following are the instruments that were used to collect the data for this dissertation. 
Among the thermal neutron instruments, we used (i) BT-2, and (ii) BT-9 triple-axis spec- 
trometers (TAS); (in) BT-4 Filter Analyzer Neutron Spectrometer (FANS); and (iv) BT- 
1 High-Resolution Powder Diffractometer (HRPD). Cold neutron instruments used were 
time-of-flight (TOF) spectrometers (v) Fermi Chopper Spectrometer (FCS), and (vi) Disc 
Chopper Spectrometer (DCS); (vii) Spin-Polarized triple-axis Spectrometer (SPINS); and 
(viii) High Flux Backscattering Spectrometer (HFBS). 



36 



Part II 
Frustrated Systems 



37 



Chapter 4 

Negative Thermal Expansion in ZrX^Og 
and Sc 2 (W0 4 ) 3 

4.1 Introduction 

In 1996, a large (a — 10~ 5 K _1 ), isotropic, Negative Thermal Expansion (NTE), over 
a wide temperature range (from 0.3 to 1050 K, linear from 50 to 420 K), was observed in 
zirconium tungstate (ZrW 2 8 ) [27]. The discovery attracted a considerable interest, both 
for its potential applicability, and because of the need for a theoretical explanation. As a 
result, other materials with large (though anisotropic) NTE were discovered soon thereafter 
[28,29,30,31]. 

Potential applications for NTE materials are numerous, ranging from industrial (diffrac- 
tion gratings, circuit boards) and medical (dentistry), to household applications (kitchen- 
ware). Composites in which NTE materials are combined with PTE (positive thermal ex- 
pansion) materials, offer a precise control over the magnitude of thermal expansion (zero 
thermal expansion often being the most desirable). 

Theoretically, the phenomenon is fascinating. Until these recent discoveries, no other 
material was known to exhibit so large a NTE over a wide temperature range. NTE is 
known to exist in many materials, but is usually confined to a very limited temperature 
range, often around a phase transition (a well known example being that of water near the 
freezing temperature). As such, NTE can be considered an anomaly, in contrast to the 
positive thermal expansion, commonly found in most materials. Figure 4.1 illustrates the 
magnitude of the thermal expansion in ZrW 2 8 and Sc 2 (W0 4 ) 3 , compared to that in some 
other materials. 

All of the recently discovered NTE materials belong to two iso- structural families, 
AB 2 8 , and A 2 (B0 4 ) 3 , where A=Zr, Sc, Hf, Y, or Lu, and B=W, or Mo (although not 
all the combinations are possible). They all share a structure in which the oxygen atoms 
are in a twofold (B-O-B, or more often A-O-B) coordination to the metal cations. Based 
on this observation, a mechanism was proposed in which the transverse thermal motion of 
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oxygen atoms causes the contraction of the unit cell [32]. (This is the origin of a negative 
thermal expansion in Cu0 2 [33], or Rbl [34].) As long as the metal-oxygen bonds are 
strong enough, so that the distance between them increases negligibly with the tempera- 
ture, the transverse oxygen motion can cause the net negative thermal expansion. Figure 
4.2 illustrates this mechanism in a simple two-dimensional model. Notice that in this model 
the motions of oxygen atoms need not be correlated. Such correlations, however, do exist 
in real NTE materials. 

The structure of all these new oxide NTE materials consists of corner sharing quasi- 
rigid units: A0 6 octahedra, and B0 4 tetrahedra. Strong metal-oxygen bonds, and the 
repulsive oxygen-oxygen forces, maintain rigidity of these polyhedra, which are coupled 
only by their corners and are able to rotate (relatively) freely. Phonon modes in which 
these quasi-rigid units rotate undistorted are called rigid-unit-modes (RUM). A theory was 
proposed in which collective rotations of these rigid units are responsible for NTE [35, 38]. 
Figure 4.3 illustrates a simple two-dimensional model in which the NTE is driven by a 
rigid-unit-mode. It is still the large transverse motion of oxygen atoms that causes the 
NTE; however, in the RUM approach such motions are strongly correlated. 

The existence of quasi-rigid units is necessary, but not sufficient, for rigid-unit-modes. 
The structure also must be under-constrained. A structure is defined as under-constrained 
if the total number of degrees of freedom is larger than the number of constraints. The 
difference between the two equals the number of the rigid-unit-modes. In three dimensions, 
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Figure 4.2: A possible mechanism for the negative thermal expansion, involving large 
transverse thermal motions of atoms with two-fold coordination (usually oxygen). The 
transverse atomic motions need not be correlated to yield the NTE effect. In ZrW 2 8 and 
Sc 2 (W04) 2 , however, the transverse motions are highly correlated, leading to so-called 
rigid-unit-modes. (Compare with Fig. 4.3.) 




Figure 4.3: Simple two-dimensional model for negative thermal expansion caused by rigid- 
unit-modes. Rotation of rigid units causes contraction of the unit cell. 




Figure 4.4: Crystal structure of ZrW 2 8 . Shown are the following quasi rigid units: Zr0 6 
octahedra in green, and WO4 tetrahedra in blue. Notice that the quasi-rigid units are linked 
through corners only. 
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Figure 4.5: Phonon density of states in alumina (A1 2 3 ). Notice the lack of intensity at low 
energies (below « 30 meV). Compare with Figs. 4.6 and 4.11. 



there are six degrees of freedom per rigid unit. Generally, there are three constraints per 
link between rigid units, although the crystal symmetry can decrease the effective number 
of constraints. While there are many structures in which one can identify quasi-rigid units, 
in most cases these are over-constrained. 

As the rigid-unit-modes cost very little energy, their presence is revealed by an unusu- 
ally large phonon Density-Of-States (DOS) at low energies (say, E < 30 meV). Figure 
4.5 shows the DOS of alumina (A1 2 3 ), in which the rigid units consist of aluminum ions 
surrounded by oxygen in octahedral configuration. Alumina has an over-constrained struc- 
ture, and accordingly lacks significant phonon DOS at low energies. In contrast, Fig. 4.6 
shows a large DOS downshift in NTE material zirconium tungstate. A large phonon DOS 
at low energies is a signature of all recently discovered NTE materials. 

Low energy {i.e. soft) phonon modes usually signal, and drive, structural phase tran- 
sitions. Such a phase transition can sometimes cause a large contraction in the volume of 
the unit cell. An example is the ferroelectric phase transition in barium titanate (BaTi0 3 ), 
which is driven by a soft optical phonon mode, where Al/l k -2 x 10~ 4 between T = 
110—135 °C, yielding an average NTE of a « -10~ 5 K 1 (see Fig. 4.7) [40, 42]. The 
presence of soft modes in barium titanate, and the recently discovered materials, suggests 
that in both the NTE may be driven by a qualitatively similar mechanism. Why, then, is 
NTE in BaTi0 3 limited to a narrow temperature range, unlike in ZrW 2 8 ? Analysis in- 
dicates that the symmetry of the low energy phonon modes in ZrW 2 O g is lower than the 
crystal symmetry [35]. This frustrates the phase transition, and consequently, the narrow 
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Figure 4.6: Phonon density of states (DOS) in negative thermal expansion (NTE) material 
ZrW 2 8 . Bullets (•): TOF; open circles (o): FANS. Notice the large DOS at low energies. 
Adapted from [39]. 
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Figure 4.7: Unit cell size, and the thermal expansion in BaTi0 3 . Notice a large negative 
thermal expansion (NTE) around the ferroelectric structural phase transition (T ~ 130 °C). 
If the soft modes leading to the transition were frustrated, the temperature range over which 
the NTE occurs may considerably broaden. Reproduced from [40]. 



NTE temperature region (as seen in BaTi0 3 ) is broadened. Thus, the essential ingredients 
for NTE in ZrW 2 8 seem to be (i) the under-constrained structure, leading to a soft rigid- 
unit phonon modes, and (ii) the frustration of a volume-reducing phase transition driven by 
such modes. 

To support the hypothesis that NTE in ZrW 2 8 is driven by low energy phonon modes, 
we compare the temperature dependence of the unit cell size with the phonon density of 
states (Fig. 4.8). We see that above T »s 50 K, the NTE is linear. This suggests a mechanism 
characterized by an energy scale of about E »s 5 meV. Above 50 K, phonons with higher 
energies are activated and start contributing to the thermal expansion. Phonons at energies 
significantly below 5 meV are apparently less important as there is a reduction of the effect 
for T < 50 K. We see that the phonon DOS rises sharply around 5 meV, and so we can 
identify these as the phonons that are relevant for the NTE [39] . 
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Figure 4.8: Negative thermal expansion (o), and phonon density of states (•) in ZrW 2 8 . 
Energy (top) and temperature (bottom) scales are equivalent (1 meV=11.6 K). Notice the 
reduction of the NTE effect for T < 50 K, and the sharp rise in the phonon density of states 
around 5 meV. This suggests that phonons at energies significantly below 5 meV are less 
important for the NTE effect. Adapted from [39]. 
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4.2 Experiments 

Phonon density-of-states at high pressure in ZrW 2 8 

To more directly investigate the role of the low energy phonons, we measured the 
phonon density of states in ZrW 2 8 under high pressure. If a conventional sample is placed 
under high pressure, the frequency of its phonon modes generally increases. This is param- 
eterized by the so-called Griineisen parameter, 7, of anharmonic lattice theory, which is 
defined as: 

7 = 5^- (4-D 

Here, q is the specific heat, and 7, the Griineisen parameter of a phonon mode with fre- 
quency uji. 

ain^ 

while V is the volume of the unit cell. In conventional systems, the Griineisen parameter is 
positive: 7 > 0. Since obviously dV/dp < 0, subjecting a sample to high pressures tends 
to increase the phonon frequencies. The following equation relates the thermal expansion 
coefficient, a, to the Griineisen parameter, 7: 

We see that the NTE materials must have 7 < (the heat capacity C v , and the bulk modulus 
B, are always positive). Consequently, there must be phonon modes the frequency of which 
decreases with increasing pressure: 

and these modes are relevant for the negative thermal expansion. Figure 4.9 shows the 
phonon density of states at ambient pressure, at p = 1.6 kbar, and p = 2.8 kbar. We see 
that at increased pressures the DOS decreases for E > 4 meV, while it increases below 
3 meV. The changes are qualitatively consistent with the spectral downshift predicted by 
Eq. (4.4), and the conclusion that phonons close to and above 5 meV cause NTE. 



Phonon density-of-states in Sc 2 (W0 4 )3 

The remainder of this section describes our experiments on scandium tungstate. Sc 2 (W0 4 ) 3 
crystalizes in an orthorhombic lattice, space group Pnca, with the unit cell size a = 
9.580 A, b = 9.672 A, c = 13.318 A. The crystal structure of Sc 2 (W0 4 ) 3 is shown in 
Fig. 4.10. For clarity, Fig. (b) shows only the top b-c layer of the structure. Notice the 
openness of the structure, which reduces the number of constraints in the system. The 
quasi-rigid units are linked only by their corners, thus enabling rigid-unit modes. (Notice 
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Figure 4.9: Phonon density of states zirconium tungstate (ZrW 2 8 ) at high pressures. Solid 
line (•) represents DOS at ambient pressure, dotted line (A) at p = 1.6 kbar, and dashed 
line (o) at p = 2.8 kbar. The density of states decreases with increasing pressure at en- 
ergies above 4 meV, while it increases below. This indicates that the low energy phonons 
have negative Griineisen parameter, and therefore are important for the negative thermal 
expansion. 




b 

Figure 4.10: Structure of Sc 2 (W04) 2 . (a) Unit cell, (b) The top b-c layer in the unit cell. 
Shown in green are WO4 tetrahedra, and in blue Sc0 6 octahedra. Compare with Fig. 4.3. 
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Figure 4.11: Density-of-states ("neutron weighted," i.e. NDOS) in powder Sc 2 (W0 4 )2. Cir- 
cles show NDOS measured in a neutron scattering experiment. Solid line is a calculated 
NDOS (instrument resolution applied). Measured using time-of-flight (TOF) Fermi Chop- 
per Spectrometer (FCS), and Filter Analyzer Neutron Spectrometer (FANS) at NCNR. 



the qualitative similarity between the Sc 2 (W0 4 ) 3 structure and the simple two-dimensional 
model in Fig. 4.3.) 

Figure 4.11 shows the phonon density of states in Sc 2 (W0 4 ) 3 , measured by powder 
inelastic neutron scattering. As in ZrW 2 8 , we see a large spectral downshift, consistent 
with the existence of rigid-unit- modes. Notice also the great similarity, even quantitative, 
between the two spectra. The spectrum is divided into a low and high energy regions by a 
gap at around 80 meV. Low energy part consists of rigid- and near-rigid-unit-modes, while 
the phonons in high energy part distort the Sc0 6 and W0 4 polyhedra. 

In powder measurements, information about phonon modes is spherically averaged, and 
therefore not very detailed in Q space. To determine the wave-vector dependence of the 
phonon scattering intensity, we used inelastic single crystal neutron scattering, and scanned 
the Q = (h k 0) reciprocal plane in the E = 6.0 ± 2.2 meV energy window. The result 
is shown in Fig. 4.15. We found a region of increased intensity localized in an irregular 
arc extending approximately from Q k, (6 2 0) to Q »s (280) (centered in the ori- 
gin). The phonon structure factor has a term proportional to Q 2 ; thus the intensity rapidly 
diminishes as Q — >■ 0. At high wave-vector transfer, the Debye-Waller factor causes ex- 
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Figure 4.12: "Neutron photographs" of Sc 2 (W04) 2 single crystal in a neutron beam. The 
sample is a disk-like object seen in the center of each photograph. Such photographs are 
made by placing a neutron camera in a well collimated direct beam immediately after the 
sample. They are used for centering the sample. 



ponential extinction of intensity, proportional to exp( — 2W) ~ exp(— Q 2 (u 2 )/3). Hence, 
the maximum intensity is confined to a "ring" in reciprocal plane, at the center of which 
Q 2 ex-p(—Q 2 (u 2 )/3) is near its maximum. The actual intensity distribution also depends 
on the structure factor, which has an interference term that varies rapidly and unevenly with 
Q, making the "ring's" appearance irregular. 

Phonon dispersion relations in Sc 2 (W04) 2 

To map the phonon dispersion relations, we measured energy spectra for fixed wave 
vector transfer in the regions of the highest intensity along the ring. The one-phonon neu- 
tron scattering cross section depends strongly on the phonon polarization. Indeed, only 
longitudinal atomic motions (parallel to the wave-vector transfer, Q) contribute to the in- 
tensity. To observe modes with different polarizations, we mapped each dispersion relation 
in two different parts of the reciprocal space. Those were chosen so that the wave-vector 
transfer is nearly perpendicular between the two. Figure 4.17 shows a constant-Q scan at 
Q = (h 8 0), h = 2.5 — 3.0. Four phonons are visible in most of the scans, and their 
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Figure 4.13: Measured neutron scattering intensity in Q = (ft k 0) reciprocal plane, 
integrated over a wide Gaussian energy window E = 5.0 ± 1.8 meV, in ZrW 2 8 . Mea- 
surement was done at room temperature, on the BT-9 at NCNR, using fixed incident energy 
Ei = 30 meV, collimation 40'-44'-52'-open, and PG filter before the monochromator. Com- 
pare with Fig. 4.14. 
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Figure 4.14: Calculated phonon intensity in Q = (h k 0) reciprocal plane, integrated over 
a wide Gaussian energy window E = 5.0± 1.8 meV, in ZrW 2 8 . The Figure is comparable 
to Fig. 4.13, except that the Q resolution convolution was not applied in the calculation. 
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Figure 4.15: Neutron scattering intensity in Q = (h k 0) reciprocal plane, integrated over 
E = 6.0 ± 2.0 meV Gaussian energy window, in Sc 2 (W0 4 )2. Measurement was done at 
room temperature, on the BT-2 at NCNR, using fixed incident energy Ef = 30 meV, colli- 
mation 60'-60'-60'-open, and PG filter before the monochromator. Compare with Fig. 4.16. 
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Figure 4.16: Calculated phonon intensity in Q = (h k 0) reciprocal plane, integrated over 
a wide Gaussian energy window E = 6.0±2.2 meV, in ZrW 2 8 . The Figure is comparable 
to Fig. 4.15, except that the Q resolution convolution was not applied in the calculation. 
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energies were extracted by Gaussian fits to the data. By performing such constant-Q scans 
at various points in the reciprocal space, we mapped the dispersion relations (with ener- 
gies up to 15 meV) along: (i) Q = (h 8 0), h = 2—3; (ii) Q = (h 1 0), h = 6—7; (Hi) 
Q = (6 k 0), h = 0—1; (iv) Q = (2 k 0), h = 8—9; and (v) Q = (4 I), h = 2—3. These 
are shown in Figs. 4.18 to 4.22. It is apparent that the ~ 5 meV peaks in the phonon DOS, 
that are related to negative thermal expansion, can be associated with the lowest energy 
optic modes in the system. 

Temperature dependence of the lattice parameters in Sc 2 (W0 4 ) 3 

By using single crystal elastic neutron scattering, we measured the thermal expansion 
in Sc 2 (W0 4 ) 3 along all three axes; the results are shown in Fig. 4.23. We found negative 
thermal expansion along the a and b axes, and positive thermal expansion along the c axis. 
As with the zirconium tungstate, we find the thermal expansion coefficient to be nearly 
constant for temperatures above 50 K, suggesting a mechanism driven by phonons with 
energies of about 5 meV and higher. 



4.3 Analysis 
Model 

To connect static and dynamic measurements, we used a crystal simulation package 
called General Utility Lattice Program (GULP, v. 1.2) [46, 47]. Given phenomenological 
interatomic potentials, GULP can relate physical properties, such as the phonon dispersion 
relations and the density of states, as well as the thermal expansion. Phenomenological 
potentials are extracted by fitting to experimental results for the (i) structure, (ii) phonon 
density of states, and (Hi) the dispersion relations. 

We use two- and three-body interatomic potentials, long-range Coulomb interaction, 
and the shell model for oxygen ions. GULP offers a number of different two- and three- 
body potentials. For two-body potentials, we adopted a so-called Buckingham potential: 

C 

U huck (r) = Aexv(-r/p) - — . (4.5) 

Here, A is the strength, and p the range, of the exponential, repulsive force between atoms, 
while C is the strength of the attractive, power-law force. We used Sc-O, W-O, and 0-0 
two-body potentials. Among the available three-body interatomic potentials, we selected a 
harmonic potential: 

u thKe = ^{e-e ) 2 . (4.6) 

Here, k is the strength of the interaction, 9 is the angle between the two metal-oxygen 
bonds, and 9 is the equilibrium angle (for Sc0 6 octahedra 9 = 90°, while for W0 4 tetra- 
hedra 9 = 104.7°). Three body potentials were used for O-Sc-0 and O-W-0 triplets. 
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Figure 4.17: Constant-Q scans for Sc 2 (W0 4 ) 2 , along the Q = (h 8 0) line in the reciprocal 
space, (a) h = 2.5, (b) h = 2.6, (c) h = 2.7, (d) h = 2.8, (e) h = 2.9, (f) h = 3.0. Lines 
are Gaussian fits. Approximately four phonons are visible in each scan. The energies of 
those phonons are marked in Fig. 4.27. By performing constant-Q scans at various points 
in the reciprocal space, phonon dispersion relations are mapped. Data obtained using BT-2 
thermal neutron triple axis spectrometer at NCNR. 
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Figure 4.20: Measured phonon dispersion relations in Sc 2 (W0 4 ) 2 , along the Q = (6 k 0) 
line. Adapted from [45]. 
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Q = (2 k 0) (r.l.u.) 

Figure 4.21: Measured phonon dispersion relations in Sc 2 (W0 4 ) 2 , along the Q = (2 k 0) 
line. Adapted from [45]. 
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Figure 4.22: Measured phonon dispersion relations in Sc 2 (W0 4 ) 2 , along the Q = (4 I) 
line. Adapted from [45]. 
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Figure 4.23: Measured (circles) and calculated (lines) thermal expansion in Sc 2 (W04) 2 . 
The measured data was obtained by elastic neutron scattering, using SPINS spectrometer 
at NCNR. The uncertainties are comparable to the symbol size. Lines show a thermal ex- 
pansion for a model based on phenomenological interatomic potentials. See also Fig. 4.30. 
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The shell model, used for oxygen ions, splits an ion into two charged entities, a massive 
core (nucleus), and a massless electronic shell [48]. The two are connected by a harmonic 
(spring) interaction, of adjustable strength. Only the shell interacts with neighboring ions. 
All the ions are charged (oxygen core and shell have different effective charges), and inter- 
act via long-range electrostatic (Coulomb) interaction. 



Structure fit 

To extract interatomic potentials, we first performed a structure fit. This consists of 
selecting a set of interatomic potentials, using GULP to relax the structure (i.e. minimize 
the potential energy by varying unit cell size and atomic positions), and then comparing 
the relaxed and the experimentally determined structure. The comparison is quantitatively 
done with the following function: 

X 2 ({Pi}) = 5^|(dex P )|-(d re i)/ + 

const, x ((a exp - a re i) 2 + (6 exp - b rei ) 2 + (c exp - c re i) 2 ) . (4.7) 

Here, {p^ refers to all the interatomic parameters, such as A, p, and C in the Buckingham 
potential, and k and 6> in the three-body harmonic potential. Next, d are the fractional 
atomic positions within the unit cell, and a, b, c unit cell dimensions; subscripts "rel" 
and "exp" refer to relaxed and experimentally determined values. The constant const. 
is adjusted as needed during the fit, to achieve acceptable agreement of all the structural 
parameters. During the structure fit, we minimize x 2 by varying interatomic potentials 

{Pi}- 

Using the structure fit, we were able to extract the two-body potential parameters, as the 
X 2 has a well-defined minimum as a function of those; see Fig. 4.24. As a function of three 
body potential parameters, on the other hand, x 2 has a very broad, nearly flat minimum, 
as shown in Fig. 4.24. Only a lower cutoff value for k could thus be determined from the 
structure fit. 



Phonon density-of-states fit 

We used the measured density of states to further refine the phenomenological potential 
parameters. The density of states measured in a neutron scattering experiment is not the ac- 
tual phonon density-of-states (DOS). Phonons generally have different scattering structure 
factors, and accordingly different weight in the measured spectrum. We refer to the mea- 
sured density of states as "neutron weighted density-of-states" (NDOS). GULP calculates 
actual phonon density of states, but can also calculate the phonon dynamical matrix. By 
diagonalizing the dynamical matrix, the phonon eigenvectors can be calculated, and from 
these the scattering structure factor, and the weight of each phonon mode in a measured 
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Figure 4.24: Quality of the structure fit as a function of the strength of the repulsive com- 
ponent of the two-body interatomic potential between scandium and oxygen ions (A Sc -o), 
inSc 2 (W0 4 ) 2 . 
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Figure 4.25: Quality of the structure fit as a function of the strength of the three body inter- 
atomic potential, involving one scandium and two oxygen ions (A; -sc-oX in Sc 2 (W0 4 ) 2 . 
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spectrum can be determined. In this manner, we were able to calculate the NDOS, and 
compare it with the measured spectrum. We found the NDOS to be strongly dependent on 
the parameters of the three-body potentials. We also found that the fit to the three body 
potentials was not unique. Rather than a point, the minimum in the error function consisted 
of large valleys in the parameter space. Figure 4. 1 1 shows the measured and calculated 
NDOS, for a particular choice of the three-body parameters. 

Phonon dispersion relations analysis 

The next step in the refinement consisted of a fit to the dispersion relations. Having 68 
atoms per unit cell, Sc 2 (W0 4 ) 3 has 204 phonon modes. Density-of- states measurements 
suggest that there are no phonons with energies above approximately 150 meV. This leads 
to a large density of modes, more than one mode per meV. Analysis suggests that 144 of 
these modes have energy below 70 meV (i.e. below the gap separating the high and low 
energy parts of the spectrum), and this implies an even higher density at low energies. 
Most of these modes are not observed in an experiment, due to the small structure factor, 
and consequently the low intensity. For illustration, figure 4.26 shows dispersion relations 
of all the calculated phonon modes, over-plotted with the experimentally measured ones. 
On the other hand, if only the most intense modes are considered, the number of calculated 
modes drops substantially, and is comparable to the number of observed modes. This is 
illustrated in Fig. 4.27. 

Regrettably, we found that we could not adjust the small number of the three-body 
parameters so as to reproduce the energies and intensities of all the observed phonon modes. 
Since the parameter space valleys (corresponding to the best structure and NDOS fits) are 
very large, disconnected, and not well-charted, it is plausible that our initial parameters 
were too far off to yield a good fit to the dispersion relations, by incremental improvements 
in the fit. A better set of initial interatomic potentials might resolve this problem. Such a 
set of potentials could possibly be obtained from the first-principles calculations [49]. 

To investigate the wave-vector dependence of the phonon modes calculated from the 
phenomenological interatomic potentials, we calculated the neutron scattering intensity of 
those phonon modes in the Q = (h k 0) reciprocal plane, at energies around 5 meV. 
Figure 4.13 shows calculated phonon intensities in ZrW 2 8 , over a wide Gaussian energy 
window E = 5.0 ± 1.8 meV, while the Fig. 4.14 shows the corresponding neutron scatter- 
ing measurement. Figures 4.16 and 4.15 show the equivalent calculation and measurement 
for Sc 2 (W0 4 ) 3 , with Gaussian energy window E = 6.0 ± 2.0 meV. When comparing the 
corresponding figures, one should keep in mind that the Q resolution convolution was not 
applied in calculation, to preserve the higher level of detail. We notice significant similar- 
ities in the wavevector dependence of calculated and observed intensities. This suggests 
that the calculated phonon modes resemble those that are observed experimentally. Figure 
4.28 shows four frames from an animation depicting one such low-energy phonon mode. 
Notice that the polyhedra fold, but do not distort, during the phonon cycle. In contrast, 
Fig. 4.29 shows a high-energy phonon mode, in which a large distortion of the polyhedra 
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Figure 4.26: Phonon dispersion relations along the Q = (h 8 0) direction, in Sc 2 (W04) 2 . 
Circles (•) represent phonons observed in the neutron scattering experiment. Lines are the 
calculated phonon dispersion relations. Compare with Fig. 4.27. 
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(h 8 0) (r.l.u.) 

Figure 4.27: Phonon dispersion relations along the Q = (h 8 0) direction, in Sc 2 (W0 4 ) 2 . 
Circles (•) represent phonons observed in the neutron scattering experiment. Lines are the 
most intense of the calculated dispersion relations. Line thickness is proportional to the 
intensity. Compare with Fig. 4.26. 
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is noticeable. 

Thermal expansion 

Figure 4.23 shows the thermal expansion in Sc 2 (W0 4 ) 3 . Circles represent the mea- 
surements obtained by elastic neutron scattering, and lines the calculation based on the 
interatomic potentials derived from the fitting to the observed structure and the phonon 
density of states, as previously described. The thermal expansion was calculated by mini- 
mizing the free energy (rather than the internal energy), as a function of unit cell size, and 
atomic positions. We notice that the current potentials cannot predict the negative ther- 
mal expansion. Still, by adjusting the three body potentials, we are able to dramatically 
change the thermal expansion, and even quantitatively reproduce the NTE along some of 
the axes. The density-of- states calculated after such adjustments, however, significantly 
deviates from the measured DOS, indicating that adjusted potentials are still inaccurate. 
We ascribe the fact that satisfactory interatomic potential parameters could not be found 
in our work to the complexity of Sc 2 (W0 4 ) 3 , rather than to a deficiency in fundamental 
aspects of the theoretical description of the material. In the simpler material ZrW 2 8 , for 
example, Figure 4.30 shows a good agreement between the measured and the calculated 
thermal expansion. 

4.4 Conclusion 

Our research indicates that the three-body interatomic potentials play an important role 
in the thermal expansion of NTE systems. In contrast, two-body potentials are largely 
responsible for the crystal structure, and do not significantly affect the thermal expansion 
properties. Furthermore, two body potentials provide the main contribution to the total 
lattice energy. This is consistent with the conjecture that the system can be represented as a 
network of semi-rigid units: cation-oxide polyhedra. Since the cation-oxide bonds are very 
strong, rigid units have very little (positive) thermal expansion. If the network of rigid units 
is sufficiently open (i.e. the number of links between them is small), and has an appropriate 
topology, rotations of these rigid units costs very little energy. In certain phonon modes, 
called rigid-unit-modes (RUM), such rotations play an important role, and can cause the 
negative thermal expansion. Although our experiments did not provide conclusive evidence 
regarding the existence of RUMs (for instance, an experimental observation of a rigid-unit 
structure factor in the dynamics spectrum), the results were consistent with the theory of 
the negative thermal expansion caused by the rigid-unit-modes. 
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Figure 4.28: Four frames from an animation depicting a phonon mode in Sc 2 (W0 4 ) 3 . The 
eigenvectors were calculated from a model based on phenomenological interatomic poten- 
tials. Shown is a top layer of the b-c plane. The phonon shown is Q = (1/4 0) phonon 
mode with E = 5.3 meV (mode #4). The phonon frequency is / = 1.3 THz, and the 
frames correspond to the following phase angles: (a) <f> = 0, (b) <$> = ir/2, (c) <f> — n, and 
(d) cj) = 37r/2. Shown in blue are W0 4 tetrahedra, and in green Sc0 6 octahedra. This a 
rigid-unit-mode; it consists of folding of the polyhedra, with negligible distortion. Such 
modes cause the negative thermal expansion in Sc 2 (W0 4 ) 3 . Compare with Fig. 4.29. 
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Figure 4.29: Four frames from an animation depicting a phonon mode in Sc 2 (W04) 3 . The 
eigenvectors were calculated from a model based on phenomenological interatomic poten- 
tials. Shown is a top layer of the b-c plane. The phonon shown is Q = (1/4 0) phonon 
mode, with E = 42.1 meV (mode #1 11). The phonon frequency is / = 10.2 THz, and the 
frames correspond to the following phase angles: (a) <f> = 0, (b) <j> = 7r/2, (c) (p = tt, and 
(d) (j) = 3ir/2. Shown in blue are WO4 tetrahedra, and in green Sc0 6 octahedra. This is 
not a rigid-unit-mode; notice a large distortion of the polyhedra during the cycle. Compare 
with Fig. 4.28. 
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Figure 4.30: Measured (symbols) and calculated (line) negative thermal expansion in 
ZrW 2 8 . Model, based on phenomenological interatomic potentials, predicts a phase tran- 
sition around T »s 350 K. Zirconium tungstate has a second order structural phase transi- 
tion at T = 428 K. Bullets (•) adapted from Ref. [39] (uncertainties are comparable to the 
symbol size); open circles (o) adapted from Ref. [50]. 
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Chapter 5 

Triangular Antiferromagnet 
RbFe(Mo0 4 ) 2 

5.1 Introduction 

During 1970s, a group of inorganic chemists embarked on a enduring task of systematic 
synthesis and study of all possible compositions of double molybdates, M + R 3+ (Mo04) 2 
[51, 52]. A family of fifteen compounds, labeled KAl(Mo04) 2 after the first compound for 
which the structure has been determined, was found to crystallize in a layered, triangular 
structure, space group P3ml. Two decades later, a member of this chemical family, ru- 
bidium iron molybdate, will pick up the interest of condensed matter physicists seeking a 
good physical realization of a two dimensional triangular antiferromagnet (2D TAF). 

In rubidium iron molybdate (RbFe(Mo0 4 ) 2 ), magnetic ions, Fe 3+ , form a simple trian- 
gular lattice, and are connected antiferromagnetically via O-Mo-0 superexchange paths. 
The unit cell contains one formula unit, thus the Fe 3+ -Fe 3+ distance within one layer equals 
the a length of the unit cell, a = 5.7 A, As the spin of iron ions is rather large, S = 5/2, 
the system can best be described as quasiclassical. 

5.2 X-ray diffraction 

To determine the structure of our sample, we sent a few small crystals to X-ray Crys- 
tallographic Laboratory at the University of Minnesota. The structure report is reproduced 
in the Appendix B.l (p. 139), while here we highlight a few relevant results. Structure 
determination has been carried out on a single crystal of approximate volume 10~ 4 mm 2 , 
at T = 173(2) K. Structure analysis was based on 1272 reflections, and confirmed that our 
sample crystalizes in a trigonal space group P3ml (as originally proposed), with lattice 
parameters a = 5.6923(6) A, and c = 7.4572(11) A. 

The structure can be described as follows. There is one iron ion per unit cell, and these 
ions reside on a perfect, stacked, triangular lattices. Each iron ion is surrounded by six 
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Figure 5.1: Structure of RbFe(Mo0 4 ) 2 . Shown are iron ions in red, oxygen ions in blue, 
Mo0 4 tetrahedra in green, and rubidium ions in pink. Unit cell is outlined, (a) Basal plane. 
Notice that iron ions lie on a triangular lattice, (b) a-c plane. Notice the lack of formal 
covalent bonds between layers. 
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oxygen atoms in octahedral coordination. Molybdenum atoms lie above and below the 
triangular lattice, and are surrounded by four oxygens in a tetrahedral configuration, three 
of which are shared by the Fe0 6 octahedra. Thus, each iron ion is connected to its six 
in-plane iron neighbors through the two Fe-O-Mo-O-Fe superexchange paths. The fourth 
oxygen in the M0O4 tetrahedra is unshared; it is only bonded to the central molybdenum 
atom, and located nearly between the planes. The rubidium atoms are not bonded, and 
reside exactly between the planes comprised of iron atoms. Covalent bonds, therefore, 
exist only within the planes, and no superexchange paths connect interplane ions. The ratio 
of interplane to intraplane distances between magnetic ions is c/a = 1.310. 

Looking only at the magnetic structure (ions and exchange paths), we see stacked planes 
containing edge sharing equilateral triangles, the vertices of which are occupied by S = 5/2 
spins, the sides representing interactions between them. If those interactions are antiferro- 
magnetic, our sample would be a good candidate for study of frustrated magnetism on a 
triangular lattice. 

5.3 Magnetic susceptibility 

DC magnetic susceptibility measurements on RbFe(Mo0 4 ) 2 were performed using a 
Quantum Design MPMS-5S SQUID magnetometer. The experiments were performed in 
two configurations: (a) with an in-plane magnetic field (denoted ||), where a single crystal 
sample of mass my = 4.119 mg was used; and (6) with a field perpendicular to the plane 
(denoted _L), with m ± = 4.495 mg. In both cases, a DC field of H = 100 Oe (H = 
7.96 x 10 3 A/m = 0.01 B) was used. The samples were mounted with a string, in the 
proper orientation, inside a non-magnetic, plastic tube weighting about 22 mg. 

The measurements were conducted over the entire accessible temperature range, from 
room temperature down to T = 1.8 K. The step size varied from 20 K at room temperature, 
to 0.05 K at the base temperature. The susceptibility was found to be a smooth function of 
temperature over the entire range. Measurements were done on both cooling and heating, 
and no hysteresis was observed. 

Looking at the raw data, a phase transition, signaled by a prominent cusp occurring at 
the same temperature for both crystal orientations, was observed at T N = 3.55(5) K. Below 
the transition, the susceptibility drops slightly (< 5%). This temperature dependence of the 
magnetic susceptibility indicates that our sample is indeed an antiferromagnet. In param- 
agnets, the susceptibility increases with cooling, and diverges as T — >■ 0. In ferromagnets, 
the divergence occurs at some finite temperature, called Curie temperature (T c > 0), below 
which a long range ferromagnetic order is established. In contrast to these, susceptibility in 
antiferromagnets attains a maximum value (but does not diverge) at the Neel temperature, 
T N , below which antiferromagnetic order is established. In conventional antiferromagnets, 
this order is characterized by an antiparallel orientation of neighboring spins, whereas in 
frustrated magnets no such configuration exists. A clue that our system is not an ordinary 
antiferromagnet is provided by the dependence of the in-plane susceptibility with temper- 
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Figure 5.2: Magnetic susceptibility in RbFe(Mo0 4 ) 2 , with magnetic field applied in the 
basal plane (A), and perpendicular to the basal plane (o). Inset: inverse magnetic suscepti- 
bility. See also Fig. 5.3. 
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Figure 5.3: Magnetic susceptibility in RbFe(Mo0 4 ) 2 around the T = 3.5 K phase tran- 
sition, with magnetic field applied in the basal plane (A), and perpendicular to the basal 
plane (o). 
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ature below the transition: as T —> it flattens at a value that is only about 3% below 
maximum. This is in sharp contrast to ordinary antiferromagnets, where X\\(T — > 0) — > 0. 
Thus far, without any numerical analysis, we have determined that our system has anti- 
ferromagnetic interactions that cannot all be satisfied by a static long range ordered spin 
configuration. 

We can use the susceptibility data to extract a quantitative information about the strength 
of the interaction between spins. The nearest neighbor exchange constant, J, is related to 
the Curie- Weiss temperature by: 

6^ = ^5(5 + 1), (5.1) 

where z = 6 is the coordination number for iron ions, and Qcw me Curie-Weiss tempera- 
ture. This is related to our measurement via Curie-Weiss law, which gives the temperature 
dependence (above the transition) of the magnetic susceptibility: 

C 

l + |fc>cw| 

Fitting our data to the Curie-Weiss law gives 6^ = -20.8(5) K, and Q^ w = -24.7(3) K. 

[Equation (5.2) is valid only for T ^> T N , and deviations — although minute — do 
show up as T — > T N . Hence, to extract ®cw, on ly the data above some cutoff temperature, 
T cut , should be used in a fit. The results just reported were obtained with T cut = 20 K. 
Different cutoff temperatures produce different Qcw, although the results agree to within 
1 K for 5 K< T cut < 50 K. Lower cutoff temperatures tend to distort the result due to the 
proximity of the transition, while higher ones lead to increased errors as the number of data 
points gets lower, and the relative uncertainty in susceptibility is higher. After examination, 
we arrived at T cut = 20 K as a good compromise.] 

Equation (5.1) now yields the strength of the exchange interaction, J m 0.1 meV. The 
values for in-plane and perpendicular orientations are: jH = 0.104(3) meV, and J L = 
0.124(2) meV, indicating a planar anisotropy. 

Finally, susceptibility measurement enables us to quantify the degree of frustration in 
our system. The ratio of the Curie- Weiss and Neel temperatures has been proposed as a 
quantitative measure of frustration, and here it is / = \®cw\/Tn ~ 6. This is quite a low 
value, as for some triangular antiferromagnets the ratio is many times higher. In NaTi0 2 , 
for instance, / > 500 [53]. Why is this so? One of the possible explanations is the high spin 
value, S = 5/2. Low spin values suppress ordering, as the spins remain in fluctuating state 
due to the quantum effects; high spins have the opposite effect. Another important factor 
is interplane interactions. Even though the structural analysis does not clearly identify any 
interplane superexchange paths, the small ratio of interplane to intraplane spin distances 
(c/a = 1.310) may enhance ordering by effectively increasing the dimensionality. Bulk 
measurements, however, are not a powerful enough probe to address these questions. 



77 



5.4 Neutron powder diffraction 

Faced with a new material, about which little is known, an experimenter must design a 
series of experiments which will enable him or her to study the interesting physics in the 
most efficient manner. One usually starts with bulk experiments, such as our susceptibility 
measurements, for which preparations are minimal, and based on the results experiments 
of increasing complexity are designed and executed. Ultimately, we wanted to use neu- 
tron scattering on single crystals — one of the most detailed probes known — to study 
the magnetic properties. However, the reciprocal space, in which single crystal scattering 
experiments take place, is a vast, four-dimensional space. Without signposts along the way, 
one can easily get lost there, just as in our direct space. An important sign in our journey 
was provided by powder diffraction experiments. 

Though not as powerful as single crystal scattering, powder diffraction experiments 
offer a panoramic view in which interesting features are easier to notice. Being isotropic, 
powders reduce dimensionality of the reciprocal space to two. Furthermore, for elastic 
scattering the space is reduced to a line — the scattering angle — thus enabling the study 
of all static features in a single sweep. This reduction in dimensionality causes much of the 
information about features to be lost. Nonetheless, a powder experiment can show where 
to look in a more detailed experiments. In our case, the powder experiment told us where 
to look for magnetic Bragg peaks, the study of which will help us decipher the magnetic 
structure. 

We conducted two powder diffraction experiments, both using the BT-1 High Resolu- 
tion Powder Diffractometer (HRPD) at the NIST Center for Neutron Scattering. For the 
first experiment, the powder sample was prepared by crushing all small remaining crys- 
tallites after the thirty-four biggest crystals were isolated. The total mass of the sample 
was m p i = 2.443 g. The sample was loaded into a cylindrical vanadium container with 
diameter d = 6 mm. Although less transparent to neutrons than aluminum, vanadium is 
commonly used for containers in powder diffraction experiments. As vanadium has a negli- 
gible coherent scattering cross section, there is no contamination of the signal by vanadium 
container Bragg peaks. The incoherent cross section is significant for vanadium; however, 
it produces a uniform, angle independent, background, which is easily subtracted. For this 
first experiment, we used a Ge(3 11) monochromator; with the scattering angle of 29 = 75°, 
it reflects neutrons of wavelength A = 2.0784 A, and flux $ = 7 x 10 5 s _1 cm~ 2 . 

Using an ILL "orange" cryostat, we ran measurements at three different temperatures: 
(i) T\ = 1.7 K, 12 hours; (ii) T 2 = 11 K, 12 hours; and (in) room temperature, 9 hours. 
The lowest two temperatures were chosen to bracket the transition temperature; then the 
difference in the spectrum between the two is caused by magnetic scattering. The purpose 
of the high temperature run was comparison with X-ray refinement and literature data. 

The experiment was very successful, as we observed a number of magnetic Bragg 
peaks; we will come to this later. First, let us address the problem that moved us to re- 
quest another sample, and repeat the experiment three months later. Structure refinement, 
initially using the highly symmetric space group P3ml, then those with lower symmetries, 
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did not yield satisfactory results. The quality of the fit was too low. While the refinement 
in space group P3ml showed the best agreement with the data, the profile of the peaks 
clearly deviated from Gaussians. (The instrument resolution of BT-1 spectrometer is such 
that all Bragg peaks have a perfect — within the experimental error — Gaussian profile.) 
Analysis suggested that the sample consisted of more than one phase, and those phases had 
similar structures — too similar to be distinguished in an experiment. 

The conclusion was reached that the culprit was the grinding process. When using a 
mortar and pestle, excessive force was used, along with the shear motion which tended to 
efficiently crush crystals. However, in layered structures such as this, where layers are very 
weakly bonded, structural strain induced by such a treatment can be non-uniform. As a 
possible remedy, it was suggested that heating the sample could possibly relieve the strain. 
Not knowing, however, the melting temperature, nor the chemical or structural stability of 
the sample at high temperatures, this approach was considered too risky. Instead, after the 
experiment, a new polycrystalline sample was requested. 

The polycrystalline sample came in the form of two cylindrical, ceramic pellets, which 
were then carefully crushed: slowly, by gentle tapping, and strictly avoiding the side- 
ways (shear) motion of the pestle. This time, the procedure took several hours. The 
total mass of this sample was m p2 = 4.854 g. The second experiment was done using 
Cu(311) monochromator, with scattering angle 28 = 90°, and neutrons of wavelength 
A = 1.5401 A. In this configuration, the rated flux is $ = 4 x 10 5 s _1 cm~ 2 , while our 
monitor recorded mrat = 7.02 x 10 5 min -1 . Collimations were the same as in the first 
experiment. [The configurations for both experiments were chosen to maximize the flux, at 
the expense of resolution. The resolution width (as a function of scattering angle), however, 
differs somewhat between the two.] 

We used ILL "orange" cryostat (70 mm, OC-8), and vanadium container of diameter 
d = 9 mm. We run: (i) T x = 1.7 K, 16 hours; (ii) T 2 = 10 K, 13.5 hours; and (Hi) 
T 3 fa 175 ± 3 K, 5 hours. This time, while the quality of the structure refinement was good, 
the results were somewhat surprising. It was shown that the sample crystallizes in space 
group P3 (less symmetric than P3ml), in disagreement with the X-ray experiment, as well 
as the existing literature [51]. Knowing, however, that this time our result was reliable, 
we focused on magnetic study. All conclusions about RbFe(Mo0 4 ) 2 being a good trian- 
gular antiferromagnet were unaltered by this finding. A year later, a report was published 
stating that RbFe(Mo0 4 ) 2 undergoes a subtle structural second-order phase transition at 
T = 180 K, corresponding to a small rotation of the Mo0 4 tetrahedra. The authors re- 
port the symmetry above the transition to be P3ml, and P3cl below; thus, corroborating 
our results, and providing an explanation for the initial disagreement. The difference in 
space groups for low-temperature phases is due to different stacking of layers, as will be 
discussed later. 

Structural issues discussed, we can now concentrate on the magnetic structure. The 
difference between the T x = 1.7 K and T 2 = 10 K spectra was due to the magnetic Bragg 
peaks, and we could distinguish — some more clearly than others — twenty-one of them. 
The number and the intensity of those peaks was too low for direct magnetic structure 
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refinement. Instead, we assumed that the the projection of the magnetic unit cell onto the 
basal plane is three times larger than for the nuclear unit cell. This is true for the ground 
state of classical spins on a triangular lattice, which form so-called 120° structure. Then, 
the magnetic Bragg peaks can appear only at Q = (h/3 k/3 I ± A), where h and k are 
integers. (Not all of those peaks need to be present: many will vanish due to the structure 
factor. However, the positions of those that are present must have the indicated form.) This 
left only one unknown parameter: A. We plotted calculated scattering angles, 29 hkA , of all 
(h/3 k/3 I ± A) peaks as a function of A, and compared them with positions of measured 
peaks. There were two close matches. The first one, A = 0.152, accounted for most, 
but not all, peaks. The second match, arrived at days later, accounted for all peaks, and 
eventually proved to be correct: A = 0.453. This was a major result: Now we knew where 
to look for magnetic peaks in the upcoming single crystal experiments. Indeed, we found 
thematQ = (/i/3 k/3 I ± 0.453). 

Before moving onto single crystal experiments, let us discuss another important result 
we extracted from powder diffraction measurements: the ordered moment. This analysis 
was done about a year after the diffraction experiments. By that time we have determined 
the magnetic structure using single crystal scattering, and could now use it in conjunction 
with the powder data. As for the nuclear structure refinement, we used the GSAS pack- 
age [54] for the magnetic data. We represented our sample as two phases: nuclear and 
magnetic. The magnetic phase contained only spins, the magnitude of which was set to 
unity. The ratio of phase fractions (magnetic to nuclear) was then fit to the data, and it 
directly gave the magnitude of the ordered moment. The details of this procedure are given 
in the Appendix. The result revealed the (static) ordered moment to be p = 4.49(5) /j,b, 
or about 75% of the available moment (25% of the moment remains in a fluctuating state). 
This corresponds to an effective spin value (S) = 1.80(2), as opposed to S = 5/2 for full 
polarization. 

5.5 Single crystal elastic measurements 

For single crystal scattering measurements, we used our largest single crystal (m^ = 
35.2 mg), mounted so that the horizontal plane coincided with the (h h I) reciprocal plane. 
Having determined the most probable locations of the magnetic Bragg peaks from the pow- 
der diffraction measurements, we were able to quickly locate them in the subsequent single 
crystal experiment. Below the transition, magnetic peaks were well defined and intense, 
and vanished right at the transition, thus confirming their magnetic origin. We found a 
number of magnetic peaks, located at Q = (h/3 h/3 I ± A). 

Using the BT-2 TAS spectrometer, we performed high resolution Q- scans through mag- 
netic Bragg peaks, along both the h and I directions, from the lowest temperatures up to the 
paramagnetic phase (to T = 3.65 K), to determine the temperature dependence of the posi- 
tion and width of those peaks. We found the width of the peak to be resolution limited at all 
temperatures (FWHM for Q = (1/3 1/3 0.453): S t = 0.008 r.l.u., and 5 h = 0.005 r.l.u.). 



83 



oj 4000 

t— 1 

3000 



2000 

Oh 



o 
u 



1000 




4.0 



Figure 5.7: Intensity of the magnetic Q = (1/3 1/3 0.453) Bragg peak in RbFe(Mo0 4 ) 2 . 
Intensity is proportional to the order parameter. Fit to M ~ (T c — T) Q yields critical 
exponent a = 0.67(2). 



The position of the peak was also found temperature independent. 

The intensity of magnetic Bragg peaks is proportional to the staggered magnetization, 
and can be used as an order parameter for this system. To extract the critical exponent char- 
acterizing the phase transition, we measured the temperature dependence of the integrated 
intensity of the Q = (1/3 1/3 0.453) magnetic peak. Below T = 3.0 K, we measured 
the intensity with AT = 0.1 K temperature step; AT = 0.01 K for 3.00 K< T < 3.40 K; 
AT = 0.005 K for 3.40 K< T < 3.65 K; and AT = 0.05 K for temperatures above 
T = 3.65 K. Close to transition, we counted for 290 s per temperature point, reaching 
about 13000 counts at T = 3.0 K. We fit the intensities (for T < T c ) to: 



M oc {T c -Tf, 



(5.3) 



and extracted the critical exponent (3 = 0.29(1), with the transition temperature T c = 
3.41(1) K. 

[Despite good statistics, the extraction of the critical exponent was complicated by the 
fact that the intensity did not drop to zero at the transition temperature; instead, a long tail 
extended to T 3.57 K, up to which the intensity gradually vanished. This additional 
intensity may be due to quasi-elastic critical scattering. We tried to measure and subtract 
the critical scattering, by measuring the intensity at the tail of the elastic magnetic peak, 
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at Q = (1/3 1/3 0.4415), counting 6 minutes per point, at temperatures 3.40 K< T < 
3.60 K, with temperature step AT = 0.01 K. Unfortunately, this scan did not return a 
critical scattering peak. A possible explanation is that mechanical strain or disorder in 
the sample may cause domains of inhomogeneity within which the transition temperature 
varies slightly, thus causing the observed tails close to the transition. We tried to account 
for the tail by adding an empirical term: 



cnt - |( T _ Tc )2 + («5T)2| 7 /2' ^ " ' 

but the fit wasn't satisfactory. Eventually, we extracted the critical exponent by using only 
intensities below some cutoff temperature, and adjusting the transition temperature to ob- 
tain the best fit to the data.] 

To determine the magnetic structure of the ordered phase, we measured the intensity 
of a large number of magnetic Bragg peaks. This experiment was done on the BT-2 spec- 
trometer, at base temperature (T = 1.7 K) in two configurations: (i) using E &x = 30 meV 
neutrons, where intensities of 114 magnetic, and 21 nuclear, peaks were measured; and 
(ii) E fix = 14.7 meV, with 56 magnetic and 12 nuclear peaks. The scattering plane was 
(h h I), and magnetic peaks were observed atQ = (h/3 h/3 I ± A), where h and I 
are integers, and A = 0.453. (It was, however, found that all magnetic peaks of the form 
Q = (h h I ± A) had zero intensity.) Nuclear Bragg peaks were scanned for normalization 
of magnetic intensities into absolute units. 

In order to determine the magnetic structure of the ordered phase, we reasoned as fol- 
lows: As the measurements were done well below the Curie-Weiss temperature (|6cw| ~ 
22 K), we can assume that the spins on each frustrated unit, a triangle, are in a minimum 
energy configuration. Spins 5 = 5/2 can approximately be treated as vectors (classical 
approach), and the minimum energy is achieved when the total spin on a triangle vanishes. 
Therefore, spins on a triangle must be coplanar, and rotated by an angle of 120° relative to 
each other. Since the triangular lattice consists of strongly connected, edge sharing, trian- 
gles, the order on one triangle will constrain the spins in the entire triangular plane, into 
a so-called "120°" spin structure. Such magnetic structure is characterized by a magnetic 
cell which is three times larger than the nuclear cell (within the plane). The ground state 
is three-fold degenerate, stemming from the three degrees of freedom in which the spins 
on the first triangle can be oriented. Two degrees of freedom determine the plane in which 
spins reside; the third an arbitrary rotation within the plane. Thus, our approach to mag- 
netic structure analysis was to determine (i) whether the intensity of the magnetic peaks is 
consistent with the 120° structure, (ii) if so, the orientation of the plane in which the spins 
reside, and (Hi) the value of the static magnetic moment. 

Clues for the first two tasks can be gained just by a qualitative analysis of the intensities. 
We have already remarked the absence of all the magnetic peaks with the general location 
Q = (h h I ± A). This is consistent with the magnetic form factor for the 120° structure, 
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Figure 5.8: Ground state for a frustrated antiferromagnet on a triangular lattice: 120° spin 
configuration. Shaded region (a) represents the magnetic unit cell, and shaded region (b) 
the nuclear unit cell. 



which is proportional to: 

4feOC Y.h'v i 5 ( h ' + h + \) x S ( k§ + k + \) 

+6(h' + h-\)x5(k' + k-\)]. (5.5) 

This suggests the affirmative answer to the first question. The second questions can be 
answered by analyzing the angular dependence of the intensities of peaks in the scattering 
plane. If the spins lie in the x-y plane, the angular dependence of the intensity would be 
proportional to: 

1(9) oc {1 + Ql) = (l + cos 2 0), (5.6) 

where 9 is the angle between the scattering wave-vector Q and the 2-axis. We see that the 
measured intensity drops faster along the (h h 0) direction than along the (0 0/). This 
suggests that the spins lie in the basal plane, perpendicular to c direction. The third ques- 
tion, about the magnitude of the static moment, cannot be answered without quantitative 
analysis. 

The details of the quantitative analysis of the intensities of the peaks are given in Ap- 
pendix (Sec. B.4, p. 151). In this section we report the results. First, we wanted to an- 
alyze the intensities of nuclear peaks, for comparison and normalization. We calculated 
the nuclear structure factors using Eq. (3.7) on p. 19 (the fractional coordinates, d^, were 
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determined from the powder structure refinement). Surprisingly, we found significant dis- 
crepancies between the observed and calculated intensities, that could not be accounted by 
small variations in atomic positions, nor extinction due to the shape of the sample. Unable 
to explain the discrepancy despite considerable effort, we reluctantly moved to analysis of 
magnetic peaks. To facilitate the analysis, we analytically calculated the structure factor for 
the stacked, helical and incommensurate along c-axis, 120° magnetic structure, and found 
it to be proportional to: 

(4«) calc <x Y,wwv [8(ti + h + ±)x8(k' + k + ±)x8(l' + l + A) 

+S(ti + h--)x S(k' + k - -) x 8(1' + / - A)] . (5.7) 
3 3 

This result also disagreed with the experiment; we observed twice as many peaks: 

(Mobs* T, W m> [8(ti + h + ±)x8(k' + k + ±)x8(l' + l + A) 

+8(ti + h + \) x 8{k' + k + \)x 8(1' + / - A) 

+5(ti + h--)x 8(k' + k - -) x 6(1' + 1 + A) 
3 3 

+5(ti + h-\)x 8(k' + k-\)x 8(1' + I - A)] . (5.8) 
After plotting the two functions, the solution became apparent: 

(I hhl) obs ~ X ((^w)calc + (I hhl) calc^j ■ (5.9) 

We tested (5.9) against the nuclear peaks, and found it to be in excellent agreement. The 
explanation, already hinted at earlier (p. 81), was that the stacking of layers is not uniform 
throughout the crystal. To better understand this, let us first analyze the high temperature 
phase. We have already mentioned that at T = 180 K, there is a structural phase transition. 
Above the transition, RbFe(Mo0 4 ) 2 crystallizes in P3ml space group. This space group 
has a [1 10] mirror plane, which bisects the Mo0 4 tetrahedra: Mo, 02, and one of 01 atoms 
sit on the mirror plane, while the remaining two 01 atoms are the mirror image of each 
other. In this phase, stacking of layers is uniform throughout the crystal. The structural 
transition consists of a slight rotation of Mo0 4 tetrahedra (approximately ip = 10.8°); all 
tetrahedra within a single domain rotate in the same direction (around c-axis). Thus, the 
crystal goes from high symmetry, into a lower symmetry P3 phase (no mirror plane). It 
appears, however, that the whole crystal does not belong to a single domain. Instead, in 
some domains tetrahedra rotate in one direction (say clockwise), and in others the opposite 
(counterclockwise). This occurs only in the low temperature phase (T < 180 K); heating 
above the structural transition restores the crystal into a single, high-symmetry domain. 
Generally, the total fraction of domains in which all tetrahedra rotated in a certain direction, 
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p^, will not be exactly equal to one-half; i.e. p^ ^ 1/2 ^ p^. (Naturally, p^ = 1.) 

To determine p^, we modified Eq. (5.9): 

(4w)obs ~ (p-r, X (4ft;)calc + * (4w)calc) , (5-10) 

and fitted fractions to measured intensities. We found that to within a few percent p^ 
1/2 p^. Thus, we adopted Eq. (5.9), rather than Eq. (5.10), for calculation. 

For the purpose of structure determination, we derived the analytical structure factor 
and scattering cross section for the 120° structure, helical and incommensurate along the 
c-axis. The derivation is somewhat tedious; the main steps are indicated in Appendix B.4. 
Here we simply present the result: 

x [5{ti + h + \) x 8(k' + k + \) x 8(1' + / + A) 

+6(ti + h--)x S(k' + k - -) x <J(Z' + / - A)] . (5.11) 
3 3 

The only unknown in this equation is the value of the static moment, (S), and so we fitted 
it to magnetic intensities (which were normalized by nuclear intensities). Using E &x = 
30 meV data, we obtained: 

(S) = 1.83(11), (5.12) 

in excellent agreement with the powder diffraction result, (S) = 1.80(2). The larger uncer- 
tainty in the result derived from single crystal measurements may be caused by the shape of 
the sample, the effect of which is an anisotropic extinction that cannot be easily accounted 
for in calculation. 

The data indicated that the magnetic structure is indeed the helical incommensurate 
120° structure, as proposed earlier. The quality of the fit is best measured by the "R factor": 



R = 





(-fj)calc _ (1 


«)obs | 


2^i=i 


(4)obs 





(5.13) 



which is R nuc = 14.8% for the nuclear structure, and R mag = 18.0% for the magnetic struc- 
ture. (R factors below 10% generally indicate well-refined structures. Sometimes, however, 
R factor cannot be reduced to 10% because of disorder, extinction, or other factors. In those 
cases, R factors below 20% may be considered acceptable. Our RbFe(Mo0 4 )2 sample has a 
very anisotropic, plate-like shape, thus making extinction effects significant. Accordingly, 
we considered our result good.) We must note, however, that it is not possible to uniquely 
determine the magnetic structure by such analysis. The Eqs. (3.18) and (3.19) cannot be 
solved in reverse; i.e. given (dcr/dfi), one cannot uniquely determine (Sf). Therefore, we 
can at most say that the measured peak intensities are consistent with the helical incom- 
mensurate 120° magnetic structure. 
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Figure 5.9: Observed and calculated nuclear Bragg peaks in (h h I) reciprocal plane in 
RbFe(Mo0 4 )2. Top (filled) semi-circles represent the measured, and lower (hollow) the 
calculated intensities. Intensity is proportional to the area of semi-circle. The uncertainty 
in the measured intensities is too small to be visible in the figure. 
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Figure 5.10: Observed and calculated magnetic Bragg peaks in [h h I) reciprocal plane 
in RbFe(Mo0 4 ) 2 . Top (filled) semi-circles represent the measured, and lower (hollow) the 
calculated intensities. Intensity is proportional to the area of semi-circle. The area between 
the hollow and the filled top semi-circle corresponds to the uncertainty in the measurement. 



90 



Finally, knowing the magnetic ordering along the c direction, and that it is probably 
due to an interlayer exchange coupling, we tried to extract the ratio of inter- to intralayer 
exchange constants. The details of this calculation are given in the Appendix B.5; here we 
present only the conclusion. The analysis revealed that at least three interlayer exchange 
constants are needed to account for the incommensurate ordering observed. Let J 2 be 
the exchange constant between nearest interlayer neighbors, and J 3 and J 4 be the exchange 
constants between diagonal interlayer neighbors (there are six such neighbors per layer, and 
they occupy two non-equivalent positions; thus the two diagonal constants; see Fig. 5.1 1). 
We conclude that: 

where <f> = A x 360° Rj 163°. We immediately see that J 3 ^ J 4 . Thus, we have more 
unknowns than equations, and therefore cannot determine the ratio we sought. However, 
the observed structure does place the constraint of Eq. (5.14) on the exchange constants. 



5.6 Magnetic field and inelastic measurements 

In this section, we report our measurements on RbFe(Mo04) 2 single crystal in high 
magnetic field, as well as inelastic neutron scattering experiments. The measurements and 
the analysis reported in this section are a part of an ongoing effort, and we anticipate more 
important results in the future. 



H-T phase diagram 

High magnetic field experiments were carried out on BT-2 and SPINS neutron spec- 
trometers. We mapped the H-T phase diagram, by measuring the intensities and positions 
of magnetic Bragg peaks. [The magnetic field is in-plane. Applying magnetic field perpen- 
dicular to the basal plane does not yield new phases for the following reason: When the field 
is applied in-plane, the minimum energy configuration has all the spins in-plane. Applying 
field perpendicular to the plane causes spins to rotate in the direction of the magnetic field, 
in an "umbrella-like" fashion; in-plane component of the spin remains the same. Thus, new 
phases are driven only by the in-plane component of the applied field.] We discovered a 
rich phase diagram, with (i) commensurate and incommensurate phases; (ii) first and sec- 
ond order phase transitions; and (Hi) low temperature hysteresis. Fig. 5.13 shows the H-T 
phase diagram. Paramagnetic phase, marked by the absence of magnetic Bragg peaks, is 
found at high temperatures, T > 3.5 — 4 K (the exact transition temperature is field de- 
pendent). At lower temperatures, the long-range ordering is field dependent. The low field 
phase, which develops below (j, H = 3 T at T = 3.0 K, is the same incommensurate phase 
that we described in the previous section. In-plane we have a 120° magnetic structure, 
while the ordering along the c-axis is incommensurate (A = 0.453, i.e. nearest interplane 
neighbor spins are rotated by <f> = A x 360° fa 163°). When the magnetic field is increased 
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Figure 5.11: Intraplane and interplane exchange constants in RbFe(Mo0 4 ) 2 . Figure (a) 
shows a perspective view, while (b) shows top view. Exchange constants in Fig. (b) refer 
to interplane interactions. Notice on Fig. (b) that the mirror plane of the unit cell does not 
coincide with the mirror plane of the M0O4 tetrahedra. J is the intraplane nearest neighbor 
exchange constant, J 2 is the interplane nearest neighbor exchange constant, while J 3 and 
J 4 are the diagonal interplane exchange constants. Analysis indicates that J 3 7^ J 4 , and 
that the two exchange constants play an important role in establishing the helical magnetic 
order. Equation 5.14 gives the relation between the diagonal exchange constants. 
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Figure 5.12: Position of magnetic Bragg peaks along Q = (1/3 1/3 I) in RbFe(Mo0 4 ) 2 . 
(a) H = 2.5 T. (b) H = 8.3 T. (c) H = 10 T. See Fig. 5.13 for phase diagram. 



above B = 3 T, the system goes through a first order phase transition into a commensurate 
configuration, with A = 1/3; interplane spin rotation is now = A x 360° = 120°). In- 
creasing the field further above /i H = 9 T, the ordering again becomes incommensurate, 
with A 0.4. This phase transition is second order. Here, unlike at lower magnetic fields, 
the position of magnetic Bragg peaks is not fixed. It is weakly field dependent; we refereed 
to them as sliding peaks. 

We found the measured H-T diagram to be qualitatively very similar to that of a tri- 
angular X-Y antiferromagnet (see Ref. [56]). This enabled us to identify the in-plane spin 
configuration in each of the phases. We know, from the results presented in the previous 
section, that in zero field we have the 120° spin configuration: spins come in three ori- 
entations, rotated by an angle of 120° with respect to each other. As the field is applied, 
the spins slightly rotate in the direction of the field. One of the spins still remains nearly 
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Figure 5.13: Magnetic field vs. temperature phase diagram for RbFe(Mo0 4 ) 2 . Adapted 
from [55]. 
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antiparallel to the field, while the two others become more aligned with it. Increasing the 
field causes a first order phase transition (at /_i B = 3 T) into a commensurate phase; see 
Fig. 5.15. In this phase, the two spins are parallel to the field, and the third is antiparallel to 
it. Since this is a first order transition, the angle between the two nearly parallel spins does 
not continuously decrease to zero, but rather abruptly when the critical field is reached. The 
higher field transition (at B = 9 T), however, is a second order; see Fig. 5.16. Above the 
transition (the second incommensurate phase), the minimum energy configuration consists 
of spins slightly and continuously rotating away from the direction of the field; the spin 
nearly antiparallel more so than the two that are nearly parallel. Theory predicts another 
transition at even higher fields, into a phase in which all three spins become aligned with 
the magnetic field. Such a transition should occur at fields much higher than experimentally 
available in a neutron scattering experiment; accordingly, we haven't observed it. 

We observed hysteresis at lowest temperatures; see Fig. 5.14. At T = 1.7 K and 
B = 6 T, the zero-field cooled phase is incommensurate (A = 0.453), while the field 
cooled system ends up in a commensurate phase (A = 1/3). We mapped the approximate 
H-T boundary of the region where hysteresis is significant. The analysis suggests that the 
two ground states (commensurate and incommensurate) are close in energy, within about 
E 0.1 meV. 



Spin waves 

Using inelastic neutron scattering experiments, we observed well-defined spin waves 
in the ordered phase. These measurements were performed using SPINS spectrometer, 
in a focusing analyzer configuration. We observed spin waves both in-plane, and along 
the c-axis, and mapped the dispersion relations. The initial analysis indicated that there is 
a large in-plane anisotropy. Accordingly, the original Heisenberg Hamiltonian had to be 
augmented: 

H = J ^ Sj • Si in plane Heisenberg 

m 

+ D y^fff) 2 in plane anisotropy 

i 

+ J 2 Si ■ S k vertical interplane 

{ik} 

+ J 3 ^ Si • Si + J 4 ^ Si • Sj<; diagonal interplane (5.15) 

{ii} {ii'} 

where {ij} denotes all in-plane nearest neighbor pairs; {ik} denotes all interplane nearest 
neighbor pairs; {il} and {il 1 } denotes all interplane next-nearest (diagonal) neighbor pairs. 
Likewise, J, J 2 , J 3 , and J 4 denote exchange interactions between these pairs, respectively 
(see Fig. 5.11. The strength of a single site anisotropy is measured by D. Some information 
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Figure 5.14: Hysteresis in the low temperature magnetic order in RbFe(Mo0 4 ) 2 . Figure 
shows magnetic Bragg peaks along the Q = (1/3 1/3 I) line, at T = 1.7 K and fj, H = 6 T. 
Circles are for the zero field cooled sample; the magnetic structure is the low temperature 
A = 0.453 incommensurate structure. Triangles are for /_i H = 10 T field cooled sample; 
magnetic structure is the A = 1/3 commensurate structure. Lines are Gaussian fits to the 
data. Adapted from [55]. 
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Figure 5.15: The first order phase transition in RbFe(Mo0 4 ) 2 , at T = 1.7 K. Open circles 
show the intensity ofQ = (l/3 1/3 0.453) magnetic Bragg peak, and bullets the intensity 
of Q = (1/3 1/3 1/3) peak. The transition, occurring around y H = 4.2 T is an 
incommensurate to commensurate phase transition. Adapted from [55]. 
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Figure 5.16: The second order phase transition in RbFe(Mo0 4 ) 2 , at T = 1.7 K. Open 
circles show the intensity of the sliding Q = (1/3 1/3 ~ 0.4) magnetic Bragg peak, 
and bullets the intensity of Q = (1/3 1/3 1/3) peak. The transition, occurring around 
fx H = 9 T is a commensurate to incommensurate phase transition. Adapted from [55]. 
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Figure 5.17: Spin waves along the c-axis, seen in a constant-Q scan in RbFe(Mo0 4 ) 2 . (a) 
T = 1.8 K, Q = (1/3 1/3 0.2). (6) T = 1.8 K, Q = (1/3 1/3 0). Measured using 
SPINS spectrometer at NCNR; £ f = 3.7 meV. Adapted from [55]. 



about these constants can be gained by studying static order; the analysis and results of this 
work are reported in Sec. B.5 (p. 158). The study of dynamics is much more complicated, 
and also a part of continuing work. Current results reveal that J = 0.096 meV, and D = 
0.024 meV. Large in-plane spin anisotropy indicates that our system is better described as 
an X-Y, rather than a Heisenberg, antiferromagnet on a triangular lattice. Incorporating 
interplane interactions is complicated, as there are three different exchange constants, all 
of which are relevant for interplane ordering. 

Kosterlitz-Thouless (KT) transition 

One of the most interesting transitions in two-dimensional magnetic systems, the Kosterlitz- 
Thouless transition, was theoretically described more than three decades ago [57, 58]. 
Good experimental systems for its study, however, remain rare [53, 59]. We attempted 
a study of KT transition in RbFe(Mo0 4 ) 2 . Let us briefly describe the KT transition; we fol- 
low the approach given in Ref. [61]. Let us define an order parameter for an X-Y magnet: 



M(r) = \m(r) | Re [(x + iy)e^ (r) ] . 



(5.16) 
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Figure 5.18: Dispersion relations of spin waves along the c-axis, in RbFe(Mo0 4 ) 2 . Adapted 
from [55]. 



For our triangular magnet, the staggered order parameter may be defined as: 

0(r) = (j>(r) ± 120° x (h + k), (5.17) 

where r = ha + kh. Thus, for a perfect 120° spin configuration, staggered order param- 
eter would be zero: 12 o° (r) = 0. Excitations referred to as spin vortices, or topological 
charges, Q, are then defined as: 

Q can take only integer values. Pairs of vortices having no net topological charge, Q net = 0, 
appear at all temperatures. Isolated topological charges, however, appear only above some 
critical temperature, called the Kosterlitz-Thouless temperature, T K t, approximately given 

^ = 2,exp( F ^). (5.19) 

Likewise, the transition is referred to as Kosterlitz-Thouless transition. The instantaneous 
spin-spin correlation length function defined via: 

(S S r ) oc — — , for r > f = 1/2™, (5.20) 
is a correlation length) has a temperature dependence: 

f oc exp(6/v / *), (5.21) 
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Figure 5. 19: Spin waves in the basal plane, seen in a constant-Q scan in RbFe(Mo04) 2 . (a) 
T = 1.8 K, Q = (0.22 0.22 0.453). (b) T = 1.8 K, Q = (0.16 0.16 0.453). Measured 
using SPINS spectrometer at NCNR; E t = 2.5 meV, with focusing analyzer. Adapted from 
[55]. 
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Figure 5.20: Dispersion relations of spin waves in the basal plane, in RbFe(Mo0 4 ) 2 . 
Adapted from [55]. 



where t = (T — T KT )/T KT , and b ~ 1.5. (Estimates for rj, which is temperature dependent, 
are: 0.25 < rj < 0.35.) At and below the transition, £ diverges (k — > 0), and we have a 
power law decay: 

{S S r ) a r-^ T) , for T = T KT - (5.22) 
We tried to detect the tail of the KT transition in our experiment. Our estimates, using: 

T KT ~ 0.28 x JS(S + 1), (5.23) 

and J = 0.096 meV, give T K t ~ 2.7 K. At that temperature, the system is magnetically 
ordered, and so the KT transition cannot occur. Above the ordering temperature, however, 
the spin-spin correlation length may follow Eq. 5.21, signaling the proximity of the KT 
transition. To experimentally measure £, we devised a particular type of scan, which we 
dubbed arc scan, designed so it maximizes the neutron intensity and resolution. While 
we were able to observe the temperature dependence of the inverse correlation length, k, 
our data had abhorrent spurious peaks which were difficult to remove, and prevented us 
from extracting reliable results. Though unsuccessful at present, we are hopeful that future 
experiments could yield exciting and important discoveries. 

5.7 Conclusions 

When we picked RbFe(Mo0 4 )2 for study, the hope was that the interlayer exchange 
coupling will be much weaker than the intralayer coupling. This hope was sustained by 
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the lack of formal chemical bonds between layers (chemical bonds in this material con- 
nect atoms only within one layer); thus we expected the interlayer superexchange paths 
to be negligible. That would make the system a good realization of a 2D magnet. On 
the other hand, the relative proximity of interlayer ions, illustrated by the small inter- to 
intralayer distance, c/a = 1.310, suggested significant interaction in the third dimension. 
The Mermin-Wagner theorem stipulates that no long-range order can break a continuous 
symmetry in two dimensions [62]. We found, however, magnetic order below T = 3.5 K; 
thus the interplane interactions are apparently significant. 

We calculated the ordered moment: (S) Rj 0.75 S, and found that a large portion of the 
available spin remains in the fluctuating state, even at the lowest temperatures. We mapped 
a rich H-T phase diagram, with (i) commensurate and incommensurate phases, (ii) first 
and second order phase transitions, and (in) low temperature hysteresis. We observed 
well-defined spin waves, and extracted significant single site anisotropy. Future research 
should include experimentally determining the high field magnetic structure, completing 
the spin wave analysis, and a more complete experimental study of the Kosterlitz-Thouless 
transition. 
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Chapter 6 

Kagome Sandwich Antiferromagnet 
QS-ferrite 

6.1 Introduction 

Classically, phase transitions are driven only by thermal fluctuations, and at T = 
such systems usually settle into a long range ordered ground state [65]. In a quantum sys- 
tem, however, fluctuations arising from the Heisenberg uncertainty principle persist even at 
T = 0, and can cause a dynamic state at T = 0. Such systems can display a quantum phase 
transition as some of the variables characterizing the system (spin, dilution concentration), 
or its environment (pressure, magnetic field, etc., but not the temperature) are varied. The 
effect of these variables is commonly parameterized by a generic coupling constant g, and 
the phase transition is defined as a point, called a quantum critical point (QCP), g = g c , at 
which the ground state energy, as a function of g, becomes nonanalytic. The T- g phase dia- 
gram around a QCP can then be constructed as follows. If, moving away from the quantum 
critical point along the T = line, the fluctuations abate, the system goes into a ground 
state with a long range order (a Neel state for magnetic systems), and by definition g < g c . 
At finite temperatures, this part of the phase diagram is called renormalized classical (RC), 
as the system is then described by the classical renormalization theory. In contrast, when 
g > g c , quantum fluctuations suppress the ordering, and the correlation length of the order 
parameter becomes finite and independent of temperature as T — > 0; we call this regime 
quantum ordered (QO). [This regime is also commonly known as quantum disordered, thus 
referring to the lack of the static long range order. The term is, however, somewhat mis- 
leading: ground state is best described as a coherent quantum paramagnet, with a long 
coherence length.] 

Since a quantum phase transition occurs only at T = 0, it cannot be directly observed in 
an experiment, which is necessarily carried out at some finite temperature. The transition, 
however, affects the properties of the system in the T > region around g = g c , which is 
called quantum critical (QC), and can, therefore, be observed indirectly. At finite temper- 
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atures, it is the free, rather than the internal, energy that determines the state of the system, 
and it can either (i) be analytic in both T > and g in the quantum critical region, in which 
case the system exhibits only a quantum phase transition (i.e. non-analyticity is confined to 
only one point, the QCP), or (ii) it can have a line of non-analytic points, corresponding to 
a finite temperature second order transitions, terminating at the QCP. The quantum critical 
region discussed in this chapter is of the second type. The transition, at sufficiently high 
temperatures, is primarily driven by thermal fluctuations, and can therefore be described by 
the theory of classical phase transitions. In the vicinity of the quantum critical point, how- 
ever, an interplay of quantum and thermal effects can yield novel properties which cannot 
be accounted for by the classical theory. 

At the QC transition, dynamics of a <i-dimensional system is determined by classical 
critical dynamics in d + 1 dimensions [66, 67]. The extra dimension is imaginary time, 
and at non-zero temperatures it is reduced to a finite thickness, ~ ch/k B T. This leads 
to a universal quantum critical scaling of the spin dynamics, huj/k^T, which has been 
observed experimentally in some materials [68, 69]. The excitation spectrum in the QC 
region is gapless, and as the system moves away from the QCP, an energy gap develops, 
A oc \g — g c \ Zv . In magnetic systems, this leads to an ordered state in the RC regime, which 
is easily detected using elastic neutron scattering, by the presence of the magnetic Bragg 
peaks. In the QO region, where the system remains in a paramagnetic state, it establishes an 
energy scale, A, independent of temperature, that can be measured with inelastic neutron 
scattering. Determining the scaling dependence can then be used to experimentally place a 
magnetic system in the QC phase diagram, and to confirm the existence of QCP. (Ref. [70] 
is an interesting and recent paper in which a crossover from QC to QO regime has been 
studied by neutron scattering.) 

Although, in principle, quantum fluctuations dominate in all physical systems at suf- 
ficiently low temperatures, magnetic systems offer ideal setting for studying such effects. 
With magnetic systems, a variety of qualitatively different models can be experimentally 
tested, and the ideas and results often extended to other, non-magnetic, problems. This 
unique role of magnetism arises from large variability of model parameters. For instance, 
the dimensionality of both the system and the degree of freedom can be independently 
selected. The cooperative behavior can be studied via magnetic and thermodynamic mea- 
surements, and microscopic insight can be gained by techniques such as neutron scattering 
and /j,SR. Disorder can be introduced nearly continuously, and its effects studied by di- 
luting magnetic ions with non-magnetic, but chemically similar, ions. Remarkably, the 
effect of quantum fluctuations can be isolated by changing the effective atomic spin, lower 
values enhancing quantum effects. It is unfortunate that spin, being a discrete variable, 
cannot be arbitrarily tuned, and often the critical spin value is not experimentally attain- 
able (i.e. S c < 1/2). Luckily, another feature, frustration, characterizing some magnetic 
systems, dramatically enhances fluctuations at experimentally accessible spins and temper- 
atures. 

A system is frustrated when all the local interactions cannot be simultaneously satisfied. 
In magnetism, it arises from the incompatibility of local interactions, usually antiferromag- 
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netic, with the symmetry imposed by the crystal structure. Frustration is commonly found 
in disordered systems, such as spin glasses; however, it is interesting to separate the ef- 
fects of quantum fluctuations from those of disorder. Luckily, there is a class of materials, 
known as geometrically frustrated magnets (GFM), on which frustration exists for a per- 
fectly ordered lattice. (Reviews of GFM are given in Ref. [2, 1,3].) The simplest and most 
common example of a GFM is a triangle the vertices of which are occupied by spins inter- 
acting antiferromagnetically, where it is impossible to simultaneously satisfy all the local 
interactions. Replicating this basic unit, which in three dimensions becomes a tetrahedron, 
large scale frustration can be achieved, without introducing any disorder. 

Frustrated magnets constitute a distinct class of magnetic materials, separate from both 
unfrustrated and disordered magnets. Like conventional antiferromagnets, frustrated mag- 
nets have a negative Curie- Weiss temperature: 

Gcw = \ S(S + 1) J/k B < 0. (6.1) 

Unlike the former, however, they do not order at an expected mean- field ordering tem- 
perature, |Ocw|> but remain in a cooperative paramagnetic state down to a much lower 
transition temperature, T c . The characteristic signature of strongly frustrated magnets is, 
therefore, that their inverse susceptibility follows the Curie- Weiss law, 

down to temperatures well below the expected mean-field transition temperature, where a 
non-Neel state develops. (The ratio of these two temperatures, / = |9 CW |/T C , has been 
proposed as a quantitative measure of geometric frustration, and magnets with / > 10 are 
said to be strongly frustrated.) The purpose of studying GFMs, in the context of quantum 
phase transitions, is then to reduce the transition temperature, while keeping the interactions 
strong, thus bringing the system into a quantum critical region. 

Although frustration certainly enhances fluctuations, the simplest frustrated system, the 
triangular lattice, nonetheless orders at sufficiently low temperatures, even for the lowest 
spin value, S = 1/2 [71]. The triangular lattice can be seen as an array of edge sharing 
triangles, which are strongly coupled (coordination number 6), and therefore spins on this 
lattice tend to develop strong correlations. The QCP in this system is at S c < 1/2, thus 
unattainable in any physical system. To suppress ordering, the basic frustrated units (trian- 
gles) need to be weakly connected with each other. One such frustrated system, displaying 
weak connectivity, is a kagome lattice, constructed by dividing a triangular lattice into four 
sublattices, and removing the atoms on one of the sublattices. The resulting array consists 
of corner sharing triangles (coordination number 4), the vertices of which are occupied 
by magnetic ions. Theoretical results for S=l/2 spins on a kagome lattice indicate that 
there is no long range order at any temperature [72], while in the classical limit, S — > oc, 
the system displays order-by-disorder as T/|0 C w — > [73, 74]. It is intuitively clear 
that fluctuations create disorder; thus the concept of order being promoted by fluctuations 
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seems paradoxical. However, when a system has a highly degenerate ground state (as is the 
case in kagome lattice [75]), without energy barriers separating different states, quantum 
or thermal fluctuations can lift the degeneracy and cause the system to order into a state 
with the highest entropy; this is referred to as order-by-disorder [76]. Further increase in 
fluctuations (as a result of rising temperature, for instance), eventually destroys the order it 
at first stabilized. 

Thus, if we take the spin quantum number as a coupling parameter, g ~ 1/5, we 
conclude that the 5 = 1/2 and 5 = oo lie on the opposite sides of a QCP, which should 
to be approached by some intermediate spin value. Furthermore, experiments on kagome 
based compounds suggest that for 5 = 5/2, kagome systems do order at low temperatures 
[77] , thus lowering an upper bound for the spin value corresponding to the QCP. This makes 
5 = 3/2 kagome system an attractive candidate for QCP studies. 

Physical systems based on a true two-dimensional kagome lattice are, unfortunately, 
rare [78]. Instead, quasi-two-dimensional compounds, in which the majority of spins reside 
on a kagome net, have been synthesized and studied. One of the best studied such com- 
pounds, SrCr 9p Ga 1 2_9 J9 0i9 (SCGO(p)) [79], consists of a kagome sandwich: two (slightly 
distorted) kagome layers separated by a triangular lattice. (Alternatively, a kagome sand- 
wich can be seen as a three-layer (111) pyrochlore slab, or a hexagonal lattice with a unit 
cell consisting of two corner sharing tetrahedra, thus containing seven sites.) Magnetic 
sites are occupied by Cr 3+ ions, with 5 = 3/2 spins. Furthermore, pairs of strongly cou- 
pled (J = 20 meV) Cr 3+ ions are located between sandwiches. Since they are only weakly 
coupled with kagome layers (J k, 0.2 meV), it was assumed that sandwiches are approxi- 
mately isolated from each others [78]. It was, however, uncertain how good approximation 
this is at low temperatures. Next, since all the samples have p < 0.92, magnetic sites 
are diluted and disorder introduced. Since the percolation threshold for kagome lattice 
is 65.27%, this dilution does not limit the correlation length [80]. However, it remained 
unclear how important disorder is for the glassy behavior that was observed in SCGO. 

Initially it was thought that the kagome lattice might be sufficient to capture the es- 
sential properties of kagome sandwich. Recent work however indicates that the two have 
some qualitatively different properties (for instance, a finite temperature chiral Kosterlitz- 
Thouless transition is predicted only for a kagome sandwich, and is absent in the single 
layer system [81]). It appears that the spin-glass-like transition observed in SCGO is repro- 
duced only on a pyrochlore slab model, and not on a true 2D kagome lattice. In the absence 
of pure 2D kagome systems, a kagome sandwich without disorder would be desirable, to 
approach the intrinsic effects of frustration and low connectivity. 

Another compound synthesized recently, Ba 2 Sn2Ga 3 ZnCr7022 (BSZGCO), with the 
"(25 ferrite" structure, resolves some of the problems just posed. Structurally similar to 
SCGO, its magnetic lattice is free from disorders, and sandwiches are well separated (> 9 
A), without any magnetic ions in between. Thus, it represents a well isolated S=3/2 kagome 
sandwich, and is the cleanest material known for its study. 

Recent work by Hagemann et al. [9] revealed strong antiferromagnetic interaction, with 
the Curie-Weiss temperature of 6 C w = —312 K. As in other frustrated magnets, it's sus- 
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(a) 



(b) 




Figure 6.1: Kagome lattice. Ground state is highly degenerate, and there are zero energy 
excitations. Examples: (a) Q = y/3 x y/Z ground state spin configuration. Any state in 
which all the spins below dotted line rotate around the vertical (dashed) line (blue spins), 
is still a ground state, (b) Q = ground state spin configuration. Rotation of spins be- 
longing to the inner hexagon (blue spins) around the vertical (dashed) line is a zero energy 
excitation. 
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Figure 6.2: Magnetic lattice of QS ferrite (E^Sr^ZnGasG^C^)- Shown are Cr 3+ ions that 
form kagome layers, with selected oxygen atoms. Sizes of atoms correspond to their ionic 
radii. Sticks represent superexchange paths; different colors indicate different exchange 
strengths. 



ceptibility follows Curie-Weiss law (6.2) to well below 6cw. and has a weak maximum at 
T sg « 1.5 K, reminiscent of a spin glass transition. 

In this chapter we describe neutron scattering experiments on BSZGCO. Using elastic 
scattering, we probed the frozen spin moment at different timescales; we found the results 
to be timescale dependent, as is the case in spin glasses. With inelastic measurements, we 
studied the excitation spectrum, and found its wavevector dependence to be independent of 
temperature, and vanishing as Q — > 0; this supports the idea that the dynamics is governed 
by an interaction between satisfied simplexes. We observed the quantum critical hu/k-oT 
scaling, and broadening of the distribution of the relaxation rates at low temperatures. 
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Table 6.1: Structure of QS ferrite (Ba2Sn 2 Ga3ZnCr 7 022) at room temperature. Space 
group: P3ml; a = 5.85217(7) A, c = 14.2390(2) A. Site Zn/Ga occupied by 0.5 Zn 
and 0.5 Ga. R p = 4.40%, R wp = 5.40%, x 2 = 1-578. Adapted from [9]. 

6.2 Sample and experimental details 

An 1 1.2 g powder sample of Ba2Sn 2 Ga3ZnCr 7 022 (QS ferrite) was prepared by R. J. Cava 
et ah, using methods reported in Ref. [9]. Figure 6.2 shows the structure of QS ferrite, 

o 

which consists of well separated (9.4 A), quasi-two-dimensional, pyrochlore slabs. Each 
slab contains triangular lattice sandwiched between two kagome layers. Ions in the kagome 
layers are surrounded by four oxygen atoms in a tetrahedral, while those in triangular layer 
are surrounded by six oxygen atoms in an octahedral coordination (as shown for selected 
Cr 3+ ion). The distances between the magnetic ions are labeled on the figure, and we note 
that the kagome lattice is slightly distorted. Different superexchange paths resulting from 
this are marked in the figure by different stick colors. It should be noted that the compound 
has chemical disorder on the 2d site shared by Zn and Ga atoms; however, as both atoms 
are nonmagnetic, the magnetic lattice remains ordered. 

The measurements presented here were carried out on the Disc Chopper Spectrometer 
(DCS), and the High Flux Backscattering Spectrometer (HFBS) at the National Institute 
of Standards and Technology in Gaithersburg, Maryland. DCS is a time- of- flight neutron 
spectrometer with 913 detectors, enabling simultaneous collection of data over a range of 
energies and wavelengths. To facilitate the measurement of the relaxation rate over a wider 
temperature range, the DCS measurements were taken in two configurations, one with the 
incident neutrons wavelengths fixed at A = 2.3 A, probing energies to about 10 meV, with 
the resolution of 1.1 meV (elastic channel). The other setting used A = 4.8 A neutrons, 
for energies up to 2 meV with 0.12 meV resolution. HFBS is a 16 detector spectrometer 
operating in a backscattering configuration, thus enabling energy resolution of only 0.80 
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6.3 Experimental results 

Color Figure 6.3 shows normalized, background subtracted, neutron intensity obtained 
from DCS at 4 selected temperatures. Before plotting, the data has been convoluted (for 
this figure only) with a simple box function, in order to diffuse the artificial sharp features 
arising from statistics. The most prominent feature present on all the images is the bright 
horizontal line of intensity arising from the incoherent elastic scattering. It is easily iden- 
tified, as it is nearly independent of the wave-vector, centered at hu> = 0, and has a width 
determined by the instrument energy resolution. Little can be learned of the dynamics of 
the system from it; however, it is indispensable for normalization of intensity into abso- 
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lute units. Magnetic Bragg peaks have not been observed even at the lowest temperature, 
T = 1.7 K <C ©cw> indicating that QS-ferrite does not develop static long range spin order. 
Therefore, we conclude that quantum effects dominate at low temperatures in BSZGCO, 
and if a QCP exists, g > g c . 

The area of increased intensity between 1 and 2 A -1 is associated with the excitations 
in the system. As the intense patch extends over a range of energies, the excitations must 
involve a macroscopic number of degrees of freedom. (The lack of intensity at low temper- 
atures at negative energies — the magnon absorption processes — is a simple consequence 
of the detailed balance principle.) These measurements should be compared with those 
performed on SCGO (Fig. 2 in Ref. [82]). The most obvious qualitative difference is the 
lack of the constant energy, hcu = 18.6 meV, ridge of intensity. Since that ridge arises 
from the isolated spin pairs located between kagome sandwiches, which are present only 
in SCGO, their absence in BSZGCO explains the difference in the spectrum. 

The nature of the spatial correlations can be inferred from the Q-dependence of the 
inelastic scattering. It displays a broad peak centered around Q = 1.5 A -1 , indicating 
short-range antiferromagnetic correlations between spins. If the correlations exist only 
within the kagome planes, the wavevector should lie in the basal plane and correspond 
to a high symmetry point in the hexagonal reciprocal lattice with lattice parameter a* = 
47r/-y/(3)a = 1.24A -1 . The two highly symmetric coplanar spin configurations on kagome 
lattice, q = and y/Z x would produce peaks at (1 0) and (| |), corresponding to 
Q = 1.24 A -1 and Q = 1.43 A -1 , respectively. The data seems to favor the latter, although 
single crystal diffraction is needed to establish this beyond doubt. In SCGO, for which the 
spectrum was qualitatively the same as for BSZGCO, and the inelastic intensity centered 
around Q = 1.4 A -1 , single crystals were available, and the experiments have shown that 
indeed \/3 x \/3 was the short range configuration [83, 84]. 

Vanishing forward scattering, I(Q — > 0) — >■ 0, suggests the absence of a net spin on 
the basic frustrated unit, a simplex. This suggests that the interactions between simplexes 
involve non-trivial degrees of freedom, which keep the total simplex spin zero. Such com- 
posite degrees of freedom have recently been proposed for a spinel lattice, a parent of 
the QS ferrite lattice, in which clusters of strongly interacting spins, rather than the spins 
themselves, govern the dynamics at low temperatures. (Such clusters, if they exist in a 
pyrochlore slab, may be different from those in the spinel lattice.) The clusters are rather 
stable at temperatures T <C |©cw| = 312 K, and interact weakly amongst themselves, the 
strength of the interaction being of the order of SG transition. It should be noted that there 
is no perceptible change in the Q dependence of the excitation spectrum over the entire 
temperature range: It remains broad and centered around Q = 1.5 A -1 . This may be due 
to the cluster form factor. The size of the cluster, inversely proportional to the width of the 
spectrum distribution, seems to be slightly above the nearest neighbor chromium distance. 
This is consistent with the classical theory, in which the cluster would be a tetrahedron, 
on which the energy is minimized when the four spins add to zero. Since the number of 
degrees of freedom for four spins is 4 x 2 = 8, and number of constraints imposed by net 
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Figure 6.4: Temperature dependence of the elastic magnetic scattering (frozen spin mo- 
ment). Circles: A = 4.8A data (DCS). Diamonds: HFBS data. Lines are guides to the 
eye. 



spin zero is 3, this yields a macroscopically degenerate ground state characterized by two 
parameters per tetrahedron (excluding 3 degrees of freedom corresponding to the rotation 
of all spins) [75]. Since in a powder the analysis is complicated by the fact that the intensity 
is spherically averaged, it is not possible to uniquely determine the magnetic structure of 
the cluster. 

Hagemann et al. discovered that BSZGCO has a very broad peak in C(T) /T at T m 
3.7 K, reminiscent of a spin glass transition, which they find around T sg »s 1.5 K. Using 
HFBS in fixed window mode, we measured the nominally elastic neutron scattering at low 
temperatures. The data in Fig. 6.4 shows the development of magnetic correlations on 
time scales r > h/AE. With the energy resolution of only 0.80 /ieV, this corresponds 
to t > 0.8 ns. As shown by solid line (diamonds) on Fig. 6.4, at this timescale spin 
correlations develop below T 7 K, and frozen moment reaches about 45% polarization 
at the lowest temperature accessed, T = 1.7 K, indicating that a large fraction of spin still 
remains in a fluctuating, unfrozen state. 

Spin dynamics was probed on the same instrument via inelastic neutron scattering, and 
Fig. 6.5 shows normalized inelastic intensity at a few selected temperatures. The excitation 
spectrum can be described by a Lorentzian of width T, corrected for the detailed balance 
principle, indicating that the autocorrelation function decays exponentially (Sj(0) • Sj(i)) = 
e~ r< . The lines through the data show the inelastic function (Lorentzian) added to the 
Gaussian describing the elastic peak at Hco = which arises from the elastic incoherent 
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Figure 6.5: Normalized hu dependence of the inelastic neutron scattering, obtained using 
HFBS at NIST. 
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scattering. (There are no Bragg peaks in this 0.6 A -1 < Q < 1.8 A -1 window.) The 
inset to Fig. 6.10 shows the temperature dependence of the relaxation rate. At the lowest 
temperatures, the corresponding relaxation time is 1/T = 0.2 ns, and rapidly decreases as 
the temperature is raised. 



6.4 Analysis and discussion 

The amount and quality of data collected on the DCS instrument warranted a more 
elaborate analysis. We have already observed in Fig. 6.3 that the wavevector dependence 
of the excitation spectrum does not change with temperature. Led by this observation, we 
proposed separating the excitation spectrum (which is related to the imaginary part of the 
spin susceptibility via the fluctuation dissipation theorem, nS(Q, u) = (n T (uj) + 

1)), into a product of the energy only, and momentum only, dependent functions: 

S(Q, uj) = x(Q) r F^~r~ 2 <M") + 1>, (6-3) 

where the Lorentzian energy dependence is apparent from the equation, while the momen- 
tum dependent function, was presumed unknown. Both functions are temperature 
dependent, as is the boson occupation number, (n T (uj) + l) = (1 — e - fi W*BT-)-i DCS 
inelastic neutron intensity was then fit to: 



d 2 a k 



f Ji 



dfl du kj 



2 



J R 5e (oj' - co) S(Q, oo') du'. (6.4) 



Here, R S e(oo') is the energy resolution, well approximated by a Gaussian of the width SE, 
which varies with the energy transfer as SE m 0.019 x (E, - E) 3 ^ 2 meV -1 / 2 . The actual 
DCS energy resolution function, however, deviates considerably from a Gaussian, having 
extensive tails with intensity of about 1% of the maximum. Though this is insignificant for 
the inelastic data, the elastic peak profile had to constructed before the fit, as the intensity 
of the tail is comparable to the inelastic signal. The fit to equation (6.4) was performed 
on the set of data points distributed symmetrically about the huj = in the Hu-Q plane, 
simplifying verification of the detailed balance principle, as well as the background sub- 
traction. 

For the purpose of the fit, x(Q) was a t first explicitly defined on a discrete set of points 
in the relevant range of Q values spaced by 0.2 A -1 apart, with interpolation in between. It 
was found that the Q dependence of the inelastic scattering, x(Q)> is well described by: 

X (Q) = C(T) J2A n (l- ) , A ee 1; (6.5) 
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Figure 6.6: Normalized Q-dependence of the inelastic neutron scattering (A = 2.3A, DCS). 



where the summation is carried out over the nearest, and the next-nearest, neighbor dis- 
tances (d 2.92 A, and di = 5.07 A, respectively). The neutron intensity of the lattice 
decreases exponentially at high Q due to the Cr 3+ magnetic form factor, F(Q), and the 
extraction of x(Q) g ets quite uncertain at momentum transfers above 2.5 A -1 ; hence the 
summation was limited only to these two terms. We found x{Q)jC to be temperature in- 
dependent (Ai = —0.25(5)), as the Figure 6.6 demonstrates for A = 2.3 A data at three 
different temperatures (the result was in agreement with the A = 4.8 A data, which, how- 
ever, spanned only lower temperatures). After ascertaining that the quality of the fit was 
unchanged, we dropped the "tabular" definition of x(Q) m favor of that given by equation 
(6.5). 

The persistence of a particular wavevector dependence, given by Eq. (6.5), over so 
wide temperature range supports the notion of a simplex scattering form factor introduced 
earlier. It should be noted that the actual Q dependence of the excitation spectrum is a 
product of x(Q), and the Cr 3+ magnetic from factor, F(Q), which causes the measured 
intensity to essentially vanish above Q > 2.5 A -1 . 

Having established the Q dependence of the excitation, we proceeded to explore its 
energy dependence, starting with the A = 2.3 A data which covers the widest temperature 
range. Figure 6.7 shows the Hco dependence of the data used in the fit at a few selected 
temperatures. We should remark, however, that the fit line in the Figure does not simply 
depict the energy dependence as expressed in equation (6.3). Instead, the intensities of all 
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Figure 6.7: Normalized hco dependence of the inelastic neutron scattering, A = 2.3A data 
(DCS). The uncertainties in the data are comparable to the symbol size. 
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Figure 6.8: Relaxation rate obtained from Lorentzian fit. Solid circles: QS ferrite; fit to 
A = 2.3A data (DCS). Squares: SCGO(p) data from [83]. The dashed line represents 
fir = k B T. The solid line is a fit to T oc T 67 , and the dotted line a fit to T oc T o n . See 
solid circles in Fig. 6.10 for a closer look at temperatures below 30 K. 
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the data points used in this fit have been "rebinned" and projected onto the Hlu axis: 

1 N 

i 

the same being done with the fit function, thus enabling a visual evaluation of the fit. 

The values of V so extracted are shown in Figure 6.8. A fit to T oc T a yields a = 
0.67(2). The dashed line represents HT = k B T, and we see that the relaxation rate fol- 
lows it at low temperatures. There is no intrinsic energy scale; instead, uj/T scaling is 
observed, demonstrating that the system is close to a QCP. Also shown are the data for 
SCGO(p=0.92(5)) [83], forwhich a = 0.71(2). Since for both systems the magnetic lattice 
is a kagome sandwich, it is interesting that so similar exponents are found. The relaxation 
rates are, however, much larger in SCGO, a result that cannot be fully accounted for by the 
larger coupling constant J. The origin of this discrepancy is uncertain, and we speculate 
that it may be caused by the disorder, which is present only in SCGO. 

With their large energy resolution (8E E=0 »s 1.1 meV), the A = 2.3 A data were 
inadequate for reliably extracting T at the lowest temperatures; accordingly we used the 

o 

more acute A = 4.8 A data (SE E =o ~ 0.12 meV). Figure 6.9 shows the data, and the fit 
(the dashed line), at three representative temperatures. At higher temperatures the energy 
dependence is well described by a Lorentzian spectral function, while at lower temperatures 
spectral weight shifts toward lower energies, and the anomalously flat distribution deviates 
from Lorentzian. Noticing the existence of a small, but distinctly visible, residual intensity 
below 0.5 meV for all the temperatures under about 11 K, we expanded the excitation 
function (6.3) by adopting two relaxation rates, I\ and T 2 . We fit the data with two different 
functions. The first is simply two distinct Lorentzians: 



S.M = + (6.7) 



FiU) r 2 o> 

The second is a uniform ("box") distribution of Lorentzians extending from Ti to T 2 



1 [ T2 co dT 



tan" 1 (r 2 /w) - tan" 1 (T 1 /uj) 



(6.8) 



With the addition of this one parameter, we could now account for the residual in the data; 
however, we could not discern between the two functions, as they both fit the data equally 
well at all temperatures. Since above about 10 K a single Lorentzian describes the data well, 
T 2 in equation (6.7) becomes redundant, and, during the fit, causes numerical instability; 
whereas in equation (6.8) both Yi and T 2 remain well defined, and converge toward a 
common value given by a single Lorentzian fit. For this reason we adopted equation (6.8) 
to fit the excitation spectrum. 
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Figure 6.9: Normalized fuo dependence of the inelastic neutron scattering, A = 4.8A data 
(DCS). Excitation spectrum is given by: Dashed line, equation (6.3); Solid line, equation 
(6.8). The uncertainties in the data are comparable to the symbol size. 
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Figure 6.10: Relaxation rate. Solid circles: Fit of A = 2.3A data (DCS) to a single 

o 

Lorentzian. Open circles: Fit of A = 4.8A data (DCS) to a single Lorentzian. Trian- 
gles, regular and inverted: Fit of A = 4.8A data (DCS) to a uniform "box" Lorentzian 
distribution between r\, and T 2 , respectively. Diamonds (inset): Fit of HFBS data to a 
single Lorentzian. Dashed line: HT = k B T. 



121 



The relaxation rates obtained for temperatures below 30 K are shown in Figure 6.10: 
solid circles for A = 2.3 A, and open circles A = 4.8 A, data, both fit to Eq. (6.3) (single 
Lorentzian). At higher temperatures, the two results agree, and it seems that hT ~ k^,T 
(dashed line). Triangles represent Yi and T 2 of the "box" distribution, as defined by 
Eq. (6.8). The quality of the single Lorentzian fit (measured statistically by x 2 > as we U 
as visually, Fig. 6.9) decreases below about 11 K, the temperature at which the deviation 
between it and the "box" distribution becomes significant. 

The data, therefore, indicates that as the system is cooled below ^11K, the relaxation 
distribution rate broadens, and remains such to the lowest temperatures accessed. Further- 
more, below ^2.5 K, as the temperature approaches the glass transition, the lower cutoff 
for the relaxation rate distribution, Ti, drops to zero. Recent work uncovered a similar 
relaxation rate broadening in ZnCr 2 4 , an S=3/2 spinel lattice which is a parent compound 
of QS ferrite [7]. The authors of that paper associated the vanishing lower cutoff of the 
relaxation rate with the spin glass transition. 

We already mentioned the measurement of the diffuse elastic intensity, with the HFBS 
instrument, signaling the development of short range magnetic correlations at timescales 
r > 0.8 ns for temperatures T < 7 K. The A = 4.8 A -1 DCS data enabled us to measure 
correlations at shorter timescales, r > 5 ps (it is unfortunate that for A = 2.3 A -1 data, 
the resolution was too coarse to separate elastic magnetic scattering from nuclear Bragg 
peaks). Figure 6.4 shows that at this timescale correlations develop at higher temperatures, 
T 15 K, and the frozen moment reaches higher values, about 70%, as T — > 0. (Frozen 
moment should, of course, be higher at shorter timescales, since fluctuations vanish for 
r — > 0.) This is a common property for spin glasses, that the response function depends 
on the timescale at which the system is probed. This sets spin glasses apart from con- 
ventional magnets, in which no significant timescale dependence is observed above some 
characteristic microscopic timescale of the system [85, 86]. 

6.5 Conclusions 

We shall now examine the effects of frustration and disorder on the anomalous, spin- 
glass-like, properties observed in SCGO and BSZGCO. The properties of spin glasses were 
traditionally attributed to both frustration and site disorder, and until recently, all real ma- 
terials displaying such behavior have had some structural disorder [7]. Geometrically frus- 
trated magnets, in which frustration is achieved without disorder, enabled the study of the 
effects of frustration alone. Structurally ordered, ZnCr 2 4 is a highly frustrated 5 = 3/2 
three-dimensional magnet, in which frustrated tetrahedra are weakly connected at the cor- 
ners (spinel lattice). ZnCr 2 4 is not a spin glass; instead it has a sharp peak in C(T)/T, 
and establishes long range order at about T = 12 K. This indicates that the spin-glass-like 
behavior does not arise from frustration alone. SCGO(p) does display spin glass proper- 
ties. Unfortunately, SCGO(p) contains site disorder (except for p = 1, which has not been 
achieved in practice), and even though the percolation threshold indicates that disorder 
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does not limit the correlation length, a cleaner system is desired. BSZGCO is a frustrated 
5 = 3/2 magnet structurally similar to SCGO, thought to have less disorder, though it still 
displays spin-glass-like properties. This points to frustration and low dimensionality, rather 
than disorder, as being essential for anomalous behavior. On the other hand, we know that 
frustrated two-dimensional systems in which frustrated units are strongly connected (for 
instance, triangular lattice, consisting of edge sharing triangles), order at some finite tem- 
peratures, for all physical spin values. Furthermore, the evidence suggests that even the py- 
rochlore lattice, populated with S > 5/2 spins orders at finite temperature (again g < g c ). 
Thus, the anomalous behavior seen in SCGO and BSZGCO seems to arise at a particular 
critical region in dimensionality and spin value, in combination with frustration and weak 
connectivity (low coordination number). 

We find that at low temperatures the excitation spectrum in BSZGCO follows cu/T 
scaling, and does not have any temperature-independent energy scale. This is a charac- 
teristic behavior of systems in the quantum critical region. At temperatures below 1 1 K, 
the relaxation rate distribution broadens, and below about 2 K its lower cutoff vanishes, 
which may be associated with the spin-glass-like transition. We studied the spin freezing 
at low temperatures, and observed that the onset and polarization of the frozen moment is 
timescale dependent. At timescales r > 0.8 ns, freezing starts around 7 K, and at the lowest 
temperatures reaches about 45% polarization, while for r > 5 ps it develops around 15 K, 
and reaches about 70% polarization. Unfrozen portion of the spin remains in a fluctuating 
state. It is characteristic of spin glasses that the onset and amount of freezing is timescale 
dependent at all timescales, as opposed to conventional magnets, where no significant de- 
pendence exists above some characteristic microscopic timescale. 

We have observed that in BSZGCO the wavevector dependence of the excitation spec- 
trum, is temperature independent, and vanishes as Q — > 0. This supports the no- 
tion that the dynamics is driven by an interaction between spin simplexes, which have no 
net spin (thus have minimum intrasimplex energy). Since intrasimplex spin interactions 
are governed by strong spin-spin interaction, simplex remains stable for all temperatures 
T <C |0 C w|> while the observed low-temperature dynamics is governed by a much weaker 
intersimplex interactions. The low temperature order so induced does not involve individ- 
ual spins, but some composite degrees of freedom, and can therefore be difficult to observe. 
Coherent fluctuations of these composite degrees of freedom at low temperatures charac- 
terize the observed glassy phases. 

As the frozen state has lower continuous symmetry than the Hamiltonian, there must ex- 
ist excitations with gapless excitation spectrum in the long-wavelength limit, the Goldstone 
modes. If these modes involve some composite degrees of freedom, rather than individual 
spins (as in spin waves), it is possible that neutrons cannot directly couple to them, but 
instead only diffuse inelastic scattering is observed. We observed no gap in the excitation 
spectrum, and therefore it might be associated with the Goldstone modes. Further study of 
the frozen state and excitations would benefit from experiments on single crystals, where 
intensity is not spherically averaged. 
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Part III 
Conclusions 
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Chapter 7 
Conclusions 



Our research focused on four different types of frustrated systems, one structurally frus- 
trated: (i) materials exhibiting negative thermal expansion (NTE); and three magnetic frus- 
trated systems: (ii) the triangular antiferromagnet, (Hi) an antiferromagnet on a pyrochlore 
lattice, and (iv) an antiferromagnet on akagome sandwich. Structurally frustrated systems, 
Sc 2 (W0 4 ) 3 and ZrW 2 8 , display a highly unusual uniform negative thermal expansion 
over a very wide temperature range. Here, soft phonon modes drive a system toward a 
phase transition, but the lattice symmetry prevents the transition from occurring [35]. In 
magnetic systems, frustration is expressed through a competition between magnetic inter- 
actions which cannot be simultaneously satisfied, leading to quantum magnetism. This 
results in novel magnetic states, characterized not by a long-range magnetic order and spin 
wave excitations, but rather by quantum order (with short range spin order), and excitations 
defined in terms of composite degrees of freedom. Frustration in magnetic systems is par- 
ticularly easy to visualize, as they are generally composed of frustrated units, triangles or 
tetrahedra, vertices of which are occupied by antiferromagnetically interacting spins. 

Effects of frustration are particularly pronounced when coupled with "weak connec- 
tivity," or "undemonstrative," a term referring to a topology in which order in one part 
of the system does not constrain the rest of the system. Systems in which elements are 
strongly connected usually prevent the emergence of composite degrees of freedom, which 
drive novel states and excitations. In structurally frustrated systems, such as NTE materi- 
als, weak connectivity manifests itself by means of the "openness" of the structure: rigid 
units are connected by their corners in such a way that their rotation does not cost much 
energy. Materials in which rigid units are strongly connected, such as Si0 4 or A10 4 , do 
not exhibit anomalous behavior such as NTE [36]. In magnetic systems, corner sharing 
of frustrated units (triangles or tetrahedra) leads to weak connectivity. When magnetically 
frustrated units share edges, and so are strongly connected, system exhibits properties char- 
acteristic of conventional magnets. This was particularly well illustrated by a progression 
of our magnetic samples. Spins arranged in a network of strongly connected triangles, as 
is found in RbFe(Mo0 4 ) 2 , order at sufficiently low temperature, and their excitations are 
spin waves; just as is the case for conventional magnets. (These effects can be weakened 
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in magnetic systems by lowering the spin length and the number of spatial dimensions.) 
Weakly connected systems, on the other hand, such as ZnCr 2 4 and QS-ferrite, display un- 
usual behavior characterized by quantum order, and dynamics driven by composite degrees 
of freedom. 

All frustrated (and, in particular, weakly connected) systems seem to have some com- 
mon properties, such as an unusually large entropy at low energies, or even a macroscopic 
ground state degeneracy [37]. Furthermore, they seem to be characterized by large spec- 
tral weight (DOS) at energies lower than the strength of interaction. These properties may 
sometimes be used as an identifying fingerprint of frustrated systems. For instance, large 
low-temperature heat capacity may be a signature of dynamics driven by composite degrees 
of freedom. Magnetic systems are, again, well suited for study, as frustration in them can 
be identified by unique susceptibility signature: their inverse susceptibility remains linear 
down to temperatures well below Curie-Weiss temperature. 

As mentioned earlier, a variety of systems can be labeled as frustrated. Therefore, al- 
though they share some common properties, it should not be surprising that there are also 
great differences between them. Generally, similarities between frustrated systems will 
be more pronounced in the temperature region in which individual degrees of freedom re- 
main constrained, and composite degrees of freedom thermally disordered. This can be 
best illustrated in magnetic systems. Here, individual degrees of freedom are spins, and 
below Curie-Weiss temperature, 6 C w> their dynamics is constrained by the other spins in 
a frustrated unit. Interactions between composite degrees of freedom, however, are much 
weaker, so they remain in a disordered state down to a much lower temperature, T c . Be- 
tween these two temperatures, the system can be described as a cooperative paramagnet; 
dynamics of such systems will be similar. For example, dynamics of frustrated magnets 
close to a quantum critical point is characterized by a temperature dependent energy scale: 
a so called co/T scaling. Below T c , differences arising from small perturbations, or other 
system-specific details, will cause diverse behavior. 

Systems studied in this thesis are certainly very diverse. Most obviously, NTE systems 
display interesting structural properties, while the others are magnetic. Even among the 
magnetic systems, variability in dimensionality, spin value, connectivity class (weak or 
strong), etc., cause a variety of low-temperature properties. The rich H-T phase diagram 
of RbFe(Mo0 4 ) 2 is unparalleled by those of other systems, while QS-ferrite, with a large 
©cw/T c ratio and no structural phase transition, seems to be ideally suited for studying 
the quantum critical point. Unique properties of these systems have been discussed in 
concluding sections of their respective chapters. Common to all of them is a presence 
of unsatisfied interactions that drive the systems toward a phase transition, as well as a 
global symmetry precluding it. Sometimes, local interactions eventually win. The system 
relieves frustration and, in the case of magnetic systems, establishes long range magnetic 
order. Such is the case with zinc chromite, which completely relieves frustration by a phase 
transition at T = 12.5 K, far below its Curie- Weiss temperature of over 400 K. Zirconium 
tungstate, on the other hand, has a phase transition at T = 428 K, which only partially 
relieves frustration. Rubidium iron molybdate achieves long-range order at T = 3.5 K, well 
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below Curie- Weiss temperature of about 22 K, though it still remains frustrated. In other 
systems, the global symmetry remains unaltered in the entire temperature range. A lack of a 
phase transition demonstrates this in scandium tungstate and QS-ferrite. In the temperature 
range in which the frustration plays an important role, all of the materials we have studied 
display anomalous and intriguing properties that set them apart from conventional systems. 
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Appendix A 

Negative Thermal Expansion - Details 



A.l Interatomic potentials 

In this section we describe the interatomic potentials that were used for the analysis of 
NTE in ZrW 2 8 and Sc 2 (W0 4 ) 3 . 

A number of two body potentials is available in GULP: Buckingham, 1 Buckingham 
four range, 2 Lennard- Jones, 3 Morse, 4 Harmonic (with higher order, an-harmonic terms), 
Polynomial, General, 5 Spline, etc. In principle, nearly any potential can be adjusted so as 
to reproduce the physical properties of the crystal. The following properties of U(r)\ ro , 
where r is the interatomic distance in the crystal, are relevant for modeling: 

(i) The first derivative of the potential with respect to distance, (dU/dr)\ ro . This deter- 
mines the generalized force on each atom. The vector sum of all such generalized 
forces on a single atom must be zero, as the atoms in a crystal do not move (neglect- 
ing oscillations). This parameter determines the static properties of a crystal; namely, 

Buckingham: exponential repulsion and r~ 6 attraction. See Eq. (A. 3). 

2 Buckingham four range is a complicated potential which assigns a different functional form for four 
different regions: 

U(r) = Aexp(-rfp), for l min < r < l cutl ; 

U(r) = a + air + a 2 r 2 + a 3 r 3 + a 4 r 4 + a 5 r 5 , for / cut i < r < r ; 

U(r) = b + 61 r + b 2 r 2 + b 3 r 3 , for r < r < Z cut2 ; 

U(r) = -C/r 6 , for Z cut2 < r < Z max . 

Here, r is the energy minimum (or at least a stationary point), and the parameters are constrained so that the 
potential and its first derivative are continuous at the boundaries. 
3 Lennard- Jones: r~ 12 repulsion and r~ e attraction. 

4 Morse: exponential attraction and exponential attraction, the two having different ranges. 
5 General: 

tt - A ex P(- r /p) C 
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the structure. [Alternatively, one could use the location of the minimum; i.e. r such 
that U (r ) is at the minimum. This determines the equilibrium distance between two 
free atoms (generally not equal to the interatomic distance in a crystal).] 

(ii) The second derivative: (d 2 £7/dr 2 )| ro . This parameter determines the harmonic dy- 
namic properties of the crystal; namely the phonon modes. Also, a range of proper- 
ties can be calculated from these: elastic constants, static dielectric constants, phonon 
density of states, electric field and gradient, heat capacity, entropy, etc. 

(Hi) The third derivative: (d 3 U/dr 3 )\ ro . This parameter controls a number of anharmonic 
properties, such as the thermal expansion. 6 If free energy, rather than the internal 
energy, is to be minimized, the third derivative proves to be critically important. 

The absolute value of the energy, U(r ), is not important for crystal properties, while the 
derivatives above the third are nearly irrelevant. Thus, one can in general use the following 
potential: 



[Or, if the atoms are fixed in their positions, so that the generalized forces become irrele- 
vant, even the first term can be omitted: K\ = 0.] This potential is available in GULP; it 
is called generalized harmonic. Even though it decouples static from dynamic and anhar- 
monic contributions, this potential is un-physical, and thus not very appealing. 
In our analysis, we used only the so-called Buckingham potential: 



This potential has a physical foundation: (i) exponential repulsion, due to overlap of elec- 
tronic clouds at close distances, and (ii) r~ 6 attraction, due to weak induced dipole-dipole 
interaction. These functional forms are correct in the appropriate limits, making Bucking- 
ham potential particularly appealing. Having three adjustable parameters, A, p, and C, 
the first three derivatives can attain any value at an equilibrium distance, r . Specifically, 
the advantage of the Buckingham potential over a generalized harmonic potential is that 
one triplet of interatomic parameters (that is, {^4, p, C}), can be used for many interatomic 
distances r ; in contrast to that, generalized harmonic parameters ({ki, k 2 , k 3 }) must be 
determined separately for each interatomic distance. 

The three body potentials available in GULP are: three body harmonic (with or without 
exponential decay), Axilrod-Teller three body, Stillinger- Weber three body, Urey-Bradley 
three body, etc. We found the simple three body harmonic potential to be adequate for our 
work: 7 



^gen.harm.( r ) = M r ~ r o) + y ( r ~ r of + ~q ( r ~ r of '• 



(A.2) 



U huck (r) = Aexp(-r/p) 




(A3) 




(A.5) 



A crystal described by purely harmonic potentials has no thermal expansion. 
Generalized three body harmonic potential has the form: 



tWee = f (0 - e ) 2 + 1 (e - e f + |i (e - e )\ 



(A.4) 
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It is reasonable to expect the three body potentials to be important for covalent bonds; much 
less so for ionic interactions. The contribution of three body potentials to the cohesion 
energy is much smaller than that of two body potentials. 

Both the two- and three-body potentials are short range potentials. GULP defines a 
cutoff parameter, Z cuto ff, such that by definition U(r > /cutoff) = 0. This parameter is 
usually of the order of a few Angstroms. 

Atoms in GULP can be charged, and interact via long-range Coulomb potential. GULP 
uses a sophisticated algorithm combining direct and reciprocal space summation, to ac- 
count for this long-range interaction. We found charges to be important for our understand- 
ing of crystal properties, and assigning formal charges to atoms (i.e., Zr 4+ , Sc 3+ , W 6+ , and 
+2 ) to give good results. 

Next, GULP allows atoms to be represented by a shell atomic model, in which the 
massless electronic cloud (referred to as shell) is considered a dynamic entity separate 
from the nucleus (called core) which contains the entire mass of the atom. Both the shell 
and the core carry charge, and are connected by a spring potential: 

TJ - Spring ft /A fi N 

u spring — 2 "core-shell- 

Here, d core -sheii is the distance between the centers of core and shell. We found it useful 
to represent oxygen atoms by a shell model, while cations were well represented as point 
particles. 

We found the potentials described above to be sufficient for description of our systems. 
For reference, Table A. 5 shows the refined interatomic potentials that we used for ZrW 2 8 . 
GULP, however, allows additional potentials, which we have not used; for instance, four 
body potentials (torsion), dipolar shells, out-of-plane potential, etc. 

A.2 Fitting procedure 

In order to determine the interatomic potentials, we defined twenty parameters, pf, i = 
1 to 20; the definition for these parameters are given in Table A. 6. We used parameters p\ 
to P13 for the structure fit. 8 Here, we describe the fitting procedure: 

Given structure and interatomic potentials, GULP can calculate the internal energy of a 

We found, however, k% and k^ to be unimportant for our work, and have set them to zero. 

8 Additional parameters were introduced later for fit to phonon measurements; they were not relevant for 
the structure fit. 
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T(K) a(A) Space group 



0.3 


9 1823(4) 


P2i3 


3.3 


9 1810(5) 


P2i3 


12 


9 1 810(5) 


P2n3 


20 


9.1797(6) 


P2i3 


80 


9 1758(5) 


P2i 3 


150 


9.1686(5) 


P2 X 3 


220 


9.1631(5) 


P2 X 3 


293 


9.1568(3) 


P2 : 3 


333 


9.1530(3) 


P2 : 3 


373 


9.1500(4) 


P2i3 


423 


9.1416(5) 




483 


9.1371(5) 


Pa3 


583 


9.1318(4) 


Pa3 


693 


9.1318(4) 


Pa3 



Table A.l: Cell dimensions of ZrW 2 8 as a function of temperature. Adapted from 
Ref. [50]. 



Atom 


Occupancy 


X 


y 


z 


U eq (A 2 ) 




T=293 K, anisotropic refinement model. 




Zrl 


1.0 


0.0004(3) 


0.0004(3) 


0.0004(3) 


0.008 


Wl 


1.0 


0.3409(3) 


0.3409(3) 


0.3409(3) 


0.010 


W2 


1.0 


0.6009(3) 


0.6009(3) 


0.6009(3) 


0.006 


01 


1.0 


0.0529(3) 


-0.2069(3) 


-0.0619(4) 


0.020 


02 


1.0 


0.0697(4) 


-0.0575(3) 


0.2132(3) 


0.017 


03 


1.0 


0.4941(4) 


0.4941(4) 


0.4941(4) 


0.022 


04 


1.0 


0.2322(3) 


0.2322(3) 


0.2322(3) 


0.035 




T=483 K, anisotropic refinement model. 




Zrl 


1.0 


0.0000(0) 


0.0000(0) 


0.0000(0) 


0.014 


Wl 


0.5 


0.3349(5) 


0.3349(5) 


0.3349(5) 


0.011 


W2 


0.5 


0.6035(5) 


0.6035(5) 


0.6035(5) 


0.011 


01 


1.0 


0.0549(3) 


-0.2089(2) 


-0.0671(3) 


0.036 


03 


0.5 


0.5055(0) 


0.5055(0) 


0.5055(0) 


0.044 


04 


0.5 


0.2322(4) 


0.2322(4) 


0.2322(4) 


0.064 



Table A. 2: Fractional atomic coordinates of ZrW 2 8 . Positions refined in space group 
P2 X 3 (198) for T = 293 K, and Pa3 (205) for T = 483 K. Adapted from Ref. [50]. 
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T(K) a(A) 6(A) c(A) V(A 3 ) 



10.0 
50.0 
100.0 
150.0 
200.0 
250.0 
300.0 
350.0 
400.0 
450.0 



9.6883 
9.6875 
9.6843 
9.6812 
9.6784 
9.6753 
9.6720 
9.6689 
9.6659 
9.6620 



13.297 
13.299 
13.303 
13.307 
13.310 
13.314 
13.318 
13.322 
13.325 
13.329 



9.5956 
9.5947 
9.5911 
9.5881 
9.5852 
9.5822 
9.5795 
9.5773 
9.5749 
9.5725 



1236.2 
1236.1 
1235.6 
1235.2 
1234.8 
1234.3 
1234.0 
1233.6 
1233.2 
1232.8 



Table A. 3: Cell dimensions of Sc 2 (W0 4 ) 2 as a function of temperature. Adapted from 
Ref. [31]. 

crystal. In our model of scandium tungstate, this equals: 

Cq-o 



U = £ ^o-oexp(^^) 
{ to } V P°-° J 

a f—rw-o\ Cw-o 
+ y Av-o exp 

^ Pw-o ' 
{W-O} KW u 

a f—rsc-o\ Csc-o 
y A Sc -o exp 

{Sc-O} 

&O- 



r-w-o 



rs c -o 



+ ^2 2 (^o-sc-o _ (#o)o-sc-o) 

{O-Sc-O} 

(*o)o-w-o)' 



{O-W-O} 

(^spring)c-s ,2 

{c-s} 

+ _Ly^i. (a.7) 

Here, {O-O}, {W-O}, and {Sc-O}, represent atomic pairs; and r -o> ?"w-o> an d fsc-o 
the distances between them, respectively. The summation is carried only over pairs for 
which r < Z cu toff- Next, {O-Sc-O} and {O-W-O} represent atomic triples (belonging 
to Sc0 6 octahedra, and W0 4 tetrahedra); #o-s c -o and 6> -w-o the angles between the 
bonds; while (#o)o-s c -o an d (#o)o-w-o are the angles corresponding to minimum ener- 
gies ((#o)o-sc-o = 90° for octahedra, and (0 o )o-w-o = 109.47° for tetrahedra). The shell 
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Atom 


X 


y 


z 


< o 100 






T=10K 






Scl 


0.46534(21) 


0.38045(15) 


0.24894(28) 


1.01(5) 


Wl 


0.25000 


0.00000 


0.4751(7) 


0.55(18) 


W2 


0.1191(5) 


0.35551(35) 


0.3969(5) 


-0.25(11) 


01 


0.0926(4) 


0.14109(32) 


0.0781(4) 


0.40(11) 


02 


0.1267(4) 


0.06006(29) 


0.3692(4) 


0.33(11) 


03 


0.0127(4) 


0.26289(32) 


0.3233(4) 


0.57(10) 


04 


0.3318(4) 


0.40929(30) 


0.0810(4) 


0.77(10) 


05 


0.0759(4) 


0.47350(31) 


0.3229(4) 


0.47(10) 


06 


0.2923(5) 


0.32959(28) 


0.3638(4) 


0.61(10) 






T=300 K 






Scl 


0.46688(22) 


0.38110(17) 


0.24937(31) 


1.33(5) 


Wl 


0.25000 


0.00000 


0.4743(8) 


1.29(22) 


W2 


0.1188(6) 


0.3557(4) 


0.3949(5) 


0.68(14) 


01 


0.0917(5) 


0.14199(32) 


0.0770(5) 


1.73(13) 


02 


0.1302(5) 


0.06284(39) 


0.3693(5) 


2.09(14) 


03 


0.0118(5) 


0.26323(38) 


0.3215(5) 


2.00(13) 


04 


0.3358(5) 


0.41072(34) 


0.0787(5) 


1.81(13) 


05 


0.0733(5) 


0.47255(36) 


0.3226(5) 


1.71(13) 


06 


0.2923(5) 


0.33169(35) 


0.3608(5) 


1.64(12) 






T=450 K 






Scl 


0.46740(23) 


0.38146(18) 


0.24892(34) 


1.58(5) 


Wl 


0.25000 


0.00000 


0.4744(9) 


1.54(23) 


W2 


0.1179(6) 


0.3555(4) 


0.3941(6) 


0.91(15) 


01 


0.0916(5) 


0.1422(4) 


0.0750(5) 


2.36(15) 


02 


0.1323(6) 


0.0636(4) 


0.3689(6) 


2.78(15) 


03 


0.0119(6) 


0.2623(4) 


0.3209(6) 


2.72(15) 


04 


0.3373(6) 


0.4119(4) 


0.0779(6) 


2.76(16) 


05 


0.0699(5) 


0.4720(4) 


0.3228(6) 


2.51(15) 


06 


0.2917(6) 


0.3335(4) 


0.3598(5) 


2.17(13) 



Table A. 4: Fractional atomic coordinates of Sc 2 (W0 4 ) 2 . Positions refined in space group 
Pnca (60:cab). Adapted from Ref. [31]. 
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Zr(core) — O(shell) 
W(core)— O(shell) 
O(shell)— O(shell) 



Buckingham potentials 
A(eV) p(A) 
9.0 x 10 6 0.14 
1305.22 0.375 
22764. 0.149 



C(eV/A 6 ) 





27.879 





Three-body harmonic potentials 






fc(eV/rad 2 ) 


0o(°) 


O(core)- 


— Zr(core) — O(core) 0.4 


90.0 


O(core)- 


-W(core)— O(core) 0.5 


109.47 



Spring potential 



^springC^V/A 2 ) 



O(core)— O(shell) 74.92 



Charges 



Zr(core) 
W(core) 
O(core) 
O(shell) 



Q(e) 
+4.0 
+6.0 

+0.84819 
-2.84819 



Table A. 5: Interatomic potentials for ZrW 2 8 . Cutoff distance for short-range potentials: 

o 

Cutoff = 2.30 A. 
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Two body O — O Buckingham potential 
Pi ^o-o [eV] 

P2 po-o [A] 

Ps C -o[eVA 6 ] 

Two body W — O Buckingham potential 
Pa A w _ [eV] 

Pu Pw-o [A] 

Pe C W -o[eVA 6 ] 

Two body Sc — O Buckingham potential 
P7 A Sc -o [eV] 

P8 PSc-0 [A] 

p 9 C Sc - [eVA 6 ] 

Three body O — Sc — O harmonic potential 
Pio ^o-sc-o [eVrad" 1 ] 

Pn (^o)o-w-o [A] 

Three body O — W — O harmonic potential 
P12 ko-w-o [eV rad^ 1 ] 

Pn (0 o )o-w-o [A] 

Oxygen charge distribution (x) 

(Qcore — ~\~X, Qshell — 2 X) 

Pu x [e] 

Oxygen spring constants 

PlS tO P20 (fespring)oi tO (fe spr ing)o6 [gV A" 2 ] 



Table A. 6: Definition of {p^ parameters used for the fit of Sc 2 (W0 4 ) 3 interatomic po- 
tentials. Parameters p x to p 13 were used for structure fits. Parameters p u to p 2 o were 
introduced later, for density-of- states fits. 
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atomic model is used only for oxygen atoms, and {c-s} refers to core-shell oxygen pairs 
(belonging to the same atom); d c _ s is the distance between the core and its shell. Finally, 
the last summation is carried out over all charged pairs, indexed with {ij}. This summation 
is not trivial, due to the long-range nature of electrostatic potential. Here, GULP imple- 
ments a sophisticated algorithm, using "Ewald summation, in which the inverse distance is 
rewritten as its Laplace transform and then split into two rapidly convergent series, one in 
reciprocal space, and one in real space" [92]. 

Using the "optimize" keyword, GULP can vary the unit cell size and fractional atomic 
positions in order to minimize the internal energy (See Eq. A.7). We refer to this structure 
as relaxed; obviously, it is a function of interatomic potentials: 

({r d }{a, 6, c}) = ({r d }{o, b, c}) ({ Pi }) . (A.8) 

V / relaxed V / 

(Notation: ({r d }{a, b, c}) refers to crystal structure; ({a, 6, c} is the unit cell size; while 
{r d }) refers to atomic fractional coordinates.) Next, to measure how close the relaxed 
structure is to actual crystal structure, we define the following goodness of fit function: 

x 2 ({pi}) = Yl |( rex p) d ~ ( T ™i.)d + 

d 

const, x ((a e xp - a r ei.) 2 + (b exp - 6 re i.) 2 + (c exp - c re i.) 2 ) . (A.9) 

Constant const, is adjusted as needed during the fit. The structure fit now consists of 
varying parameters {p^ in order to minimize \ 2 - 

GULP receives all its instructions via an input file, and can also store the output into a 
file. We coded an IDL driver function, "sumof sq," which takes {p^ as an input parameter, 
and then (i) writes a properly formatted GULP input file, inserting {pi) for interatomic 
parameters; (ii) calls GULP for structure relaxation, 9 and has the output saved into a file; 
(Hi) reads the structure from the output file, and then calculates and returns the value of 
X 2 . Using IDL built-in optimization functions, we minimize "sumof sq," and arrive at the 
optimal interatomic parameters. 

We found some of the parameters to be very important for the stability of the structure, 
while others much less so. Therefore, we often performed fits using a reduced, rather than 
the entire, parameter space. 

A.3 Phonon eigenvector properties 

We start with equation (6.22) from Dove [48]: 

D aP (jj', k) = 1 ]T $ Q/3 ( Jjl ) exp (ik ■ [r(j'l') - r(jO)]) , (A.10) 
\/ m j m j' // v / 

9 A single "optimize" call to GULP takes about 5-10 seconds. 
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where j is atom index, and / unit cell index. Substituting k with k + K we get: 

B °"W' k + K » = ^^^(M') eXp ( i(k + K) - |r(iT) - r(i0)1 ) 

= 1 V $ J ii'} e ^Wn-r(jO)] JK-[r(j'l')-r(jO)] (A U) 

V"w v v 0/ v ' - ' 

The underbraced factor can be developed as follows: 

e iK-[p(j'l')-p(jO)] = e iK.[(p(j'0)+H,,)-p(jO)] 
e iK-R ; , e iK.[r(j'0)-r(jO)] 



=1 

o iK-[p(j'0)-p(j0)] 



(A.12) 
Thus (A.l 1) reduces to: 

A* k + K) = e «-[rO'0)-rO0)] x J- ^ e*WW* 

= e iK WV- r ^D a p(jj',k). (A.13) 
Equation (6.19) from Dove gives: 

^ 2 (k^)U(j,k,.) = Y,**>(oi) ■U(j',k jI /)exp(ik.[r(j'0-r(jO)]) 

= 5^ u o -, . k i v)v m 3 m j' 

" v ' 

=D(jj'M) 

= ^v^W^0'i'> k )U(i',k^). (A.14) 

i' 

Equation (A.14) is valid for any k; therefore we can replace k by k + K: 

m j o; 2 (k + K,^)U(j,k + K,^) = y/mjmfD(jj', k + K)U(j', k + K, v). 
v v ' i' 

=OJ 2 (k,i/) J 

(A.15) 

Using equation (A.13) right-hand side can be written as: 

u; 2 (k,i/)U(j,k + K» = ^V7^/J(jj',k) e lK -[ r ^- r ^U(j',k + K^) 



m 3 



e- iK <^ £ ^/frqmpDUf, k) e iK ' r(/0) U(j', k + K, u); 
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mj uj 2 (k, u) U(j, k + K, z,)e iK '^ > = V v ^nJrniD(jj', k) U(j', k + K, . 

v ' s, ' 

(A.16) 

If displacement vectors U are properly normalized, so that equation (6.27) from Dove is 
satisfied, then from equations (A. 14) and (A.16) it follows that: 

U(j,k + K» = [U(i,k^) e - iK ^°)] e ^), (A.17) 

where (/>(k, v) is some arbitrary phase factor, independent of atom index j. Since eigenvec- 
tors e(j, k, v) are related to displacement vectors U(j, k, v) by Dove (6.21): 

e(j, k, v) = ^m~U(j, k, v), (A. 18) 

it follows that eigenvectors satisfy: 

e(j, k + K, v) = e(j, k, v ) e -^m e ^(^)_ (A. 19) 

For the special case of an atom at the origin of the unit cell, r(j 0) = 0, it follows that: 

e(j ,k + K» =e(j ,k,iy)e i ^\ (A.20) 

Once (j)(k, v) is determined (either from (A. 19) or (A.20)), equation (A. 19) can be used to 
relate the eigenvectors at points k + K and k. 10 



10 Equation (A. 19) did not hold for eigenvectors calculated by GULP. This prompted a continued investi- 
gation, yielding a conclusion that the eigenvectors calculated by GULP v. 1 .2 are incorrect. 
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Appendix B 
RbFe(Mo0 4 ) 2 - Details 

B.l X-ray structure report 

Following pages reproduce 1 the Crystal Structure Report on RbFe(MoO '4)2, prepared by 
the X-Ray Cry stallo graphic Laboratory at the University of Minnesota. The conclusions 
drawn from this report are presented in Sec. 5.2 (p. 71). The content of Tables 3 (bond 
lengths and angles) and 5 (torsion angles) have been shortened, as the data contained 
therein can easily be calculated from the structure information. In the original report, 
Table 3 had 151 entries (27 lengths and 124 angles), while Table 5 had 250 entries. 



1 Original report has been written in RTF format (file "01115B.RTF")- The version presented here is 
written in I4Tj5C, and has been slightly modified to conform to the printed format of this thesis. 
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CRYSTAL STRUCTURE REPORT 
RbFe(Mo0 4 ) 2 

Report prepared for: 
G. Gasparovic and Prof C. Broholm - Johns Hopkins University 
18 April 2001 Victor G. Young, Jr. 




X-Ray Crystallographic Laboratory 
Department of Chemistry 
University of Minnesota 

207 Pleasant St. S.E. 
Minneapolis, MN 55455 
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Data collection/ 

A crystal (approximate dimensions 0.094 x 0.074 x 0.015 mm 3 ) was placed onto the tip of 
a 0.1 mm diameter glass capillary and mounted on a Bruker SMART system for a data 
collection at 173(2) K. A preliminary set of cell constants was calculated from reflections 
harvested from three sets of 20 frames. These initial sets of frames were oriented such that 
orthogonal wedges of reciprocal space were surveyed. This produced initial orientation 
matrices determined from 108 reflections. The data collection was carried out using 
MoKa radiation (graphite monochromator) with a frame time of 30 seconds and a 
detector distance of 4.9 cm. A randomly oriented region of reciprocal space was surveyed 
to the extent of 1.5 hemispheres and to a resolution of 0.84 A. Three major sections of 
frames were collected with 0.30° steps in u at 3 different settings and a detector position 
of —28° in 26. The intensity data were corrected for absorption and decay (SADABS). 1 
Final cell constants were calculated from 1272 strong reflections from the actual data 
collection after integration (SAINT 6.01, 1999). 2 Please refer to Table 1 for additional 
crystal and refinement information. 

Structure solution and refinement 

The structure was solved using SHELXS-86 3 and refined using SHELXL-97. 3 The space 
group P 3 ml was determined based on systematic absences and intensity statistics. A 
direct-methods solution was calculated which provided most non-hydrogen atoms from 
the E-map. Full-matrix least squares / difference Fourier cycles were performed which 
located the remaining non-hydrogen atoms. All non-hydrogen atoms were refined with 
anisotropic displacement parameters unless stated otherwise. All hydrogen atoms were 
placed in ideal positions and refined as riding atoms with relative isotropic displacement 
parameters. The final full matrix least squares refinement converged to Rl = 0.021 1 and 
wR2 = 0.0548 (F 2 , all data). 

Structure description 

The structure is found as proposed. All atoms are on special positions. The Mo atom is 
bonded to two oxygens and is nearly tetrahedral. These Mo04 subunits are linked by Fe 
atoms to form sheets. The Rb atoms reside at the cell origin. These have contacts with 
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oxygen in the 3.0 to 3.3 Arange. There is one formula unit per unit cell. The space group 
and atomic positions have been checked for possible higher symmetry with the IUCr 
CHECKCIF program. 4 A packing drawing (cover) and the contents of the asymmetric 
unit are provided. 

Data collection and structure solution were conducted at the X-Ray Crystallographic 
Laboratory, 160 Kolthoff Hall, Department of Chemistry, University of Minnesota. All 
calculations were performed using SGI INDY R4400-SC or Pentium computers using the 
current SHELXTL suite of programs. All publications arising from this report MUST 
either l)include Victor G. Young, Jr. as a coauthor or 2)acknowledge Victor G. Young, Jr. 
and the X-Ray Crystallographic Laboratory. Drafts must be supplied to the coauthor prior 
to journal submission. Additional drawings will be supplied upon request. 



1 An empirical correction for absorption anisotropy, R. Blessing, Acta Cryst. A51, 33 - 

38(1995). 

2 SAINT V6.1, Bruker Analytical X-Ray Systems, Madison, WI. 

3 SHELXTL-Plus V5.10, Bruker Analytical X-Ray Systems, Madison, WI. 1272 

4 CHECKCIF, International Union of Crystallography, 2000. 



143 



RBM) 



012) 



Mod) 



0(11 



Fed) 



Table 1 . Crystal data and structure refinement for 01 1 15b. 



Identification code 
Empirical formula 
Formula weight 
Temperature 
Wavelength 
Crystal system 
Space group 
Unit cell dimensions 



Volume 



01115b 

Fe Mo 2 8 Rb 

461.20 

173(2) K 

0.71073 L 

Trigonal 

p3 m l 

a = 5.6923(6) A, 
b = 5.6923(6) A, 
c = 7.4572(11) A, 
209.26(4) A 3 



a = 90°. 
/3 = 90°. 
7 = 120 c 
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Z 

Density (calculated) 
Absorption coefficient 
F(000) 

Crystal habit and color 
Crystal size 

Theta range for data collection 
Index ranges 
Reflections collected 
Independent reflections 
Observed Reflections 
Completeness to theta = 25.02° 
Absorption correction 
Max. and min. transmission 
Refinement method 
Data / restraints / parameters 
Goodness-of-fit on F2 
Final R indices [I>2sigma(I)] 
R indices (all data) 
Extinction coefficient 
Largest diff. peak and hole 



1 

3.660 Mg/m 3 
10.442 mm" 1 
211 

Yellow, Plate 

0.094 x 0.074 x 0.015 mm" 1 
2.73 to 25.02°. 

-5 = h = 0,0 = k = 6,0 = l = 8 
165 

165 [R(int) = 0.0420] 
165 

100.0 % 

Face-indexed 

0.80967 and 0.43738 

Full-matrix least- squares on F 2 

165/0/18 

1.006 

Rl =0.0211, wR2 = 0.0548 
Rl =0.0211, wR2 = 0.0548 
0.061(7) 

0.523 and -0.786 e.A- 3 



Table 2. Atomic coordinates (x 10 ) and equivalent isotropic displacement parameters 
(A 2 x 10 3 ) for 01115b. U eq is defined as one third of the trace of the orthogonalized Uy 
tensor. 



x y z U(eq) 

Rbl 18(17 

Fel 5000 9(1) 

Mol 3333 6667 2681(1) 8(1) 

01 1621(3) 3243(6) 3412(4) 21(1) 

02 3333 6667 368(8) 19(1) 
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Table 3. Bond lengths [A] and angles [°] for 01 1 15b. 



Rbl-Ol 
Fel-Ol 



3.005(3) 
1.990(3) 



Rbl-02 
Mo 1-02 



3.2979(6) 
1.725(6) 



Ol#5-Fel-Ol#10 
Ol#10-Fel-Ol#3 
01#14-Mol-01#15 



180.0 

91.81(13) 

111.01(10) 



01#5-Fel-01#3 
02-Mol-01#14 
Mol-Ol-Fel 



88.19(13) 

107.89(10) 

161.35(19) 



Symmetry transformations used to generate equivalent atoms: 
#1 -x,-y,-z #2 x-y,x,-z #3 -x+y,-x,z #4 y,-x+y,-z 
#5 -y,x-y,z #6 -x,-y+l,-z #7x,y-l,z #8 -x+l,-y+l,-z 
#9x-l,y-l,z #10 y,-x+y,-z+l #1 1 x-y,x,-z+l 
#12 -x,-y,-z+l #13 x,y,z+l #14 -x+y,-x+l,z 
#15 -y+l,x-y+l,z #16x,y+l,z #17 x+l,y+l,z 

Table 4. Anisotropic displacement parameters (A 2 x 10 3 ) for 01115b. The anisotropic dis- 
placement factor exponent takes the form: 

-2vr 2 [h 2 a* 2 U n + ... + 2 h k a* b* U 12 ] 

Un U22 U33 U23 U13 U12 

Rbl 22(1) 22(1) 12(1) 11(1) 

Fel 7(1) 7(1) 13(1) 3(1) 

Mol 8(1) 8(1) 10(1) 4(1) 

01 22(1) 12(2) 25(2) 5(1) 3(1) 6(1) 

02 20(2) 20(2) 17(3) 10(1) 
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Table 5. Torsion angles [°] for 01 1 15b. 



01#1-Rbl-Fel-01#5 

01-Rbl-Fel-01#5 

01#2-Rbl-Fel-01#5 



-60.0 
120.0 
60.0 



02#8-Rbl-02-Rbl#17 0.0 
02#9-Rbl-02-Rbl#17 180.0 



Symmetry transformations used to generate equivalent atoms: 
#1 -x,-y,-z #2 x-y,x,-z #3 -x+y,-x,z #4 y,-x+y,-z 
#5 -y,x-y,z #6 -x,-y+l,-z #7x,y-l,z #8 -x+l,-y+l,-z 
#9x-l,y-l,z #10 y,-x+y,-z+l #1 1 x-y,x,-z+l 
#12 -x,-y,-z+l #13x,y,z+l #14 -x+y,-x+l,z 
#15 -y+l,x-y+l,z #16x,y+l,z #17x+l,y+l,z 
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Atom Sym. Mult, x y z 100 x U\ 



Rb 


-1 


1 











0.54(6) 


Fe 


-3 


1 








1/2 


0.26(5) 


Mo 


3 


2 


1/3 


2/3 


0.2690(3) 


0.55(5) 





1 


6 


0.2192(3) 


0.3238(4) 


0.3435(3) 


0.70(3) 





3 


2 


1/3 


2/3 


0.0370(3) 


0.72(6) 



Table B.l: Structure of RbFe(Mo0 4 ) 2 at T = 1.7 K. Lattice constants: a = 5.5971(1) A, 
c = 7.4436(2) A. Space group P3. (R p = 6.88%, R wp = 8.67%, %? educed = 3.135.) 



B.2 Nuclear structure refinement (powder) 

In this section, we report the details regarding the analysis of the neutron diffraction 
measurements on powder RbFe(Mo0 4 ) 2 . 2 Structure has been refined in P3 space group 
at three different temperatures, using GSAS package (POWPREFF program). Figures 5.4, 
5.5, and 5.6 show the spectra, and Tables B.l, B.2, and B.3 the derived structures at T = 
1.7 K, T = 10 K, and T = 180 K, respectively. 

The results presented here are derived from the second diffraction experiment, con- 
ducted on July 25-27, 2001 (the ceramic sample). The measurements from the first ex- 
periment, conducted April 24-26, 2001, and judged to be unsatisfactory due to strain from 
grinding, 3 are not reported. Both experiments were done on BT-1 diffractometer at NCNR. 

Initially, isotropic thermal motion constants, U\ so , were used for refinement at all tem- 
peratures; values so obtained are reported in Tables B.l, B.2, and B.3. Thereafter, anisotropic 
thermal motion tensors, Uij, were used for T = 10 K and T = 180 K data; values obtained 
are reported in Tables B.4, and B.5. For T = 1.7 K data, the introduction of an anisotropic 
tensor did not significantly improve the fit, probably due to (i) low temperature (thus small 
thermal motion); and (ii) the presence of the magnetic peaks. 4 (U iso = (Un + U 2 2 + U ss )/2) 



B.3 Magnetic structure (powder) 

We used the low temperature (T = 1.7 K) powder diffraction measurement to deter- 
mine the magnitude of the static magnetic moment in RbFe(Mo0 4 ) 2 . This analysis was 
completed in September 2002; by that time we have determined the magnetic structure 
from the single crystal measurements. This was necessary because both the number and 
the intensity of magnetic peaks in powder diffraction spectrum were insufficient to do a 
direct magnetic structure refinement. 



2 See Sec. 5.4 (p. 77). 

3 See Sec. 5.4 (p. 81) for details. 

4 See Sec. B.3 for magnetic analysis. 



148 



Atom 


Sym. 


Mult. 


X 


y 


z 


100 x U iso 


Rb 


-1 


1 











0.49(5) 


Fe 


-3 


1 








1/2 


0.24(4) 


Mo 


3 


2 


1/3 


2/3 


0.2689(2) 


0.46(4) 





1 


6 


0.2186(2) 


0.3232(3) 


0.3428(2) 


0.64(2) 





3 


2 


1/3 


2/3 


0.0371(2) 


0.67(5) 



Table B.2: Structure of RbFe(Mo0 4 ) 2 at T = 10 K. Lattice constants: a = 5.5980(1) A, 
c = 7.4428(1) A. Space group P3. (R p = 5.46%, R wp = 7.00%, Educed = 1-716.) 



Atom 


Sym. 


Mult. 


X 


y 


z 


100 x U iso 


Rb 


-1 


1 











1.52(8) 


Fe 


-3 


1 








1/2 


0.61(5) 


Mo 


3 


2 


1/3 


2/3 


0.2697(3) 


0.79(5) 


O 


1 


6 


0.2022(3) 


0.3222(4) 


0.3429(3) 


1.46(3) 


O 


3 


2 


1/3 


2/3 


0.0377(3) 


1.60(7) 



Table B.3: Structure of RbFe(Mo0 4 ) 2 at T = 180 K. Lattice constants: a = 5.6338(1) A, 
c = 7.4620(1) A. Space group P3. (Rp = 6.04%, R wp = 7.46%, Educed = 1-175.) 



Atom 


100 x U n 


100 x Um 


: 100 x C/ 33 


100 x U 12 


100 x U 13 


100 x U 23 


Rb 


0.48(8) 


0.48(8) 


0.57(10) 


0.24(4) 








Fe 


0.30(6) 


0.30(6) 


0.26(7) 


0.15(3) 








Mo 


0.39(5) 


0.39(5) 


0.75(10) 


0.19(3) 








O 


0.59(6) 


0.59(6) 


0.78(5) 


0.24(6) 


0.30(6) 


0.26(8) 


O 


0.92(6) 


0.92(6) 


0.14(10) 


0.46(3) 









Table B.4: Atomic thermal motion tensor, Uu, for RbFe(Mo0 4 ) 2 at T = 10 K. 



Atom 


100 x U u 


100 x U 22 


100 x U 33 


100 x U 12 


100 x U 13 


100 x U 23 


Rb 


1.78(11) 


1.78(10) 


1.21(14) 


0.89(5) 








Fe 


0.51(8) 


0.51(8) 


0.80(10) 


0.25(4) 








Mo 


0.84(7) 


0.84(7) 


1.21(13) 


0.42(4) 








O 


1.54(11) 


0.91(8) 


1.91(9) 


0.39(11) 


0.32(10) 


0.51(11) 


O 


2.04(10) 


2.04(10) 


0.60(16) 


1.02(5) 









Table B.5: Atomic thermal motion tensor, Uu, for RbFe(Mo0 4 ) 2 at T = 180 K. 
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1 x y z 

2 -y x-y z 

3 y-x -x z 

4 -x -y -z 

5 y -x+y -z 

6 -y+x x -z 



Table B.6: The equivalent positions for space group P3. Description: centric, primitive 
trigonal. Laue symmetry: 3. Multiplicity of a general site: 6. 

The calculation was done using POWPREFF program (part of GSAS package). Since 
POWPREFF cannot be used for non-periodic crystal structures, and here the magnetic 
structure is not periodic along the c-axis, we used the following approximation. Within 
each layer, a 120° structure is easily represented by a magnetic unit cell that contains three 
spins (thus three times larger than the nuclear cell), at and angle of 120° to each other. In 
the neighboring layer, spins are rotated by <f> k, 163° ; 5 thus the structure is either incom- 
mensurate, or has a very large unit cell in the c direction. Assuming, however, the rotation 
angle to be ^ = 162°; 6 the structure becomes periodic, repeating itself after 20 layers. 
Thus, we enlarged the magnetic unit cell 60 times (3 times in-plane, times 20 layers), and 
manually fixed in it 60 spins of unit magnitude. This approximation causes the magnetic 
peaks to be shifted by a small amount; too small to be relevant in this case. 

Next, we refined the data using two powder phases: (i) nuclear phase, Pi, and (ii) 
magnetic phase, P 2 . As both the nuclear and magnetic structures have already been deter- 
mined, there were only two fit parameters: the fraction of each phase. These were found to 
be (units are arbitrary): 

Pi = 27671. 

P 2 = 5981.5 (B.l) 

The magnetic moment in the phase P 2 was set to unity (unity in the sense (gS) = 1). 
Calculated intensity of nuclear peaks is proportional to the fraction of phase Pi. 

I N = c 1 xP 1 ; (B.l) 

while that of magnetic peaks is: 

I M = c 2 x 60P 2 {gS) 2 . (B.3) 

where ci and c 2 are complicated functions of the structure, and are calculated by POW- 
PREFF. For our real, physical sample, Pi = 60P 2 (as the magnetic unit cell is 60 times 

5 (f> = 360° xA;A« 0.453. 

6 This would correspond to a magnetic Bragg peak at Q = (1/3 1 /3 0.45), i.e. A = 0.45. 
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larger than the nuclear one), and (S) is the actual value of the static moment; thus: 

Ci X P x Ci 



(-) 



I M Jobs c 2 x 60P 2 (gS) 2 c 2 x {gS) 2 
while the model sample has {gS) = 1, and so: 

Ci X Pi Ci X Pi 



(B.4) 



(-) 



A//caic c 2 x 60P 2 {gS) 2 c 2 x60P 2 ' ) 

When a fit is good, the calculated ratio of intensity of nuclear and magnetic peaks is the 
same as the observed one: 

(£) =(r L ) • (M) 

\ I M / obs vijf/ calc 

Follows: 

Ci Ci X Pi 



c 2 x(^) 2 ' c 2 x60P 2 

<*> - ^ 

Since in RbFe(Mo0 4 ) 2 the magnetic ion is Pe 3+ , we have g = 2, 7 and so: 



Substituting values obtained from the fit, Eq. (B.l), we get: 

(S) = 1.80. (B.9) 

This is about 72% of the available moment, which for Fe 3+ is equal to S = 5/2. As 
discussed in Sec. 5.4, the remaining fraction of the spin remains in fluctuating state. 

Unfortunately, POWPREFF does not calculate the uncertainties in fitted fractions, which 
are needed to determine the uncertainty in (S). Luckily, it does calculate x 2 , the goodness 
of the fit: 



x 2 



N - 1 



N .calc _ , .exp 2 



which we can utilize as follows. The fitting procedure consists of minimizing x 2 , and here 
the minimum value is xLin = 1-5067, attained for P 2 = 5981.5. Rather than fit, we can 
now fix the value of P 2 , and calculate x 2 \ we can then regard x 2 a function of P 2 : 

X 2 = X 2 {P2). (B.ll) 



7 Because the orbital moment is quenched: L = 0; therefore J = S, and g = 2.0023 s=s 2 (see Eq. (3.21) 
on p. 22). Ground state term for Fe 3+ ion is 6 S 5 / 2 , basic electron configuration 3d 5 , and p = 5.91 [93]. 
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Calculating \ 2 f° r 28 different values of P 2 , and fitting to a polynomial, we obtained the 
following dependence: 

X 2 (P 2 ) = 2.4864 - 3.2754 x 1(T 4 P 2 + 2.7379 x KT 8 P 2 

= 1.5067 + 2.7379 x 10" 8 (P 2 - 5981. 5) 2 (B.12) 

If SP 2 is the uncertainty in P 2 , the the following equation holds: 

X 2 (P 2 + SP 2 ) = x 2 ni „(l + Nd * N ) = xLn x 1-000340, (B.13) 

where N d = 2998 is the number of data points, and N p = 61 the number of fitted param- 
eters. Comparing Eqs. (B.12) and (B.13), we get SP 2 Rj 140. From Eq. (B.8) we then 
derive: 

(S) = 1.80(2). (B.14) 
[The uncertainty in Pi is much smaller than SP 2 , and can safely be ignored.] 

B.4 Magnetic structure (single crystal) 

The first sample we received consisted of numerous small, green, transparent, thin 
(plate-like) crystals, with total mass of about 2.7 g. A few small crystals have been sent 
to a specialized laboratory for X-ray structure determination. Thirty-four largest crystals, 
those with mass m > 3.5 mg, were put aside for single crystal experiments; their total mass 
was msGs = 246.4 mg. The largest single crystal, labeled "A," had mass = 35.2 mg, 
and was used for all TAS neutron scattering experiments. (Masses of a few other largest 
crystals were: m B = 28.5 mg, m c = 6.8 mg, m D = 6.8 mg, and m E = 6.5 mg.) 

In this section we present the details of magnetic structure analysis of single crystal 
neutron scattering measurements; it complements Sec. 5.5 (p. 82). First, we will derive 
the Eq. (5.11) (p. 87), the scattering cross section for the proposed magnetic structure in 
RbFe(Mo0 4 ) 2 . We start from the general expression for elastic magnetic scattering: 8 



/dax 



1 



2' 



2 



) =( 7 r ) 2 -gF(k) e- 2W I, (B.15) 



where 



I = Y1 (<W - kh) ^(5,?)<5f ) exp (ik • (1 - 1')). (B.16) 

a/3 IV 

The definitions and values of constants used in these equations are given in Table of Con- 
stants on p. 184, and of symbols in the List of Symbols on p. 182. The magnetic form 
factor, F(k), is given by Eq. (3.25) (p. 23), and numerical constants for Fe 3+ ion in Table 
3.5 (p. 23). 

8 These are Eqs. (8.20) and (8.21) from Squires [24]. 
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We define ion positions as follows: 9 

1 = 1/ = l mno = ma + nh + oc, (B.17) 

where the unit vectors along the principal crystal axes, a, b, and c, are expressed in their 
Cartesian form as: 



a = ax 

x # 

2 2 - 



b = a 



c = cz. (B.18) 

The helical incommensurate 120° spin configuration is represented with the following spin 
orientation: 

(S m „ ) = (S) (xcos0 mno + y sin0 mno ) 

= (S) Re [(x - iy) exp(i0 mno )] , (B.19) 

where 

(771 ~\~ 71 \ 
— g— + Aoj (B.20) 

is the in-plane rotation angle of spin (mno), and A = 0.453 controls the incommensurate 
helical ordering along the c-axis. 

The momentum transfer, k, is in Cartesian coordinates equal to: 

k = ha* + kb* + lc*, (B.21) 

where a*, b*, and c*, are the unit vectors in reciprocal space: 

a = T< x + l- y) 

b* = 

a 3 
2tt 

c* = — z (B.22) 

c 

First, let us derive some equations useful for carrying out the first summation in (B.16), 
over q, (3 = x. y, z. Absolute value of the wavevector k is equal to: 



k 



hkl — \ K hkl\ 



. I (h\ 2 fVS, 2 V / 3 I \ 2 1 (l\ 2 



9 Index / stands for the triple (mno), for convenience. An equivalent equation holds for primed variables 
and indices. 
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Our measurements were done only in (hhl) plane, and so for h = k Eq. (B.23) simplifies 
to: 

= 21 V(t)' + (;)'■ (B - 24) 

The Cartesian components of the unit vector k/j W = k hhl /k hk i are then: 

ft 



{k h ki) y 



(k hkl ) z = 1 (B.25) 



y/(2h) 2 + (aZ/c) 2 
/ 

^(2hc/a) 2 + / 2 
The equation (B.16) can now be rewritten as: 

a/3 

= ^ ] 4a ~~ ^ ] k a kpl a p 

a a/3 



T T T \ / h 

1 xx 1 xy ± xz \ I n a 



( k x ky lt z ) [ Iy X Iyy Iy Z J I ky J (B.26) 

Izx Izy Izz ) \ k z 

Calculating all nine I a p terms is a nice exercise in algebra and persistence: 

U = E^'>< 5 ?> ex P ( ik • ( ! " 1# )); (B.27) 
one has to use nearly all the equations presented so far in this section, and keep careful track 
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of summations, indices, and delta functions. The result is, however, startlingly simple: 

/= ((S^N/4){(2nf/v ) 



1 
1 




-i 
+ | i 




We note that I xx = I yy , and I, 
into (B.26), and obtain: 



8(ti + h + -)x 5(k' + k + -)x 8(1' + / + A) 
3 3 



+5{ti + h--)x 8(k' + k - -) x 6(1' + I - A) 
3 3 



8(ti + h + \) x S(k' + k + \) x <S(Z' + / + A) 



-<S(ft' + h - -) x <J(jfc' + k--)x 6(1' + I - A) 
3 3 

(B.28) 

— — I yx ; all other terms are zero. Then we insert (B.28) 



Noting that: 



we continue: 



Finally: 



h 2 J — h 2 T —\kk\T I -I- it h \T 

xx y yy lrv x rv y \ 1 xy \ < lr * j v' Xj x\ J -'i 



v x \ ± yx | 



Ixx (2 k x ky^j . 



k — 1 — k x + ky + k z , 



i = i xx (2-k 2 + k 2 z ) 

= Ixx (l + k z) ■ 



{S) 2 N (2tt) j 



v 



(1 + *, 2 ) 



x [6{ti + h + \) x 5{k' + k + \) x 5(1' + / + A) 
+<S(ft' + ft - \) x 8(k' + k-\)x 8(1' + 1-A)], 



(B.29) 
(B.30) 



(B.31) 



(B.32) 
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102.04 


(±8.8656) 


112.99 


4559.7 




1 


7 


503.73 


(±20.026) 


691.11 


37430. 


2 


2 





7525.8 


(±185.42) 


9376.8 


2.4003e+05 


2 


2 


1 


93.065 


(±7.8113) 


37.739 


1001.9 


2 


2 


2 


74.933 


(±6.2075) 


114.02 


3354.5 


2 


2 


3 


2068.6 


(±86.790) 


2015.9 


69101. 


2 


2 


4 


813.01 


(±35.064) 


812.94 


33508. 


2 


2 


5 


206.08 


(±13.205) 


123.96 


6245.0 



Table B.7: Measured and calculated intensities of nuclear Bragg peaks in RbFe(Mo0 4 ) 2 , 
at T = 1.7 K. R factor: 14.7%. Note the effects of anisotropic extinction: the observed 
intensities are relatively weaker close to the basal plane, due to the irregular, plate-like 
shape of the crystal. l calc is calculated using Eq. (5.9). 



and: 



VdQ/el 



[5{ti + h + \) x 8(k' + k + \)x 8(1' + 1 + A) 

+6(ti + h--)x S(k' + k - -) x 6(1' + 1 - A)] . 
3 3 



(B.33) 



Note that this is the scattering cross section for a single domain only. Equation (5.9) 
(p. 86) has to be applied before we can compare the calculated and observed intensities. 
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u 
ll 


K 


( 


J 

-'obs 


XT ) 


J 

i calc 




0.000 


0.000 


2.453 


1.1085 


(±0.44340) 








0.000 


0.000 


2.547 


2.8260 


(±1.1304) 








0.000 


0.000 


3.453 


3.1418 


(±1.2567) 








0.000 


0.000 


3.547 


0.89587 


(±0.35835) 








0.000 


0.000 


4.453 


<0.066679 










0.000 


0.000 


4.547 


<1.5558 










0.000 


0.000 


5.453 


<0.35766 










0.000 


0.000 


5.547 


2.5744 


(±1.0297) 








0.000 


0.000 


7.453 


<0.64281 










0.000 


0.000 
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<0.20063 










0.333 


0.333 


0.453 


715.71 


(±286.28) 


711.17 


1457.6 


0.333 


0.333 


0.547 


689.11 


(±275.64) 


713.36 


1530.6 


0.333 


0.333 


1.453 


426.26 


(±170.50) 


491.01 


1817.0 


0.333 


0.333 


1.547 


387.86 


(±155.15) 


459.93 
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Table B.8: Measured and calculated intensities of magnetic Bragg peaks in RbFe(Mo0 4 ) 2 , 
at T = 1.7 K. i? factor: 18.0%. Omitted are peaks which had a large contamination from 
nuclear peaks. See also Table B.7, and its caption. 
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B.5 Interlayer coupling 

As described in Sec. 5.5 (p. 82), knowing the incommensurate helical magnetic or- 
dering, we tried to determine the ratio of inter- to intralayer exchange constants. In this 
analysis, we will consider the interactions between a spin, and its nearest and next-nearest 
neighbors on one neighbor layer. 10 Each magnetic ion has one nearest neighbor in each 
neighbor layer; let us denote the exchange constant between them with J 2 . There are six 
next-nearest interlayer neighbors, and they come in two non-equivalent locations; let their 
exchange constants be J 3 and J 4 . (See Fig. 5.11.) 

We assume that the spins in each layer reside in a 120° configuration. Spins on a 
neighboring layer are rotated by an angle <fi. The interlayer interaction energy is then: 



between layers 



+ J 2 COS (j) 

+3 J 3 cos(60° + </>) + 3 J 4 cos(60° 
3 



J2 + -(J S + Ji) 



COS (ft 



3x/3 



-Js + Ji) 



sin i 



(B.34) 



Since nearest interlayer spins are nearly antiferromagnetically ordered, we conclude that 
J 2 < 0. When the magnetic order is established, £ b .i. is at a minimum; so we now minimize 
it: 



dE hl . 



dcj) 



0. 



J2 + + Ji) 



3x/3 



(-J3 + Ji) 



sin 1 



cos </) E 



(B.35) 



We get: 



tan 1 



¥(-^ + j 4 ) 

^2 + 1(^3 + ^4) 



(B.36) 



We know that <j) min = A x 360° ~ 163°. Thus, we can use Eq. (B.36) to relate J 2 , J3, and J4. 
Unfortunately, however, this is insufficient to determine the relative strength of interlayer 
interactions, as we were initially planning. We can only conclude that J 3 7^ J 4 ; otherwise, 
the interlayer ordering would be either ferro- or antiferromagnetic, not incommensurate 
helical. 



10 There are, of course, two neighbor layers. It is, however, sufficient to consider the interactions with only 
one. 
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B.6 Arc scan 

Arc scan is a specialized type of scan, which we devised to study the Kosterlitz-Thouless 
transition, described in Sec. 5.6, p. 98. The experiment was carried out using SPINS cold 
neutron spectrometer. The motivation for this scan was to maximize the neutron intensity. 
The spectrometer is used in the two axis mode, by placing the detector before the analyzer. 
In the scan, the I component of the wavevector transfer, Q = (hhl), is undetermined; only 
the h component is known. As the KT transition occurs in two-dimensional systems, we 
only need to know the component of the wavevector in the basal plane, h. We achieve this 
by keeping the crystallographic (00/) axis aligned with the final wavevector, k/: 

foot = c*= k f . (B.37) 

This is achieved by keeping spectrometer angles ^43 and A4 equal: 

A3 = A4. (B.38) 

Scan is then done by using increment buffer (I-buffer) in Instrument Control Program (ICP), 
by scanning angles A3 and A4. With A4 = 29, we can write: 

Q = ki{l - cos 20) (B.39) 

Analysis shows that: 

k- 

h = — sin 29 

TllO 

I = — (1- cos 29), (B.40) 

7"001 

where 

Tno = y/3 a* = 2.2602 A _1 

t 001 = c* = 0.8533 A -1 . (B.41) 

We see that the scan plots a circle, which passes through the origin, in the reciprocal plane. 
Thus the name: arc scan. 

We made measurements in two configurations, with the following fixed energies (i.e. 
incident neutron energies): 

(i) E = 6.68 meV (ki = 1.7957 A" 1 ): 

29 = 24.8072° 

29 = 57.0487° (B.42) 



( 5 J 0.1942); 
(| I 0.9598); 
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(ii) E = 13.7 meV (jfej = 2.5716 A" 1 ): 



^ 0.1322); 2# = 17.0361° 
(| | 0.5715); 20 = 35.8700°. 



(B.43) 



We used G-40'-40'-open collimation with PG-before filter ("G" stands for "guide," i.e. 
the collimation is defined by the neutron guide, not collimators). 
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Appendix C 
QS-ferrite - Details 



C.l Disc-chopper-spectrometer (DCS) energy resolution 

Here, we derive the energy resolution of Disk Chopper Spectrometer (DCS), a time-of- 
flight (TOF) cold neutron spectrometer at NCNR. This spectrometer was used for inelastic 
powder measurements on QS-ferrite (Ba2Sn2Ga 3 ZnCr 7 022), and the analytical expression 
for energy resolution as a function of energy, AE(E), was needed for analyzing the data. 

We used (i) A = 2.3 A neutrons, where incident neutrons have the following wavevec- 
tor, energy, and speed: 





=2.3 A ~ 


27T o -1 

— = 2.73 A , 

A« 




=2.3 A 


h 2 k 2 

1 = 15.5 meV. 
2m 




=2.3 A ~ 


\ — - = 1720 m/s 
V m 



(C.l) 



and (ii) A = 4.8 A neutrons, where: 





o 

! A 


= 1.31 A \ 




o 

! A 


= 3.55 meV 




o 

iA 


= 826 m/s. 



DCS uses a series of eight chopper aluminum disks to monochromate the incoming 
continuous beam of cold neutrons. The resulting beam is not only monochromatic, but also 
pulsed; this time structure is used to determine the neutron energy transfer. The distance 
between the last, ninth disk, the chopper disk (the one that determines the time-width of 
the pulse, At), and the sample, is: 



li = 1.50 m; chopper- sample distance, 



(C3) 
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while the sample to detector distance is: 

// = 4.01 m; sample-detector distance. (C.4) 

In our experiment, the choppers were rotating at: 

/ = 20 000 rpm = 333 rps. (C.5) 

To eliminate contamination from higher order neutrons (such as 2A), we used every other 
pulse; 1 thus: m = 2. The time between pulses was therefore: 

T = ^ = 6.0 ms. (C.6) 

We already defined Vi and E { as the incident speed and energy. Similarly, for scattered 
neutrons we have: 

= vf • 

Next, let ti be the time it takes a neutron to get from the chopper disk to the sample; tf the 
sample-detector time; and t the total chopper-detector time. Then: 

t = ti + t f = - + (C.8) 

Vi v f 



Follows: 



and 



«/ = rir> (C9) 



m l f 



E f = — ? . (CIO) 

E = Ef — Ei, (C.ll) 

E = — f - x-Ei. (C.12) 

2 (t-U/vt? 

The uncertainty in energy, AE, is primarily controlled by uncertainties in (i) the time of 
arrival of neutron in detector, At; and (ii) the incident energy, AEi. We have: 



Inserting this into: 
we get: 



dvi 



(«■) At (be, 



+ 1 



AE?. (C.13) 



lr This is achieved by rotating the chopper disk (the ninth disk) at angular frequency /chopper = f I'm- 
The first eight disks (they monochromate the beam) rotate at angular frequency /. The energy is selected by 
adjusting the rotation phase, not the frequency, of the first eight choppers. 
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The uncertainty in t is equal to: 

At = ^(At z y + (At / )2; (C.14) 

here, Ati is the time- width of the pulse coming from the chopper disk, and Atf is the 
uncertainty in timing of the neutron detection. Analyzing the elastic scattering in our data, 
we found that A^ 79 //s; independent of energy. We should note that the pulse profile 
is not rectangular; instead, it is something between triangular and Gaussian. The value for 
Ati reported is the FWHM for a Gaussian fit. The value of Atf is equal to the width 
of the detector time channel. Those are energy dependent, and for elastic channels are: 

o o 

Atf Ri 5 /is for A = 4.8 A neutrons; and Atf 2 /is for A = 2.3 A. We conclude that 
Ati 3> Afy, and so: 

At « At, « 79 yus. (C.15) 

Since we found the width of the elastic scattering peak to be energy independent (At 
79 /us), we conclude that the uncertainty in incident energy, AEi, is negligible; i.e. the beam 
if very monochromatic. [For elastic scattering, Vi = Vf, and so: 

Avi = Av = l -i-^At. (C.16) 
t 2 

Analyzing the widths of elastic peaks at different energies, we estimate that Avi < 0.1 m/s, 
and AEi < 0.04 meV.] This greatly simplifies Eq. (C.13), as we can ignore the term inside 
the square bracket. After carrying out the derivative, we get: 

AE = SfLlJtll At (c. 17) 

lfy/m 

Numerical evaluation yields: 

AE = 0.0172 x (E + Eifl 2 ; AE, E in meV. (C.18) 
We tested this result against experimentally measured elastic energy widths: 

( A ^xp) £=3 . 55meV = 0.126 meV 
( Alcaic) £=3 . 55meV = 0.116 meV 

( A ^x P )^ 15 . 5meV = 1.17 meV 

( A ^aic) £=15 . 5meV = 1-04 meV, (C.19) 

and found measured widths to generally be about 11% larger than calculated, though the 
functional dependence was correct. The discrepancy may be due to non-Gaussian profile 
of the neutron pulse: We used the standard error analysis, which is strictly valid only for 
Gaussian error distributions. Accordingly, we modified Eq. (C.18): 

AE = 0.0191 x (E + Ei) 3/2 ; AE, E in meV, (C.20) 
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and used this for our analysis. 

For lower incident energies (Ei <C 2 meV), one could instead use: 

AE = \J [0.0191 x (E + Eifl 2 ] 2 + (0.04) 2 ; AE, E in meV. (C.21) 

For lowest energies, however, the term in square bracket in Eq. (C.13), which was negli- 
gible for our energy range, becomes dominant, and the analysis presented thereafter is not 
applicable. Instead, one should use use the full resolution equation: 



AE 



2 S (E + E i ) s o 



l 2 f m 



At 2 + 



l KE i + EY 
P f Ef 



1 1/2 



+ 1 AE- 



(C.22) 



Before using this equation, one should determine the actual intrinsic energy width of the 
incident beam, AEf, the value we reported earlier is merely an upper limit estimate. On 
DCS instrument, however, one never enters a regime where this would be important. 



C.2 Normalization of magnetic scattering intensities 

Introduction 

There are many factors that determine the actual number of neutrons that are detected 
in an experiment. The intensity will generally be proportional to the sample size (for small 
samples), the flux of the incident neutron beam, and detector sensitivity. There are innumer- 
able other factors that are chiefly instrument specific: reflectiveness of monochromator(s), 
analyzer(s), collimators, and neutron guides; container (mainly aluminum or vanadium) 
and air scattering; divergence of the beam; higher order scattering (X/n neutrons for TAS 
instruments, nX for TOF); monitor sensitivity; filter efficiency; etc. If one could model all 
of these, it would be possible, in principle, to directly relate neutron counts to some phys- 
ical property of the sample; the procedure would, however, be highly sensitive, and likely 
unreliable. Instead, measurements are commonly done on two samples; 2 one for which the 
result is known in advance, and which serves for calibration, and then on the sample of 
our interest. The instrument dependent factors are then the same, and the differences in 
intensity are due only to physical properties of the samples. The principal normalization 
techniques are: 

(i) Vanadium normalization. Sometimes, the calibration sample is made of vana- 
dium. Vanadium has a large incoherent elastic, and negligible coherent, scattering. 
Incoherent elastic scattering is nearly perfectly isotropic, nearly proportional to the 

2 The two samples need not really be distinct. For instance, magnetic scattering can be calibrated against 
the nuclear scattering off the same sample. 
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mass of the sample, 3 and very easy to measure and manage analytically. This makes 
vanadium a very good candidate for calibration measurements. Specifically, vana- 
dium measurements are done on multidetector spectrometers; detectors generally 
have different sensitivity, 4 which are determined by a vanadium run. Before ana- 
lyzing data from DCS, which has 913 detectors (or FCS, which has 63 detectors), 
a detector sensitivity correction, based on vanadium run, must be applied. Detector 
sensitivity correction is, however, a relative correction, and is usually not used for 
normalization. 

07) Normalization by incoherent scattering. Rather than using a separate vanadium 
sample, we can utilize the incoherent scattering from the sample of interest. Most 
samples have considerable incoherent elastic scattering, 5 which is isotropic, easy to 
measure, and easy to calculate once the chemical (and isotope) composition is known 
(the structure need not be known). The disadvantage is that one must be certain that 
the incoherent scattering is really coming from the sample, rather from some outside 
sources; this can sometime be difficult to ascertain. 

(Hi) Normalization by nuclear scattering. Normalizing magnetic scattering intensities, 
from neutron counts (or "arb. un.") into physically meaningful units, is commonly 
done by comparison with nuclear scattering intensities. This technique has several 
advantages: (i) nuclear structure is usually better known than the magnetic one, and 
can even be determined by other techniques, such as X-ray scattering; (ii) nuclear 
peaks are generally much more intense than magnetic ones; (in) the intensity of 
nuclear peaks is nearly temperature independent; 6 (iv) contamination in the spectrum 
is usually easier to notice and isolate. Nuclear scattering, therefore, provides a useful 
yardstick against which the magnetic scattering is measured. 

Analytic expression for neutron intensity is (both nuclear and magnetic, coherent and 
incoherent): 



where C is an unknown instrument dependent constant, a spectrometer constant. We need 
to determine this constant in order to normalize our data. It is important for our analysis 
that the resolution function, Rq^{Q — Q',cu — cu 1 ), is normalized to unity: 



3 For moderately small samples. In large samples absorption becomes significant, but for isotropic samples 
it is easy to track analytically. 

4 Typically, these can vary from one detector to another by a factor of 2 or more. 

5 Though usually not as large as that of vanadium. Aluminum is a notable exception having negligible 
incoherent scattering. 

6 Temperature dependence of nuclear scattering is mainly controlled by the Debye-Waller factor, which 
varies very weakly with temperature, and has a simple functional form. 




(C.23) 




(C.24) 
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Scattering cross section for magnetic scattering, is: 

The scattering function, spherically averaged for a powder sample, is: 

S(Q,co)= f^.^(^p-QaQp)S aP {Q,uj). (C.26) 

J 71 a/3 

The dynamic spin correlation function: 

S ^ u ) = jfThf d*e-']T<S£(*))<S£,(0)) exp(-iQ-(R-R')). (C.27) 

^ J RR' 

Normalization by incoherent scattering 

We shall first describe our attempt to normalize the magnetic intensity using incoher- 
ent scattering. The approach presented here follows that given in Ref. [89]. 7 Knowing 
the chemical composition of QS ferrite, Ba2Sn2Ga 3 ZnCr 7 022, we can easily calculate the 
elastic incoherent scattering cross section: 8 

1 / d(T \ calc 
N V df2 / el,incoh 



d 

Because the incoherent scattering is (i) isotropic, (ii) elastic (A; = k'), and (Hi) the resolu- 
tion function is normalized (Eq. (C.24)), equation (C.23) gives: 

(Afj \ calc 
— ) = NC 1.09 x 10- 28 m 2 . (C.30) 

fli l / el,incoh 

We experimentally measure the incoherent elastic scattering by integrating the elastic 
peak over the wavevector windows which contain no nuclear Bragg peaks: 

ijLoh = f-^^j^ 1 ^®'^ = 1.32 x 10- 5 meV (cts/time); (C.31) 

here (cts/time) stands for number of counts per unit of time; an arbitrary unit that will 
cancel out later. Follows: 9 

rcalc 2 

(NC) incoh = = 8.25 x 10- 24 r . (C.32) 

V yincon xp , tg / time 

eljincoh 



7 See also Ref. [90]. 
8 Notation: 



(Si =f^>M- (C28) 



9 For practical calculation, we actually adopted the following constant: NC/r^ = 2.85 x 
10 5 [meV (cts/time)] -1 . 
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Thus determining the spectrometer constant C, we can use it now to relate theoretically 
meaningful values in equations (C.25)-(C.27) (such as (S)) to neutron counts. 10 

Normalization by nuclear scattering 

Let us now describe normalization using nuclear scattering. We follow the approach 
given in Ref. [91]. For a single crystal, nuclear scattering cross section is: 

/daya, (2nf £ _ r)|Fjv(Q)|2 . (C 33) 

and the structure factor, F N (Q), is: 

F N (Q) = b d exp(iQ • d d ) exp(-W d ). (C.34) 

d 

As before, having a powder sample, we spherically average the scattering cross section: 11 

tr = ,M E fc), iM1 (C35) 

Vdf2/coh,el D 47TT 2 

T 

We turn again to Eqs. (C.23) and (C.24). Assuming that the spectrometer resolution (i.e. the 
width of Rqu{Q — Q',uj — uj')) is much smaller than the integration range, 12 we have: 

/ duj dQ I(Q,uj) = cf duj[ Q2 dQ(^-) 

J-e JQ! J e VdQ/cohxl 

= £ iif^n, 1^(0)1'. (C.36) 

Q!<r<Q 2 

Here, subscript "el" refers to nominally elastic scattering, i.e. a band of nearly elastic scat- 
tering, the width of which is resolution limited; the energy integration should be slightly 
wider than the band. For integration, we used Q 1 = 0.70 A and Q 2 = 2 AO A, while for the 
elastic width integration we used e = 1.20 meV. 13 Knowing the structure, we used (C.34) 
to calculate the sum in (C.36): 

— E * ( ?~ 2 T) "t I^(Q)| 2 = 3.18x10- 8 m. (C.37) 

Qi<t<Q 2 



10 Voila! Note that we do not need to know TV, as both the nuclear and the magnetic intensity is proportional 
to TV. Still, TV is easy to calculate knowing the sample mass: TV" = N a m/M. 
n Here, n T is the multiplicity of the nuclear Bragg peak t. 
12 In practice, the instrument resolution is always such that this is the case. 
13 We tried a number of different cutoff values, and arrived at these as the best choices. 



168 



The numerical integration of measured intensity, in (C.36), gives: 

f duj dQ /^ eI (Q, u ) = f duj dQ fc p aLel (Q ; u) - ILZUQ, a,)) 

= 8.53 x 10 4 A" 1 meV (cts/time). (C.38) 

[Note that I(Q,uj) in Eq. (C.36) refers only to nuclear coherent scattering (i.e. Bragg 
peaks); accordingly, we used: 

C P a,e.(^ ") = i:ZM <*) + Coh,e>(^ ( C - 39 ) 

where /f n x C Q h el (Q, u;) is, in practice, simply the elastic intensity between Bragg peaks.] In- 
serting these numerical values into (C.36), we get: 14 

(NC) coh = 3.73 x 1(T 23 meV (cts/time). (C.40) 

Now, comparing results (C.32) and (C.40), we couldn't escape noticing the big dis- 
crepancy between the two numbers. We put our calculation to scrutiny to verify that no 
mistakes we made. Then, we decided to accept the second result, (iVC) CO h, for our analy- 
sis. As mentioned in the introduction, although we try to shield against it, the incoherent 
scattering can originate not only from the sample, but also from outside sources. Thus, 
we cannot be certain that all the intensity we integrated over in (C.31), f Ie\fncoh(Qi ^ s 
indeed the incoherent scattering that comes from the sample. 15 The coherent scattering is, 
on the other hand, well localized in the (Q, E) space, and the sharp, strong Bragg peaks 
are easy to identify. Therefore, the normalization procedure based on coherent scattering 
should be more reliable than that based on incoherent. 



14 We used NC/rl = 1.28 x 10 6 [meV (cts/time)] -1 in our calculation. 

15 For instance, it could come from an impurity in the sample, or elastic background due to incomplete 
shielding. Hydrogen isotope 1 H is among the strongest incoherent scatterers; thus moisture in the sample 
or spectrometer could cause increased incoherent scattering. All these contribution should, however, be 
negligible in practice. 
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Appendix D 
Thesis proposal 



Following pages are a reproduction of an essay submitted as a part of my Graduate 
Board Oral Examination (GBO), taken on April 23, 2001. Since this essay is, essentially, 
considered a thesis proposal, as well as a summary of the research performed up to that 
time, I bring it here unchanged 1 in its entirety, for completeness and collation with the 
thesis. 



lr The formatting, however, had to be slightly modified, in order to conform to the printed format of this 
thesis. 
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I. INTRODUCTION 

Recently, a new classes of materials with novel and unusual properties have been found in seemingly 
unrelated areas. While the mainstream research continued within each area without regard to the other, there 
seems to be a link of new physics between the two. Structurally and magnetically underconstrained systems 
may share deeper connection through the concept of geometrical frustration. The physics behind negative 
thermal expansion could be related to low temperature behavior of geometrically frustrated antiferromagnets. 
The research proposed here attempts to reveal some of these connections. 

II. UNDERCONSTRAINED SYSTEMS 

A. Negative Thermal Expansion (NTE) 

Recently, a large (a s=s 10~ 5 K -1 ) and isotropic negative thermal expansion has been observed in 
ZrWiO%, over a wide temperature range of 0.3 to 1050 K (linear between 50 and 420 K)[l]. The dis- 
covery has motivated further study of the phenomenon, since no other known material exhibits all of these 
properties. 

The discovery of NTE has attracted considerable interest, both for its applicability, and need for theoret- 
ical explanation. Potential applications range from industrial, such as optical mirrors and electronic circuit 
boards, to more practical (cookware and dentistry)[5-16]. By incorporating NTE materials into composites 
with PTE materials, precise control over thermal expansion can be achieved (with zero thermal expansion 
often being most desirable). 

The discovery prompted a need for theoretical explanation. NTE has been known to exist in many 
materials, but was confined to a very limited temperature range, mostly close to a phase transition (a well 
known example being water near freezing temperature), and as such has been considered an anomaly, in 
contrast to positive thermal expansion. 

Soon, discoveries of other materials with large (though anisotropic) NTE followed [2-5], such as Sc2{WOa)z 
[5]. Common to all NTE materials is the structure of corner sharing quasi-rigid units, suggesting existence 
of rigid unit modes (RUM). Rigid units consist of AO e octahedra (A = Zr, Sc, Y, Lu), and BO4 tetrahedra 
(B = W, Mo), where strong metal-oxygen bonds are necessary to maintain rigidity of these units, which 
are coupled only by their corners and thus able to (relatively) freely rotate. Theory was proposed in which 
collective rotations of these rigid units is responsible for NTE [6,25]. 
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B. NTE Model 

While RUM are essential part of the model, their existence does not necessarily imply NTE. Indeed, rigid 
unit modes have been found to exist in SiO^ and ^4/04 [7], neither of which exhibits NTE. 

Instead of enforcing rigid units in the crystal model, our approach has been to use the interatomic po- 
tentials to study crystal dynamics. Interatomic potentials can be either (a) calculated from the first principles 
[8], or (b) phenomenologically derived by fitting to various experimental data. Since these crystals are very 
complex, and have large unit cells, direct calculation is very difficult. 

Interaction between atoms can be described by various potentials that can be grouped according to (a) 
number of atoms involved (two, three, or four body), (b) functional dependence (harmonic, exponential, 
etc.). Atoms can be charged, and interact through long range Coulomb potential. Further, each atom can be 
split into massless electron shell, and nucleus that contains all the mass of the atom, thus modeling atomic 
polarization [13]. 

Consistent fitting of interatomic parameters to (a) crystal structure, (b) experimental density-of-states, 
and (c) experimentally mapped phonon dispersion relations, yields detailed description of various crystal 
properties. Although the fit has not yet converged, we anticipate that once completed it will be able to 
elucidate NTE in the material of interest. 

While it appears that this method may be applicable to other complex systems as well, and provide 
detailed description of various crystal properties [9] (along with NTE), the results obtained are difficult to 
generalize. Though it may reproduce thermal expansion in various materials, the method makes it difficult to 
qualitatively answer why only certain materials have NTE. 

C. NTE Generalization and Its Difficulties 

Rigid unit modes, thought to be responsible for NTE, are made possible by crystal structure consisting 
of an underconstrained network of corner sharing polyhedra. However, a number of materials that are also 
well described by rigid units do not exhibit NTE over wide temperature range [7]. RUM in these materials 
develop soft modes that drive displacive phase transitions [13], an example being BaTiOz which undergoes 
a cubic to tetragonal phase transition at 400 K [1 1,21,23]. Such phase transitions are characterized by a large 
contraction of volume in a narrow temperature range around the transition temperature. 

It has been proposed that NTE is driven by the same mechanism as those phase transitions; however, 
in NTE materials crystal symmetry inhibits the phase transition, thus resulting in qualitatively new behavior 
[11]. While phonon dynamics try to drive the system into a favorable lower energy (and volume) state, it is 
opposed by structural symmetry that is not compatible with such a transition. 

Although this approach seems to be qualitatively correct, it is difficult to quantify the incompatibility 
between symmetry and phonon dynamics in a simple manner (symmetry being un-tunable variable). Without 
a quantitative measure of this competition it is difficult to evaluate and generalize this approach to all rigid 
unit systems. 

In order to grasp only the most essential features of such competition between incompatible tendencies, 
let us analyze other systems seemingly unrelated to NTE, in which similar competition exists. 

D. Geometrical Frustration 

Geometrical frustration is a generic term for all systems where geometry of the system (global symmetry) 
is incompatible with local interactions between individual system members, therefore preventing all local 
interaction energies from being simultaneously minimized. 

Frustration is attainable (a) through site disorder, as is common in glass systems, or (b) through frustration 
due to intrinsic lattice symmetry, as in the case of geometrical frustration. Geometrically frustrated systems 
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have advantage of having exactly known Hamiltonian. Also, by removing disorder the effects of frustration 
can be better singled out. 

Some striking features of frustrated systems are thermodynamically large number of ground states (i.e. 
number of states grows with number of atoms in the system), and finite entropy at T=0, indicating significant 
underconstraint of these systems. Particular realizations of frustrated systems have unique signatures in 
density of states, susceptibility, specific heat, etc., that set them apart from other systems. 

E. Geometrically Frustrated Magnets (GFM) 

One realization of frustrated systems that is particularly suitable for study are Geometrically Frustrated 
(Antiferro) Magnets (GFM). Compared to other systems, its main advantages are: 

(i) Choice of size and dimensionality of degrees of freedom (spins), 

(ii) Various types of spin-spin interactions, 

(iii) Suitability for experimental study. 

All known geometrically frustrated magnets come as networks of either triangles (2D) or tetrahedra (3D). 
Various arrangements are possible, such as corner sharing (2D Kagome lattice in SrCrgGa^iQ [17], 3D 
spinel in ZnCr 2 4 [18]), or edge sharing (2D triangular lattice in LiCr0 2 [19], 3D FCC in K 2 IrCl 6 [20]). 

One (dimensionless) measure of frustration in magnetic systems is the ratio of Curie-Weiss temperature, 
9cw, which indicates the magnitude of spin-spin interaction, and Neel temperature, Tq, at which the phase 
transition driven by these interactions occurs. Mean-field theory for two-sublattice antiferromagnet (not 
frustrated system) yields the ratio of one [27]. However, frustration tends to suppress the critical temperature, 
and in GFM the ratios of ten to few hundred have been found. 

The low and medium energy regions of GFM can display novel behavior, such as spin-glass and spin- 
liquid effects. Experimentally, they seem promising search areas for new phases, and possible new super- 
conducting mechanisms. Theoretical understanding of these systems is very challenging both analytically 
and computationally, compared to conventional magnets. While interesting in their own right, they may also 
provide better understanding of other frustrated systems, such as those with structural frustration. 

F. Generalization 

To connect frustration in different systems, namely magnetic and structural, we should find thermo- 
dynamically equivalent variables in two systems. Since NTE is characterized by a = d\nV/3dT = 
(dS/dV)T /B < (entropy increases with decreasing volume), magnetic analogies would be (dM/dT)n > 
and (dC I8H)t > 0, both of which are unique properties of antiferromagnetism [11]. Frustrated soft mode 
would then correspond to a state in geometrically frustrated antiferromagnet, thus further extending analogue. 

Since magnetic Hamiltonians are easier to express and coupling constants easier to measure, they offer 
much better ground for studying various effects of frustrations. It is to be expected that many features found 
in magnetic systems will not have its analogue in structural ones. Therefore, study of magnetic systems 
should not be limited only to properties already known in, or applicable to, other systems. 

II. RESEARCH TECHNIQUES 
A. Neutron Scattering 

The power of neutron scattering technique as a probe of condensed matter is based on the following 
properties: (a) energy and wavelength of thermal neutrons are comparable to those of most interesting features 
in solids, (b) neutrons have no net charge and probe bulk of a solid, and (c) neutron has a magnetic moment, 
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which can probe local magnetic fields within the solid. No other probe (photon and electron being the most 
common ones) has such a unique and desirable combination of properties. 

Elastic neutron scattering can be used to probe crystal structure, static properties, phase transitions, etc., 
while inelastic scattering can map dispersion relations, density of states, excitation spectra, etc. Experiments 
can be performed in various environments, variables being temperature, magnetic field, pressure, etc. Un- 
like thermodynamic measurements, which measure certain macroscopic variable, neutron scattering provides 
direct access to microscopic structure and dynamics. 

B. Thermodynamic Measurements 

Measurements of macroscopic thermodynamic properties, such as susceptibility or heat capacity, give 
us important information about the state of the system. Susceptibility can provide us with a magnitude of 
spin interactions (through Curie-Weiss temperature), and temperature of phase transition, and thus give us 
a measure of the frustration present in a system. Heat capacity data can be related to density of states and 
entropy of the system. The nature of many systems can be immediately recognized through their macroscopic 
signatures. 

C. Computational Techniques 

In order to gain a detailed insight into system dynamics, it may be necessary to use various computational 
techniques. Modeling of crystal dynamics based on phenomenological interatomic potentials has been shown 
as a useful tool in understanding microscopic behavior of complex systems. Other computational methods 
for simulation of spin systems may be necessary to enhance understanding of frustration. 

D. Other Techniques 

Various other techniques may be (and some were already) necessary in the course of research progress, 
such as structure determination, crystal growth, Raman spectroscopy, etc. 

IV. RESEARCH PLAN 
A. NTE Experiments 

A number of NTE experiments has already been performed at neutron scattering facilities at NIST [12]. 
These include temperature dependence of lattice parameters, density of states measurements, and mapping 
of phonon dispersion relations. Various environments were used, variables being temperature and pressure. 
Mapping of phonon dispersion relations (using inelastic neutron scattering on single crystals) constitutes a 
critical part of the research, as they are necessary for unique determination of interatomic parameters. Since 
large single crystals were unavailable for ZrWiO%, our focus shifted to Scz^WOijz, another NTE material 
for which experiments were easier to carry out, though modeling was complicated by larger unit cell and 
lower crystal symmetry. Further experiments may be needed if fitting to current experimental data does not 
converge or does not give unique interatomic parameters. 

In addition to obtaining high quality data by neutron scattering, achieving mastery of, and fluency in, 
neutron scattering techniques is considered of utmost importance. Thorough understanding of neutron instru- 
mentation and design, use of various environment devices (such as cryostats, pressure cells, and magnets), 
and of trends in experimental neutron research, is deemed extremely beneficial. 
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B. NTE Modeling 

The process of fitting interatomic parameters to experimental data seems to be nearing its convergence. 
The attempts were marked by many trials and errors, while searching for the best way to use the experimental 
data. As the method is judged to be applicable to other complex systems, certain conclusions are to be reached 
regarding the most optimal process of extracting parameters from experimental data, and regarding the most 
desirable type of experimental data for such an analysis. 

As an important side product, a program has been made that can calculate, by using interatomic po- 
tentials, neutron intensities as would be observed in an actual experiment (with certain approximations) for 
particular instrument configurations. It appears that a general program for such a calculation has not been 
made before. Therefore it may be useful to continue work on it even after desired research results are ob- 
tained, and convert it into a distributive and user-friendly format. 

C. GFM Experiments 

Some GFM experiments were already performed, both neutron scattering, and susceptibility, though it is 
perceived that bulk of experiments is yet to be done. As the focus on various GFM materials rapidly changes 
in search of the best system for proposed studies, it is impossible to accurately predict which material is 
most promising for future research. Current research includes Ba2Sn<iGa?,ZnCr70<i2, a two-dimensional 
spinel based Kagome antiferromagnet, and RbFe{MoOi)2, a two dimensional Heisenberg triangular anti- 
ferromagnet (on which susceptibility measurements were recently performed). Future experiments should be 
tailored to discovery of new and unknown properties of GFM, and closely coordinated with current research 
developments. 

D. GFM Theory and Modeling 

It may be necessary to apply computational methods for analyzing experimental data, or understanding 
some of the observed properties. Since the current theoretical description of GFM is incomplete and rapidly 
evolving it will be important to follow the literature in this area closely. 

E. Generalization 

Generalization of GFM results to other systems, namely systems with structural frustration, and NTE in 
particular, is considered desirable. However, it is possible that such a generalization shall not be attainable, 
due to the significantly different nature of the systems. Even such a result would be significant, as it may lead 
to resolution of some uncertainties currently present in the field. 

V. SUMMARY 

Geometrically frustrated systems are characterized by the incompatibility between global symmetry and 
local interactions. Competition between two opposing tendencies is followed by underconstraint in dynamic 
system, resulting in novel and unusual properties. By studying two distinct classes of systems, structurally and 
magnetically frustrated, comprehension of the most essential features of geometrical frustration is expected. 
Use of the experimental technique of neutron scattering is considered indispensable. 
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List of Symbols 

B Magnetic field. (Also called magnetic induction.) Unit: Tesla T. 
b Coherent scattering length (except on p. 17). 
frincoh Incoherent scattering length. 

d Index for an atom in the unit cell; used in summations, 
d Fractional coordinate of atom d in the unit cell. 

E Energy (excitation, transfer). 
Ei, Ef Initial, final neutron energy, 
e Elementary charge. 

e, e ds Polarization tensor, vector, for atom d and phonon mode s. 

F M Magnetic form factor. 
F N Nuclear structure factor. 

g Lande gyromagnetic ratio. 

H Magnetic field due to free currents. Unit: Ampere/meter [A/m]. 
h Reduced Planck constant. 

J Total angular momentum operator. 

J Eigenvalue of the total angular momentum operator. 

J Exchange constant. 

j n (kr) Spherical Bessel function of order n. 
J n (k) See eq. (3.22), p. 23. 



k-B Boltzmann constant. 

kf, kf Final neutron wavevector. 
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kj, ki Incident neutron wavevector. 

L Orbital angular momentum operator. 

L Eigenvalue of the orbital angular momentum operator. 

1, 1' Vector of the crystal lattice. 

/, /' Index of the unit cell located at 1, 1'; or index of an atom in the magnetic lattice. 

M Magnetization (magnetic moment per unit volume). Unit: Ampere/meter [A/m]. 
m d Mass of the atom indexed with d. 
m e Mass of the electron. 
m n Mass of the neutron. 
m p Mass of the proton. 

N Number of unit cells in a crystal. 

(n) Bose-Einstein boson occupation number. 

S Spin operator. 
S Spin 

s Excitation (phonon, magnon) index. 

(Sf) An a component of the average dipole magnetic moment on the atom / (magnetic 
lattice). 

v Volume of the unit cell. 

a,P = x,y,z. 

a Total scattering cross section. 

(^), (df^fe) Differential scattering cross section. 

8(lo) Step function. 

i[) Wave function (electron). 

t Vector in the reciprocal lattice. 

X Magnetic susceptibility. 

dfl Element of the solid angle. 
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Table of Constants 



a = 0.529 177 2083(19) x 10~ 10 m, Bohr radius (47re ft 2 /m e e 2 ). 

c = 299792458 m/s, speed of light in vacuum. 

e = 1.609 176 462(63) x 10" 34 C, elementary charge. 

h = 6.626 068 76(52) x 10" 34 Js, Planck constant. 

h = 1.054 571 596(82) x 10~ 34 Js, reduced Planck constant. 

k B = 1.380 6503(24) x 10" 23 J/K, Boltzmann constant. 

m d = 3.343 583 09(26) x 10" 27 kg, deuteron mass. 

m e = 9.109 381 88(72) x 10~ 31 kg, electron mass. 

m n = 1.674 927 16(13) x 10~ 27 kg, neutron mass. 

m p = 1.672 621 58(13) x 10" 27 kg, proton mass. 

N A = 6.022 141 99(47) x 10 27 , Avogadro constant. 

R = 8.314 472(15) J/mol K, molar gas constant. 

u = 1.660 538 73(13) x 10 -27 kg, atomic mass constant. 

7 = 1.913 042 72(45), neutron ^-factor (fjiJfj, N ). 

e = 8.854 187 81710~ 12 F/m, electric constant (l/> c 2 ). 

fj, = 4n x 10~ 7 N/A 2 , magnetic constant. 

/i B = 927.400 899(37) x 10" 26 J/T, Bohr magneton (eh/2m e ). 

A t n = 5.050 783 17(20) x 10" 27 J/T, nuclear magneton (eh/2m p ). 

$o = 2.067 833 636(81) x 10~ 15 Wb, magnetic flux quantum (h/2e). 
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(Source: [16] and [102].) 

Neutron energy equivalents: 



E 


= 1.609 x 1(T 22 J x £[meV] 


A 


= 9.045 A x E- 1/2 [meV] 


k 


= 0.6947 A -1 x £' 1 / 2 [meV] 


V 


= 437.4 m/s x E l/2 [me V] 


T 


= 11.604 K x E [me V] 


f 


= 0.2418 THz x E[meV] 




= 1.519 rad/s x E[meV] 


v/c 


= 8.0655 1/cm x E[meY] 
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Glossary 

120° A magnetic configuration on a triangular lattice. Represents a ground state for clas- 
sical Heisenberg and X-Y spins. See Sec. 5.5 (p. 84). 

\/Z x \/2, A magnetic configuration on a kagome lattice. Represents a ground state for 
classical Heisenberg and X-Y spins (this ground state is not unique; see q = 0). 

Arb. un. Arbitrary units. Commonly used on y axis before the data is normalized into 
absolute, physically meaningful units. 

BNL Brookhaven National Laboratory. See NSLS. (http://www.bnl.gov) 

BT Beam tube. Designation for guides of thermal neutron (confinement building) at 
NCNR. See NG. 

BT-1 Neutron diffractometer at NCNR. 

BT-2, BT-4, BT-9 Thermal triple axis neutron spectrometers at NCNR. 
CDF Composite degree of freedom. 

DCS Disc chopper spectrometer. A TOF (q.v.) spectrometer at NCNR. 
DFT Density functional theory. 
DOS Density of states. 

FCS Fermi chopper spectrometer. A TOF (q.v.) spectrometer at NCNR. 
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Frustration A topology in which the system cannot simultaneously satisfy all the local 
interactions, even in the ground state. Often leads to highly degenerate ground state, 
especially if combined with weak connectivity (q.v.). 

FWHM Full Width at Half Maximum. The full width of a Gaussian peak, at one-half of 
its peak value. 

Geometrical frustration Frustration without structural disorder. Can often be attributed 
to incompatibility of local interactions and global symmetry. 

GULP General Utility Lattice Program. A computer program for the symmetry adapted 
simulation of solids. See Ref. [46]. 

HFBS High flux backscattering spectrometer. A neutron spectrometer at NCNR, with 
extremely high energy resolution. 

HWHM Half Width at Half Maximum. One-half of FWHM (q.v.). 

Kagome A two-dimensional lattice of corner sharing triangles. 

MACS Multi Axis Crystal Spectrometer at NCNR (NG-0). 

NCNR NIST Center for Neutron Research. A part of NIST (q.v.), focused on provision of 
neutron measurement capabilities to researchers. All neutron scattering experiments 
described in thesis have been conducted using NCNR facilities, (http://www.ncnr.nist.gov) 

NG Neutron guide. Designation for guides of cold neutron at NCNR. See BT. 

NDOS "Neutron weighted" density of states. 

NIST National Institute of Standards and Technology. See NCNR. 
(http://www.nist.gov) 

NSLS National Synchrotron Light Source. A part of BNL (q.v.), providing researchers 
with X-ray and UV scattering facilities. Synchrotron X-ray scattering experiments 
described in this thesis have been conducted at NSLS. 
(http://www.nsls.bnl.gov) 
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NTE Negative thermal expansion. 
PTE Positive thermal expansion. 

q = A magnetic configuration on a kagome lattice. Represents a ground state for 
classical Heisenberg and X-Y spins (this ground state is not unique; see \/2, x \/3). 

Ramping (magnetic field) Increasing or decreasing (magnetic field). 

RUM Rigid unit modes. A type of phonon mode that consists of rotation of metal-oxygen 
polyhedra, which remain nearly undistorted. Characterized by low energies. 

Shell In shell atomic model, a massless, charged electronic cloud, connected to massive, 
charged nucleus, called core, via spring potential. 

SPINS Spin polarized neutron spectrometer. A triple axis neutron spectrometer at NIST. 

Spin vortex A type of excitation in magnetic systems, specifically in two dimensions. See 
p. 99. 

TAS Triple axis spectrometer. A type of neutron spectrometer, which uses monochroma- 
tor and analyzer for neutron energy determination. 

TOF Time of flight. A type of spectrometer in which the neutron beam is pulsed, and this 
time structure is used for energy determination. 

Topological charge See spin vortex. 

TRASH A measurement mode of the BT-4 spectrometer at NCNR (before the upgrade). 
Suitable for DOS measurements. 

Weak connectivity A topology in which the order in one part of the system does not 
constrain the ordering in other parts. See frustration. 
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Index 



Analyzer, 40 
Arc scan, 138, 215 

Backscattering, 49 
Bose-Einstein distribution, 30 
Bragg law, 38 
BT-2, 46 

Cadmium, 26, 39 
Collimator, 39 
Cross section, 29 

DCS, 48, 49 
Debye-Waller factor, 29 
Detector, 41 

FANS, 46 
FCS, 45, 48 
Filter, 43 

Fluctuation dissipation theorem, 157 
Focusing 

Monochromator or analyzer, 50 
Form factor, magnetic, 33 

Gadolinium, 26 

General Utility Lattice Program (GULP), 
82 



HFBS, 47, 48 

Interatomic potentials, 82, 177 

Magnetic form factor, 33, 35 
Maxwellian distribution, 23 
Monitor, 42 
Monochromator, 38 

NCNR, 23, 46-48 
Neutron 

discovery, 19 

photograph, 70 

properties, 20 

reactors, 23 

scattering length, 26, 28 

sources, 22 

spallation, 23 
Neutron scattering 

history, 19 

theory, 25 
Normalization, 223 

incoherent, 224 

nuclear, 225 

vanadium, 224 
NSLS (BNL), 16 



Helium atom scattering, 22 



Openness, structure, 9, 66 
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Phonon, 30 

Position sensitive detector (PSD), 44 
Proposal, thesis, 230 
Pyrolytic Graphite (PG), 38 

QS-ferrite (Ba 2 Sn2Ga 3 ZnCr 7 022), 14 

Resolution function, 41, 44, 226 
Rubidium iron molybdate (RbFe(Mo0 4 ) 2 ), 
14, 140 

Rubidium iron molybdate (RbFe(Mo0 4 ) 2 ), 
96 

Scandium tungstate (Sc 2 (W0 4 ) 3 ), 14, 55 
Scattering, 27 

coherent, 26 

cross section, 26 

elastic coherent, 29 

elastic incoherent, 30 

experiment, 37 

incoherent, 26 

inelastic coherent, 30 

inelastic incoherent, 3 1 

length, 26 
Scattering triangle, 41 
Spin vortex, 135 
SPINS, 47, 48 
SQUID, 17 
Structure factor 

nuclear, 29 
Synchrotron, 16 

Time of flight (TOF), 44, 45 



Topological charge, 135 

Triple axis spectrometer (TAS), 37 

Two axis spectrometer, 40 

Vanadium 

normalization, 224 

X-rays, 15, 21 

synchrotron, 16 
tube source, 16 

Zinc chromite (ZnCr 2 4 ), 11, 14 
Zirconium tungstate (ZrW 2 8 ), 14, 55 
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